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I. SUMMARY AND INTRODUCTION 


This report describes a computer program HADY-I for calculating the linear 
incompressible or compressible stability characteristics of the laminar boundary 
layer on swept and tapered wings. The stability analysis and computational pro- 
cedures used in this program are outlined in details in references 1 and 2. 

The eigenvalue problem and its adjoint arising from the linearized distur- 
bance equations with the appropriate boundary conditions are solved numercially 
using a combination of Newton-Raphson iterative scheme and a variable step size 
integrator based on the Runge-kutta-Fehl burg fifth-order formulas. The inte- 
grator is used in conjunction with a modified Gram-Schmidt orthonormalization 
procedure. More details concerning the integrator is found in reference 3. 

The computer program HADY-I calculates the growth rates of crossflow (CF) 
or streamwise Tol Imien-Schl ichting (TS) instabilities. It also calculates the 
group velocities of these disturbances. The program incorporates all methods 
of calculation outlined in reference 1, namely, they are: 

(1) MMSGR, method of maximum spatial growth rates 

(2) MFCWL, method of fixed spanwise component of wavelength 

(3) MFWL, method of fixed wavelength 

HADY-I is restricted to parallel stability calculations, where the boundary 
layer (meanflow) is assumed to be parallel. The meanflow solution is an input 
to the program, see Appendix I. A nonparallel stability computer program HADY-II 
is another part that is under preparation. 

Figure 1 shows a definition of the coordinate system used in this analysis, 

X is the chordwise direction, z is the spanwise direction, and y is the normal 
to the x-z plane. 




II. PROGRAM STRUCTURE 


A. FLOW CHART 



Figure 2.- Overall flow chart for HADY-I program. 





















B. SUBROUTINES, PURPOSE AND OUTPUT 


See Section C for the definition of the Fortran variable names. 

SUBROUTINE FULLP 

PURPOSE = To read from TAPEIO, the compressible parallel meanflow profiles, and 
modify these profiles if incompressible stability calculations are needed. 

OUTPUT = The variables NS, XC, MACH, R, CHL, UFS, ETA, PRA, ND, IE, and the 
arrays YI, U, UP, UPP, W, WP, WPP, T, TP, TPP, MU, MUP, MUPP, ALFA, ALFAP, PR 

SUBROUTINE CFA 

PURPOSE = To calculate the angle of the velocity vector with respect to x at 
the y-location of an inflection point of the crossflow velocity profile 

OUTPUT = EPSICF 

SUBROUTINE CDMINV 

PU R POSE = To calculate the transpose of the eigenvector matrix 
OUTPUT = The transpose matrix A 
SUBROUTINE JOCK 

PURPOSE = To integrate the disturbance equations and to orthonormal ize the 
soTution whenever needed 

OUTPUT = I FLAG, the array PH 

SUBROUTINE FMAT 

PURPOSE = Provides the disturbance equations to be integrated in JOCK 
SUBROUTINE RNITN 

PURPOSE = To iterate on the complex eigenvalues KX. Used only when ICASE = 1 or 
ICASE = 3 

OUTPUT = Updated value of KX 
SUBROUTINE CFITN 

PURPOSE = To iterate on the real eigenvalues ALFAR and EPSI. Used only when 
ICASE = 4 or ICASE = 5 

OUTPUT = Updated values of ALFAR and EPSI 
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SUBROUTINE GIMIR 


PURPOSE = Iteration to maximize the disturbance growth rate, at each chordwise 
location. The iteration is performed on BETAR when ICASE = 3, and on WN when 
ICASE = 4 

OUTPUT = Updated values of BETAR or WN 
SUBROUTINE SONG 

PURPOSE = Numerical itegration of a complex quantity from y = 0 to y = y (the 
edge of the boundary layer) 

OUTPUT = The value of the integral at yg 

SUBROUTINE PROF 

PURPOSE = During the integration of the disturbance equations, meanflow quan- 
tities are needed at a specific y-location. This subroutine uses the real 
function INTER to interpolate the meanflow profiles 

OUTPUT = Meanflow values at the required y-location thev are U. UP. UPP. W. 
WP, WPP, T, TP, TPP, MU, MUP, MUPP, ALFA, ALFAP, PRANDL 

SUBROUTINE GRVEL 

PURPOSE = To calculate the ratio of the components of the group velocity and 
the group velocity angle 

OUTPUT = VIGl, VIG2, TATA, EPGR 
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C. FORTRAN VARIABLE NAMES 


Variable Name 

Description 

ALFA 

d(MU)/dT, where MU and T are the viscosity and temperature 
of the meanflov/ 

ALFAP 


ALFAR 

wave number component in x-direction 

ALFA I 

growth rate component in x-direction 

BETAR 

wave number component in z-direction 

BETA I 

B^-, growth rate component in z-direction 

CHL 

1 /2 

L* dimensional characteristic length = (U^* x*/Ug*) ' 

ETA 

(y*/L*) 

EPS I 

’ ■ (8r/«r) 

EPS IP 

(6,/«,) 

EPGR 

= tan~^ 1’ group velocity 

gr V P ^ y 

angle 

EPSICF 

wave angle for the crossflow instabilities 

IFLAG 

Flag from the integrator to indicate the status of 
the solution 

IE 

Number of points in the boundary layer where the 
meanflow quantities is given 

KX 

= (a^, a^-). complex wave number in x-direction 

KZ 

k^ = (b^, 6^-)» complex wave number in z-direction 

MACH 

M„, freestream Mach number 

MU 

Array of meanflow viscosity 

MUP 

d(MU)/dy 

MUPP 

d^(MUi/dy^ 

NO 

tan <)) = 
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NS 

Station number along the chord 

PH 

Array of solution vector from JOCK 

PR 

Array of Prandtl number 

PRA 

The boundary layer edge value of Prandtl number 

R 

Reynolds number = U^*L*/Ug* 

T 

Array of meanflow temperature 

TATA 

group velocity ratio 

TP 

dT/dy 

TPP 

d^T/dy^ 

U 

Array of meanflow velocity in x-direction 

UFS 

Dimensional freestream velocity 

UP 

dU/dy 

UPP 

d^U/dy2 

W 

Array of meanflow velocity in z-direction 

WP 

dW/dy 

WPP 

d^W/dy^ 

WN 

wave number k or wavelength x/c 

XC 

X location as percentage of the chord 

YI 

Array of y-location in the boundary layer 

yiGi 

Proportional to group velocity in XT.direction 

VIGZ 

Proportional to ujg, group velocity in z-direction 

PHI 

(j), potential flow angle in degrees 

OMEGA 

Complex circular frequency 

FREQ 

Nondimensional frequency F (if f is inputl, or 
dimensional frequency f Ctf F is input) 

CR 

Have speed 


8 



PGR 

D 

DKX 

IWl 

IW2 

IFL 

DEPSI 

DALFAI 

N IT 


Parallel growth rate calculated as 
-a. tan 

The wall value of the eigenfunction which we iterate on 
after convergence Cit should be zero) 

Incremental variation of kx during the iteration 

Number of performed orthonormal tzati on 

Maximum number of allowed orthonormalization 

IFLAG 

Incremental variation of EPS I during the iteration 
Incremental variation of ALFAI during the iteration 
Number of iteration 
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III. COMPUTER PROGRAM USAGE 


The program calculates the stability characteristics at one or more 
chordwise locations. 

A. Machine Requirements 

HADY-I executes on a computer CDC CYBER 175 

B. Storage Allocation 

The program executes in 125000 octal words of central memory 

C. Timing 

Timing for the job depends on several different options available for 
running HADY-I. In a single execution, the CPU required depends on how near 
is the guessed eigenvalue to the exact value, number of chordwise locations 
where output is needed, . . . etc. . . . The CPU requirement for one chordwise 
location using a reasonable guess (not very far from the exact value) for one 
disturbance frequency and one wavelength is about 30 seconds. 

D. Input/Output Files 

The program card is 

PROGRAM HADY (INPUT, OUTPUT, TAPE6 = INPUT, TAPE6 = OUTPUT, TAPEIO) 

TAPEIO is an input file that contains all meanflow profiles at different 
chordwise locations on the wing, starting from the leading edge to the trailing 
edge. At each location, meanflow quantities are stored from y = yg to y = 0. 

E. Control Cards 

The following control cards can be used to execute the program 

JOBN, Tt, CM 

USER, USERNO, PASSWRD. 

CHARGE, CHARGNO, LRC. 

GET, HADY I. 

GET, TAPEIO = TAPEN. 

REWIND, TAPEIO. 

GET, NBLHLIB/UN = 35781 IN. 

LDSET, LIB = NBLHLIB, PRESET = INDEF. 
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HADYI 

EXIT 

7/8/9 end of record 
Input Cards 
e/7/8/9 

NBLHLIB is a special file that contains the integrator 
F. Program Input 

All input variables are read inside the program in a NAMELIST 
statement 

NAMELIST/SHEE/ISPTM, INCOMP, NSTATN, NSEND, NSTEP, ICASE, ICIZER, 

FR, IFR, NFR, DFR, EPSI, lEPSI, NEPSI, DEPSI, EPSIP, lEPSIP, NEPSIP, 
DEPSIP, ALFAR, ALFAI, BETAR, IBETAR, BETAI, WN, IWN, NWN, OWN, CHORD, 
IPRINT, RE, AE, EPS, ITR, OPEGAR, OMEGAI 

A description of these variables follows 

ISPTM Flag for spatial or temporal stability calculations 

ISPTM = 1 spatial 
ISPTM = 2 temporal 

INCOMP Flag for incompressible or compressible stability calculations 
INCOMP = 1 incompressible 
INCOMP = 2 compressible 

NSTATN Chordwise station number where stability calculations start 

NSEND Chordwise station number where stability calculations end. 

NSEND = NSTATN if calculations is needed at one station 

NSTEP Stability calculations are performed every NSTEP chordwise 

stations 

ICASE Flag for the method of calculation 

ICASE =1 TS calculations with MMSGR or MFWL 
ICASE = 3 TS calculations with MFCWL 
ICASE = 4 CF calculations with MMSGR or MFWL 

ICASE = 5 CF calculations with MFCWL 

ICIZER Flag for the method of calculation 

ICIZER = 1 iteration on the eigenvalue for maximum growth 
rate 

ICIZER = 2 no iteration to maximize the growth rate 

FR The disturbance frequency, dimensional (f) in CPS, or non- 

dimensional (F = 2irfu/U^ ) 

IFR Flag for the input frequency FR 

IFR =1 F is an input 
IFR =2 f is an input 
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NFR 

DFR 

EPSI 

lEPSI 

NEPSI 

DEPSI 

EPSIP 

lEPSIP 

NEPSIP 

DEPSIP 

ALFAR 

ALFAI 

BETAR 

TBETAR 

BETA I 
WN 


Number of frequencies to be used in calculations at each 
chordwise station 

Incremental variation in FR if more than one frequency 
is used 

FR (new) = FR (old) + DFR 
The wave angle v in degrees 

Flag to indicate how t is calculated, lEPSI = 2 and lEPSI 
are used for CF calculations 
lEPSI =1 y is an input 
lEPSI = 2 y = <|> +90° 

lEPSI =3 y is calculated from subroutine CFA 
lEPSI = 4 y = (j. 

Number of wave angles to be used in calculations at each 
chordwise station 

Incremental variation in y, if more than one wave length 
analyis needed 

EPSI (new) = EPSI (old) + DEPSI 

Angle of the growth rate vector, y = tan"^ (g^/a-j) 

Flag to indicate how y is calculated 
lEPSIP =1 y is an input 
lEPSIP = 2 y = <(> 

Number of angles y to be used in calculations at each 
chordwise station 

Incremental variation in y, if more than one angle is 
needed 

EPSIP (new) = EPSIP (old) + DEPSIP 

a^, wavenumber component in chordwise direction 

-a., growth rate component in chordwise direction 

wavenumber component, or x^*/C*, wavelength component 
in the spanwise direction. C is the chord length 

Flag for the input BETAR 
IBETAR =1 gp is an input 
IBETAR = 2 is an input 

-g,-, growth rate component in spanwise direction 

k, wave number or x*/C* wavelength 



IWN Flag for the input WN 

IWN =1 k is an input 
IWN = 2 \*/C* is an input 

NWN Number of wavenumbers or wavelengths to be used in the 

calculations at each chordwise station 

DWN Incremental variation in WN if more than one is needed 

WN (new) = WN (old) + DWN 

CHORD The dimensional chord length normal to the leading 

edge 

IPRINT Flag to control output printing 

IPRINT = 1 long print 

IPRINT = 2 short print 

RE Relative error tolerance used by the integrator 

AE Absolute error tolerance used by the integrator 

EPS Relative error tolerance used in calculating the 

eigenvalue 

ITR Maximum number of iterations for the eigenvalue 

calculation 

OMEGAR Dimensionless circular frequency 

OMEGA I Temporal growth rate 

G. PROGRAM OUTPUT 

The output is controlled by iPRriT, see last section. The long print 
contains the short one added to it, the eigenvalues, eigenfunctions of the 
homogeneous problem, and adjoint eigenfunctions at different y-locations. 

The output contains 

N S XC R CML 

NIT, KX, D, DkX, IFL, IWl , IW2 (for TS calculations) 

or NIT,KX, KZ, EF5I, DEF5I, DALFAI, IpL, IV/1 , IW2 (for CF calculations) 

VIGI, VIG2, TATA, EPGR 

KX, KZ, OMEGA, FR, FREQ 

CR, MACH, ETA, IE, WN , PGR, D 

PHI, EPSI, EPSIP 

See section C for the definition of the FORTRAN variable names 
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V. SAMPLE CASES 


The meanflow used for the following sample cases is the same as in reference 2. 

The boundary layer with suction on a 23° swept infinite span wing is 
calculated (see Appendix I). The airfoil section is supercritical with normal 
chord c = 1.98 m including a trailing edge flap. The freestream conditions are; 

Mach number = .82, total temperature = 393 K, and total pressure = .46 atm. 

Figure 3 shows the geometry, the distribution of the pressure coefficient Cp, 
and the distribution of the suction coefficient Cg = -PoV^/pJLI^, where the sub- 
script 0 denotes wall condition. 

Variables needed to be input are marked in the following sample cases. 
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Case. 5.- Calculation of the spatial. growth_rate_pf stationary 
£F disturbance at a .specified location .and r .. 
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Case 6.- Calculation of the maximum spatial growth rate of stationary 
CF disturbance at a specified location with ^ = c- 
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VI. THE CODE HADY-I 
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♦ TIDY* 

PROGRAM HADY(INPUT>0UTPUT> TAPES- INPUT»TAPE6«0UTPUT>TAPE10) 


C 


c 

c 

c 

r 

r 


CO 

-vl 


PROGRAM HADY{ INPUT#QUTPUT,TAPE5-INPUT#TAPE6-0UTPUT»TAPE10) 


* t** «**««** 4141**** ************************** ******* 


* * 

* * 

* HA0Y-I> STABILITY ANALYSIS FOR * 

* THREE-DIMENSIONAL BOUNDARY LAYER FLOWS * 

* * 

* BY ♦ 

* * 

* NABIL M. EL-HAOY * 

* * 

* DEPARTMENT OF MECHANICAL ENGINEERING G MECHANICS ♦ 

* OLD DOMINION UNIVERSITY, NORFOLK, VIRGINIA * 

* * 

* * 


************************************************************ 


REAL MACH,KRR,MU,MUP,MUPP,KFF,NO 

COMPLEX XI,D,DTR,DF,COR,ALAM,TATA, A(8,8 8>,EV(B ) 

COMPLEX KX,KXX,KZ,KZZ,OA, A37, A48, A84, A85, A87,EVV,0MEGA 
COMPLEX API, A2A,A25, A31, A3<», A35,G,AA2,AA3,AA6, A6A, A65, HI, H2 
COMPLEX Z1 (101),Z2(101),Z3(101),Z4(101) ,Z5(101 ),Z6(101),Z7(10l),Z8 
1 ( 101 ) 

COMPLEX Wl(101),W2(101),W3(101),W<t(101) ,W5(101),W6(101),W7(101),W8 
1 ( 101 ) 

DIMENSION. PH(16, 101 ), BMI(8,16), BMF(8,16), BCIV(8), BCFV(8) 
DIMENSION WORK(17000), IW0RK(300) 

COMMON /AAA/ XSAVE, KL, INDEX, NT T 

COMMON /BBB/ KX, KXX,KZ,KZZ, OMEGA, XI, P,Q, GAMMA, OB 
COMMON /B8I/ XC,CHL,UFS,R, MACH,ETA,PRANDL,NO,NS,NY 
COMMON /CCC/ Y(101),U(101),UP(101),UPP( 101),W(101),WP(101),WPP(101 
1 ) ,T(101), TP( 101 ), TPPdOl ), MUdOl ) ,MUP(101 ),MUPP( 101), ALE A (101 ), ALF 
2AP(101),PR(1C1) 

COMMON /FFP/ Z1,Z2,Z3,Z4,Z5,Z6,Z7,Z8 

namelist /SHEE/ ISPTM, INCDMP,NSTATN,NSEND,NSTEP, ICASE, ICIZER,FR, IF 
1R,NFP,0FR,EPSI,IEPSI,NEPSI,DEPSI,EPSIP,IEPSIP,NEPSIP,DEPSIP,ALFAR, 
2ALFAI,BETAR, I BE T AR, B ET Al , WN , I WN , NWN , OWN ,CHORO, IPR I NT, R E , AE , EPS , I TR 
3,nMEGAR,0MCGAI 
JL-1 

XI»(0. ,1. ) 
gamma-1, A 
F-0.8 
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*TinY* 

PROGRAM HADY( INPUTS OUTPUT^ TAPE 5- INPUT,TAPE6«0UTPUT» TAP ElO) 


r 


1 

? 


5 

b 


7 


8 


P« (1. +2.*E ) /3. 

0*?.*(E+2.)/3. 

4c 4c 4t 4c 4: j«r * * 4c * * * * » * :(c * 4: * * * * 4; * 4r « * « 4r « 4: * <t 4: 4: 4 4c « * 4c « 4r 4r 4c 4t 4c * 4c 4r « 4c 4c 4t * # 4c 4i 

reading TN'OIJT data and profiles 

4c4c4'4>4'4c4<4<4cic4c4c4'>f4‘4<4>4c4‘4>4c4c4c4c4<4c4<4'4>*4>4'4<4'4>4c>f4‘4c4<4c>»4c4'4>#4<4:4<4>4c4’4c4>4c4c4>*<t<4c4‘4c4i4c4> 
READ (54SHEF) 

WRITF f 6 ^SHFE) 

CALL PULL“ ( NSTATN» INCOMP) 

0B»GAMMA*MACH*MACH 

FPST«EPSI/57. 29577 

DEPSI»0EPSI/57. 29577 

FPSIP»EPSIP/57. 29577 

DFPSIi>»0EPSTP/57. 29577 

IF (ICASE.EO,3.0P.ICASE.E0.5) GO TO 1 

IF nWN.FQ.2) WN«6. 2831853*CHL/ (WN*CHORD) 

GO TO 2 
CONTINUE 

IF (IBETAR.R0.2 ) BE T AR -6 . 2 831 8 53*C HL / ( B E T AR*CHORD ) 

CONTINUE 

PNI«ATAN(NO) 

IF (IEPSTP.E0.2) EPSIP-PHI 
GO TO lEPSI 

EPSI«PHI+1.57 
GO TP 6 

CALL CFA (EPSI^NY) 

GO TO b 

EPSI«PHI 

CONTINUE 

NIT-0 

NITB-0 

INDEX-1 

TNDOX-l 

IF (ICASE.EQ.3.0R.ICASE.E0.5) GO TO 7 
Al FAR=WN*COS ( EPSI ) 

PETAR-WN*STN(EPSI) 

GO TO 8 
CONTINUE 

IF (EPSI.EO.O.) GO TO 8 
IF EPSI-0. ALFAR IS INPUT 
ALFAR-BETAP/TAN(EPSI ) 

PFTAI-ALFAI*TAN(EPSIP) 

KX=CMPLX{ ALFAR 4 ALFAI) 

K Z-CMPLX ( BETAR 4 BETAI ) 

IF (ISPTM.EQ.2) GO TO 10 


A A6 
A A7 
A A8 
A A9 
A 50 
A 51 
A 52 
A 53 
A 5A 
A 55 
A 56 
A 57 
A 58 
A 59 
A 60 
A 61 
A 62 
A 63 
A 6A 
A 65 
A 66 
A 67 
A 68 
A 69 
A 70 
A 71 
A 72 
A 73 
A 7A 
A 75 
A 76 
A 77 
A 78 
A 79 
A 80 
A 81 
A 82 
A 83 
A 8A 
A 85 
A 86 
A 87 
A 88 
A 89 
A 90 



A T I D Y ♦ 

PROGRAM had Y( IN PUT, OUTPUT* TAPES- INPUT# TAPE 6- OUTPUT, TAPE 10) 


IF (IFP.EQ.l) GO TO 9 A 91 

0MEGAR»6.2831853*FR*CHL/UFS A 92 

OMEGAI-0. A 93 

OMEGA»CMPLX(OMEGAR,OMEGAI ) A 9A 

FPEQ-OMEGA/R A 95 

GO TO 11 A 96 

9 nMEGAR«FR*P A 97 

FPEQ-0MEGAR*UFS/(6.2831853*CHL) A 98 

PMEGAI-0. A 99 

QMEGA»CMPLX(0ME6AR,aMEGAI) A 100 

GO TO 11 A 101 

10 CtHEGA-CMPLX(0MEGAR,0ME6AI ) A 102 

C A 103 

11 WRITF (6,56) NS,XC,R,CHL A 10^ 

IF (ICASE.E0.1.0R.ICASE.E0.3) WRITE (6,53) A 105 

IF (ICASE.E0,A,0R,ICASE,EQ.5) WRITE (6,54) A 106 

12 CONTINUE A 107 

DO 13 I«l,8 A 108 

DO 13 J-1,16 A 109 

BMI(T,J)«0. A 110 

13 BMF(I,J)»0. A 111 

C BOUNDARY CONDITIONS AT THE WALL , HQMOGEN FOUS PROBLEM A 113 

C t****1f***t*^i************^li:¥*************************M*:**M*******'^ A 114 

BMF(1,1)-1. A 115 

BMF(2,2)-1. A 116 

BMF(3,5)«1. A 117 

BMF(4»6)-T. A 118 

0MF(5»9)-1. A 119 

BMF(6,10)»1. a 120 

BMF{7,13)»1. A 121 

BMF(8,14)«1. a 122 

C A 123 

DO 14 I-l.P A 124 

BCFV(I)«0. A 125 

14 BCIV(I)»0, A 126 

r ^,*^:^l^^^**^^^l^Hl^*l^:*^^^f^l^,^^^^^^i^^it****■¥*¥******************************* A 127 

C CALCULATION OF EIGEN VALUES A 128 

0 ^itL^^^ii^t****************************************’^**************** A 129 

15 CONTINUE A 130 

KXX-KX+KX A 131 

KZZ-KZ4KZ A 132 

OA»KX+ND*KZ-OMEGA A 133 

A21«XI*R*0A + KXX + KZZ a 134 

A24-X I*KX*P-P*OB*KX*QA A 135 


* T I D Y * 

pi?OGRAM HADY(INPUT^0UTPUT»TAPE5-INPUT»TAPE6«0UTPUT>TAPE10) 
o 

A25*P*KX*0A A 136 

A31»-XI*KX A 137 

A3A.-XI*OP*oA a 138 

A35*XI*OA A 139 

A37«-XI*KZ A lAO 

r, = P + XI*QB*0*QA A lAl 

4A?.-XI*KX/G a 1A2 

AA3»(-KXX-KZZ-XI*R*QA)/G a 1A3 

AA6*XI*0*0A/G A lAA 

AA8=-XI*KZ/G A 1A5 

A6A--XI*P*PRAN0L* (GAHMA-1. )’*-MACH*MACH*QA A 1A6 

A65-XI *P*ppA NDL*OA+KXX+KZZ A 1A7 

A8A-X I+R*KZ-P*KZ-*-Q3*QA A 1A8 

A85«P*KZ*QA A 1A9 

AR7-A21 A 150 

H1«AA2*A2A + AA3>'‘A3A + AA6*A6A + AA3*A3A A 151 

H2-AA2+A25+AA3*A35+AA6*A65+AA0+A85 A 152 

EV(1)«-CG0PT<A21) a 153 

EV(2)*-CS0RT(0.5>MH1 + A65)+CS0RT(0.25X'(H1-A65)**2 + H2*A6A) ) A 15A 

EV(3)»-CSQRT(0.5*(H1»A65)-CS0RT(0.25*(H1-A65 )**24-H2*A6A) ) A 155 

EV(A)«EV(1) a 156 

EV(5)--EV(l) A 157 

EV(6)»-EV(?) A 158 

EV(7)«-EV(3) A 159 

EV(8)*-EV(A) A 160 

C ^^,^^^,^,^ii^^L^L^i^i*^^^c:^^,^^t********************************^^M■¥**^^**i‘*■¥*** A 161 

C CALCIlLATinN OF FIGEN VECTORS»HQMOGENEOUS PROBLEM A 162 

; *^,,^ii^:^i^i^^L:^^L^*^^:^H:^:l^^i^f^^i^L^L^i^i^t1HLi^^H:***********^‘******************* A 163 

DO 21 J«l,8 A 16A 

EVV«EV(J)*EV( J) A 165 

GO TO (16>18>18»17»16,18»18,17), J A 166 

16 B1(1»J)-1. A 167 

R1(2»J)*0. A 168 

Bl(3iJ)«0. A 169 

Rl(A,J)-0. A 170 

GO TO 19 A 171 

17 B1(1»J)*0. A 172 

Bl(2»J)-0. A 173 

B1(3»J)»0. A 17A 

B1(A,J)»1. A 175 

GO TO 19 A 176 

18 Bid, J )«( (EVV-A66)*A24 + A25*A6A)/(A21-EVV) A 177 

Bl(2, J )«A65-EVV A 178 

B1(3,J)«-A6A A 179 

B1(A,J )-(A6A*A85-A8A>MA65-EVV))/(A21“EVV) A 180 



* T I D Y + 

PROGRA^I HADY( INPUT# OUT PUT. TAPES* INPUT>TAPE6»0UTPJT.TAPE10) 


19 CONTINUE A 101 

IF (J.EQ.4.0R.J.E0.8) GO TO 20 A 182 

A(1.J)*1. A 183 

A(2.J)=EV( J) A 104 

A (3.J )»( A31*B1(1. J )+A34*01 (2. J ) + A3 5*Bl{ 3. J ) + A37*Bl(4, J ) ) /(EV(J )tBl A 105 

1( 1. J) ) A 106 

A(4.J)*31I2. J)/B1(1.J) A 107 

A(5.J)*B1(3» J)/B1(1»J) A 108 

A (6. J)»EV( J) *B1 (3. J )/Bl(l. J ) A 189 

A(7,J)*B1(4, J)/B1(1» J) A 190 

A(8,J)»(A34+B1(2.J)+A85*B1(3»J)+A07*B1(4,J))/(EV{J)+B1(1,J)) A 191 

GO TO 21 A 192 

20 A (1, J )»81 (1. J ) A 193 

A(2,J)«EV(J)*01(1.J) A 194 

A(3»J)»( A31*B1(1»J)+A34*B1(2» J)+A35*B1{3»J)+A37*B1(4» J)) 7EV( J) A 195 

A(4, J )«B1 (2, J ) A 196 

A ( 5. J )«81 n. J) A 197 

A(6,J)»EV(J)*01(3.J) A 198 

A(7»J)-B1(4.J) A 199 

A (8. J ) «( A84*B1 (2. J ) +A05*B1 (3. J ) +A87*B1( 4, J) ) /E V( J) A 200 

21 CONTINUE A 201 

GO TO 29 A 202 

22 CONTINUE A 203 

no 23 T»1.8 A 204 

on 23 J*1.16 A 205 

BMni,J)=0. A 206 

23 BMF(I.J)»0. A 207 

:tfii^,i,*t********************************************************** A 208 

C BOUNDARY CONDITIONS AT THE WALL. ADJOINT PROBLEM A 209 

r ^,4^f^,^4:ll:^*i,*t:ti^,^:^,ft^,*^t:^ii^it4^,fft*************t***************f*** A 210 

BMF(1,3)«1. a 211 

RMF(2.4)«1. A 212 

BMF(3.7)»1, A 213 

BMF(4,B)»1. A 214 

BMFC^.ll ) = 1. A 215 

BMF(6.12)*1. a 216 

RMF(7.15)=1. A 217 

BMF(8.16)«1. A 218 

» ^if:H,^,f:^ii4,**^,^HL************************************************* A 219 

: CALCULATION OF EIGEN VEC TORS, AO JOINT PROBLEM A 220 

C ^^:l^4ti^^^^i^,f1^44i^f^^r|^if^i■i^iHi^^if^^t^t**4****************************** A 221 

DO 28 J*l,8 A 222 

FVV»EV(J)+EV(J) A 223 

GO TO (^4, 25,25, 26. 24,25,25,26), J A 224 

24 Bm,J)»l. A 225 



^ * T I 0 Y * 

PROGRAM HADYdNPUT, OUTPUT, TAPE5»INPUT*TAPE6«0UTPUT»TAPE10) 

Bl( ?, J )=• ( + (A6i-EVV)-A6A*A25 ) / ( H2 * A6A- { Hl-EV V) * ( A65-E V V ) ) A 22 6 

PI (3, J )*( A??* (H1-EVV)-A2A*H2) / (H2*A6A-( H1-EVV)»{ A65-EVV) ) A 227 

Bl(A,J)-0. A 228 

GO TO 27 A 229 

25 B1 (1, J )-0. A 230 

P1(2,J)=>-A6A a 231 

R1 (3, J VV A 232 

B1 (A, J )«0. A 233 

GO TO 27 A 23A 

?6 Rl(l,J)-0. A 235 

B1 (2, J )- ( ABA* ( A65-EVV)-A6A*A85) /(H2*A6A-(H1-EVV)*(A65-EVV) ) A 2 36 

B1 (3, J ). ( AB5* (Hl-EVV )-A8A+H2 ) /( H2*A6A-( Hl-EVV) *( A65-EVV) ) A 237 

B1 (A, J )*1. A 230 

27 CONTINUE A 239 

A ( 1, J ) »B1 (I, J ) A 2A0 

A ( ?, J ) »(-Rl( 1, J i-AA2*Bl ( 2, J) ) /EV( J ) A 2A1 

A(3,J)=-AA3*B1(2, J)/EV(J) A 2A2 

A (A,J )-Bl (2, J ) A 2A3 

A ( 5, J ) »B1 (3, J ) A 2AA 

A(6,J)»(-AA6*B1(2,J)-B1(3,J))/EV(J) A 2A5 

A (7, J ) -B1 (A, J ) A 2A6 

28 A ( 8, J ) »(-AAP*Bl (2, J)-31( A, J ) ) /EV( J ) A 2A7 

2R CONTINUE A 2A8 

C ^*]»t4'*****t******«***>»******************************************* A 2A9 

C BOUNDARY CONDITIONS AT THE EDGE OF THE BOUNDARY LAYER A 250 

; t***************>«^******«*4i***********i**«******i:»***i*************** A 251 

CALL CDMTNV (A,8,DTR) A 252 

►'“A A 253 

DO 30 1 = 1.7, 2 A 25A 

L»0 A 255 

N.M+1 A 256 

no 30 J=l,15,2 A 257 

L»L+1 A 258 

BMKT, J) = A(M,l) a 259 

30 BMI(I,J+l)=yi*A(M,L) A 260 

DO 31 I»2,R,2 A 261 

DO 31 J=1.15,2 A 262 

RMI(I,J)»-RNI(I-1,J+1) A 263 

31 RNKI, J + l)aBM(I-l,J) A 26A 

: - A 265 

IWORK(ll)«l A 266 

IW0RK(l)=NY/6 A 267 

II«IW0RK(1) A 260 

DO 32 1=1, TT A 269 

32 WORK ( I )=Y (6*1-1 ) A 270 



* T I D Y + 

PROGRAM HADY( INPUT>0UTPUT.TAPE5»IK'PUT»TAPE6»0UTPUT*TAPE10) 


XSAVE»0. A 271 

KL»1 A 272 

IFLAG«0 A 273 

C **t***********i*^***********i{‘***4*HL*****M***********'¥*'¥*******1f*4** A 274 

C INTEGRATION AND ORTHONORMALIZ ATION A 275 

; ^^i^t***t*********’¥**¥**********************************t*^f**irM*** A 276 

CALL JOCK (PH>16>16,Y,NY,BMI,8^BCIV»8>BMF,8>BCFV,8>0,RE» AE*IFLAG»W A 277 

10PK»17000»IW0RK»300#0) A 278 

C A 279 

DF«CMPLX( PH( 3>NY)#PH(4>NY) ) A 280 

COR-l./OF A 281 

D«rMPLX(PH{13»NY),PH(14,NY))#C0R A 282 

MNB«TWGi?K(l) A 283 

MVC«IW0RK(2) A 284 

GO TO (33.33»37>50)» INDEX A 285 

33 CONTINUE A 286 

IF (ICASE.E0.4.0R.ICASE,E0.5) GO TO 34 A 287 

CALL RNITN ( E P S I , E P S I P > I F L AG> MN B , M VC > D, ISP TM, I C AS F » WN , E P S » ITR ) A 288 

GO TO 35 A 289 

34 CALL CFTTN ( E PS I > E PS I P, I F L AG# MN3 , M VC » 0» WN, E PS» IC ASE » BE T AR » I TR ) A 290 

35 IF (AIHAG{KX) ,GT. 0.0010. AND. lEPSI.NE.l) GO TO 46 A 291 

IF (INDEX .EQ.l. OP. INDEX. EO. 2) GO TO 15 A 292 

CP»PFAL(nMFGA ) /REAL (KX ) A 293 

DO 36 T«1>ny a 294 

71(1) «CMPLX( PH( 1» I ) » PH(2» I ) )*COR A 295 

Z?(I)«rMPLX(PH(3»n.PH(4»I))*C0R A 296 

Z3(I)«CMPLX(PH(5^ I).PH(6»I ))*COR A 297 

Z4(I )»CMPLX( PH(7, I ) . PH(8» I ) )*COR A 298 

Z5(I)*C«PLX(PH(9W),PH(10,I))*C0R A 299 

Z6(T)»CMPLX(PH{11#I)*PH(12#I))*C0R A 300 

Z7( I) .CMPLX{ PH( 13^ I ) . PH( 14. I ) ) +C OR A 301 

ZR(I)»CMPLX(PH(15,I).PH(16.I) )*COR A 302 

36 CONTINUE A 303 

; ^^*^|li^:^l^*^I^l^i^^;^;^^;^l^^Lt******^‘***************l^**^^******^^*****Al*****^*M A 304 

C START THE ADJOINT PROBLEM A 305 

t4i^^it*****t**i*^<¥***m******’t‘******************^^*****^^**********^^** A 306 

INDEX. 3 A 307 

GO TO 22 A 308 

37 CONTINUE A 309 

DF.CMPLX( PM( 1,NY) .PH(2.NY) ) A 310 

COR. 1. /DC A 311 

DO 38 I.l.NY A 312 

W1(I).CMPLX(PH(1,I).PH(2.I))*C0R A 313 

M2 ( I) «C^PLX( PH(3» I ) . PH(4» I ) )*C0R A 314 

W3 (I ).CMPLX(PH( 5. I ). PH(6. 1 ) ) + C0R A 315 


I 


♦ TIDY* 

P?OG‘?ft'< HADY( INPUTS OUT PUT, TAPE 5* INPUT »T APE 6«0UTPJT,TAPE10) 


W4(T)-CMPLX(PH(7,I),PH(3,I))*C0R A 316 

W5(I)»C^^PLX(PH(9,I),PH(10,I))*C0R A 317 

W6(I) *CMPLX( PH(ll,n,PH(12,I) )*COR A 318 

W7 (I ) •Cf'Pl X(PH (13, I ) ,PH(14,I ) )*COR A 319 

38 WP (I) »CMPLX( PH(15, I ), PH(16,I ) )*CQP A 320 

C A 321 

ALFAI»-XI*KX A 322 

PFTAI--XT*KZ A 323 

ALFAP-KX A 324 

PFTAP*K7 A 325 

WM«SOPT( ALFAR**2 + BETAR**2 ) A 326 

TF ( TC AGE .ED. 2) E PS I P» AT A N ( BE T A I / AL F A I ) A 327 

IF (ICASF.FQ.3) EPSI»ATAN(BETAR/ALFAR) A 328 

“ t*t*>)<t***4'***t4'*>tMf4'««t***t***«4>************************4'********* A 329 

: CALCULATION OF THE RATIO OF THE GROUP VELOCITY COMPONENTS A 330 

; A 331 

IF (CABS(EV(2)).LT.CABS(EV(3))) GO TO 39 A 332 

ALAN.EV(3) a 333 

IA*3 A 334 

GO TP AO A 335 

39 ALAM=FV(2) a 336 

IA*2 A 337 

40 CALL GRVEL ( A L A M, W2 , W3, W4 , W6 , W 8, T AT A , E P GP , SM ) A 338 

IP (I5PTM.E0.2) GO TO 41 A 339 

PGP--ALFAI-8ETAI+SM A 340 

GO TO 42 A 341 

A 1 PGPrnMFGAI A 342 

FP«OMFGA/R A 343 

PPFQ*OMEGAS*UFS/ (6.2031853+CHL) A 344 

42 WRITE (6,55) KX , K Z , OME G A, F R, F RE Q, C R , M AC H, E T A ,N Y, WN, PG R , Z 7 ( N Y ) , PH I* A 345 

157, 29577, FPSI*57. 29577, EPSIP*57. 29577 A 346 

GIN.-XI+TATA A 347 

r i^^t************************************************************** A 348 

C ITERATION FOR MAXIMUM AMPLIFICATION PATE A 349 

; t^L4,^!^^***********^f*********************■********i^*******^^*^^*^^^^*¥**^^ A 350 

IF (ICIZFR.NE.l) GO TO 43 A 351 

IF (A3S(GIM) .LE.I.E-05) GO TO 43 A 352 

CALL GIMTR ( G I M, E PSI , E PSI P , IC A SE , N I TB , I NOOX, WN, E PS ) A 353 

MIT.O A 354 

INDEX. 1 A 355 

IF (INDOX.NF .10) GO TO 12 A 356 

43 CONTINUE A 357 

IF nPRINT.FQ.2) GO TO 46 A 358 

WRITE (6,58) A 359 

WRITE (6,60) EV(1),EV(2),EV(3),EV(4) A 360 



* T I D Y ♦ 

PRQGO HADY(INPUT>0UTPUT,TAPE5=*INPUT,TAPEb*0UTPUT»TfcPE10) 


WPITE (6,61) A 361 

on AA I-1,NY,5 A 362 

WPITE (6,57) Y(I),Zl(n,Z3(I),Z5(I),Z7(I) A 363 

WRITE (6,57) Y(NY),Z1(NY),Z3(NY),Z5(NY),Z7(NY) A 364 

WRITE (6,59) A 365 

DO 45 I»1,NY,5 A 366 

45 WRITE (6,57) Y ( I ) , W2 ( I ) , W4 ( I ) , W6 ( I ) , W8 ( I ) A 367 

WRITE ( 6 , 57 ) Y(NY),W2(NY),W4(NY),W6(NY),W8(NY) A 360 

46 :PNTIMUE a 369 

r *******^>t******<if********’^*********t*******************t*4t**f^*** A 370 

C AMPTHER FREOUENCY/WAVE ANGLE /WAVENUMBER /DOWNSTRE AM STATION A 371 

IF (MPR.GT.1) go to 47 A 373 

IF (NFPSI.GT.l) GO TO 48 A 374 

IF (NWN.GT.l) GO TO 49 A 375 

NSTATN»NSTATN+MSTEP a 376 

IF (NSTATN.GT.NSENO) stop a 377 

KRP-WN/CHL A 378 

RKP.BETAR/CHL a 379 

CALL FULLP ( NSTATN, INCOMP ) A 380 

WN=KPR*CHL*0.99 A 381 

IF (ICASE.EQ.4) WN-KRR*CHL A 382 

TF (ICASF.FO,3,OR.ICASE.EQ.5) 3ETAR-RKR tCHL A 383 

GO TO 2 A 384 

47 CONTINUE A 385 

JL-JL+1 A 386 

TF (JL.GT.NFR) STOP A 387 

KFF-WN/PP A 388 

FR.FPfOFR A 309 

WN=KFF*FR*0.97 a 390 

GO TO 6 A 391 

48 CONTINUE A 392 

Jl « JL f 1 A 393 

IF (JL.GT.NEPSI ) STOP A 394 

FPSI=FPSI+DFPSI A 395 

GO TO 6 A 396 

49 CnNTTNUE A 397 

JI.»JL+1 A 398 

IF (JL.GT.NWN) STOP A 399 

WN*WN+nWN A 400 

GO TO 6 A 401 

WPITE (6.52) A 402 

^ STOP A 403 

un 50 WPITF (6.51) A 404 

STOP A 405 



^ * T I 0 Y + 

O' PROGRA'^ HAPY(INPUT>OUTPUTiTAPE5»INPUT,TAPE6«OUTPUT,TAPE10) 


C A A06 

51 FORMAT (IM flOHHA HA HA) A A07 

52 format (1H >26H N0RM«2) A A08 

53 FORMAT (IH i lX^3HNIT,13X,2HKX,24X,lHD.22X,3HDKXn2X>3HIFL»lX»3HIWl A A09 

1,1X,3HIW2/) A 410 

54 format (1H , lXi3HNIT,13X ,2HKX,24X,2HKZ» 17Xf 4HEPSI»8X, 5HDEPSI.0X, 6H A 411 

1DALFAI,5X»3HIFL^1X#3HIW1>1X>3HIW2/) A 412 


55 FORMAT (IH » 5 X^ 3HKX » , 2 ( E 1 1 . 5» 1 X ) , 3HK Z«» 2 ( E 11 . 5 » 1 X ) , 6H0ME G A« , 2 ( E 1 1 . A 413 
15,1X) » 3HTP = , F11.5> IX, 5HFRE0=» Ell .5//»5X,3HCR«> eil.5,lX,5HMACH»»Ell A 414 
2. 5,1X , 4HETA=, E11.5. 1X,3HIE I5,1X, 3HWN=,E11.5, 1X,4HPGR>,E11.5,1X,2 A 415 
3Hn«,?(E11.5,lX)//,5X,4HPHI»»F9.4,lX»5HEPSI-»F9.4,lX,{>HEPSIP».F9.4) A 416 

56 format (IH , ////,5X,3HNS-, 15, 5X,3HXC-,E11.4,5X#2HR»,E11.4,5X,4HCHL A 417 


1=,F]1.4/) A 418 

57 FORMAT (IH , 1 X, F 5 . 2 , IX , 6 ( 09 . 3 , 1 X , E 9 . 3i 1 X ) ) A 419 

53 FORMAT (IH ,1?HEIGEN VALUES/) A 420 

59 format (IH ,1X,42HEIGEN F UMC TI ONS- AD J 01 NT ( Wl» W2» W3» W4, W6» W8 ) ) A 421 

60 FORMAT (IH , 5 X , 10 ( 0 1 1 . 4, 1 X ) / / / ) A 422 

61 FORMAT (IH ,42HEIGEN FUNCT lONS-REGUL AR ( Z 1» Z2» Z 3» Z5, Z 6, Z 7 ) ) A 423 

END A 424- 


* T I D Y * 

SUBROUTINE FULL? (NW, INCOMP) 

SUBROUTINE FULLP <NW»INCQMP) B 1 

REAL MACH^MU# MUP^ MUPPjND B 2 

COMMON /BBl/ XC»CHL,UFS* R»MACH,ETA,PRA,ND»NS»IE B 3 

COMMON ACCC/ YI ( 101 ) , U ( 101 ) i UP ( 101 ) , UPP ( 101 ) , W (101 ) , WR ( 1 01 ) » WPP ( 10 B U 

11 )>T(101 ),TP(101 )iTPP(101 )»MU(101 )#MUP( 101)>MUPP(101)> ALFA(101),AL B 5 

2FAP(101)^ PR(lOl) B 6 

1 CONTINUE B 7 

RFAD (10) N'>,XC^MACH,R,CHL>UFS^ET A,PRA»ND»IE B 8 

READ (10) (YI (I )#I«1,IE ) B 9 

RFAD (10) (U( I )# I-l» IE ) 3 10 

READ (10) (UP(I)»I«1.IE) B 11 

READ (10) (UPP( I ) , 1*1, IE ) B 12 

READ (10) (W(I),I-1,IE) B 13 

READ (10) (WP(I),I»1,IE) B 1<» 

read (10) (WPP(I),I*1,IE) B 15 

RFAD (10) (T( I), I-l, IE ) B 15 

RFAD (10) (TP( I ), I»l, IE) B 17 

READ (10) (TPP(I),I*1,IE) B 18 

READ (10) (MU(I), I«1,IE) B 19 

READ (10) (MUP(I),I*1,IE) B 20 

READ (10) (MUPPd ), I=1,IE) B 21 

READ (10) (ALFA(I), I«1,IE) 8 22 

read (10) ( ALFAPd ), I*1,IE ) 8 23 

READ (10) (PR(I),I*1,IE) B 

IF (NS.LT.NW) GO TO 1 B 25 

ND*W(1) B 26 

IF (INCOMP. NE.l) RETURN B 27 

MACH.0.0001 8 28 

DO 2 1*1, IE 3 29 

Td)-1. B 30 

TP(T)*0. B 31 

TPPd)=0. B 32 

MU(T)«1. B 33 

MUP(T)*0. B 34 

MUPP(I)»0. B 35 

ALFA(I)-O.B B 36 

2 AlFAP(I)*0. B 37 

RETURN B 38 

FMO B 39- 





-> 

00 


♦ TIDY* 


SUBROUTINE FMAT ( X, Y» YP » IGOF X, S# SP ) 


SURRnUTINE FMAT ( X, Y» YP> I GQFX > S > S P ) C 1 

RFAL '^ACh, MM, MUP> MUPP» NO C 2 

CnMPLFX Gl,G2#G3^GA,G5,G6^G7>G8,Zl,Z2»Z3f Z4,Z5»Z6»XI>KX»KXX,KZ,KZZ C 3 

lt0A,0C»0'^FGA,Z7>Z8 C A 

COMPLEX A?l,A23^A2A.A25,A31>A3‘t,A35,G»A41,A^2»A<t3»AA«r, AA5» A«»6 C 5 

COMPLEX A63/ A6<»/ Afe5^ A37/A'f3»A3^,A35/A87/A47/ A83 C b 

DIMENSTON Y( 1), YP(i), S(l)> SP(1) C 7 

CO'^MON /AAA/ XSAVE»KL> INDEX, NIT C 8 

COMMON /PRP/ KX, KXX,KZ,KZZ, OMEGA, XI, P,g, GAMMA, OB C 9 

COMMON /BPl/ XC,CHL,UFS,R, MACH, ETA, PRA, ND,NS>N C 10 

C C 11 

IF (X<?AVE.FO.X) GO TO 1 C 12 

IF (XSAVE.LT.X) KL»1 C 13 

XSAVE = X C !•* 

CALL PROF (X,U1,U1P,UPP,W1,W1P,'VPP,T,TP,TPP,MU,MIJP,MUPP,ALFA,ALFAP C 15 

1, PPAMDL,KL , INDEX ) C 16 

OA»KX*Ul+KZ*wl-OMEGA C 17 

0P»Ky*UlP+K7*WlP C 18 

OD»PP ANnL*(GAMMA-l. )*MACH*MACH C 19 

A?1»XI+P*0A/ (T*MU)+KXX+KZZ C 20 

A2?=-MUP/MU C 21 

A?3-P*U1P/ (T*MU)-XI*KX*(MUP/MU+P*TP/T) C 22 

A2A»XT*R*KX/MU-P*Q3*KX*QA C 23 

A25»P*KX*GA/T-(ALFA*UPP+U1P*ALFAP) /MU C 2A 

A26«-ALFA*U1P/MU C 25 

A31*-XI*KX C 26 

A33*TP/T C 27 

A3A.-XI*QB*0A C 28 

AB'j-XI+OA/T C 29 

A37=-XT*K7 C 30 

G»P/MU+XI*0B*0*0A C 31 

A A1.-XI*KX*( 2 .*MUP/MU+0+TP/T) /G C 32 

A42.-y[+KX/G C 33 

AA3«{-KXX-KZZ+(0*MUP*TP/MU+Q*TPP-XI*R+0A/MU)/T)/G C 3^ 

/I A<,a-xr*0*on*( OA*(MUP/MU + TP/T )+QC ) /G C 35 

A45«X T*( ( ALFA/MU+O/T )*QC+Q*MUP*QA/ (T*MU ) ) /G C 36 

AA6»XI*0*0A/ (T*G) C 37 

AA7»-XI+KZ*(2,*MUP/MU+Q+TP/T)/G C 38 

AA8.-XI+K7/G C 39 

A62--? .*On*IJlP C AO 

A63»P*PRAN0L*TP/( T*MU)-2.*XI*QD*(3C C A1 

A/A*-X I*R*QD*QA/MU C A2 

A65»XT*R*PRANDL*QA/(T*MU)-ALFA*0D*(U1P*U1P + W1P*W1P)/MU4-KXX4-KZZ~MUP C A3 

IP/MU C AA 

A6f«-?. + M(|P/MU C A5 



* T I D Y * 

SUBROUTINE FMAT ( X> Y» YP » IGOF X» S» S P ) 

A68»“2.*Q0*W1P C 46 

A63*P1'W1P/(T*-MU)-XI*(MUP/MU+P*TP/T)*KZ C 47 

A84*K7+(XI+P/MU-P*QB*QA) C 48 

APf5«-( ALFA*WPP + WlP*ALFAP)/MU + P*KZfQA/T C 49 

AP6»-ALFA*W1 P/MU C 50 

AP7»A?1 C 51 

AB8--MIJP/MIJ C 52 

1 :nNTINUE C 53 

G1*CMPLX( Y(1 ),Y(2) ) C 54 

G?*CMPLX(Y(3)^Y(A) ) C 55 

G3«CMPLX ( Y(5 ) , Y (6 ) ) C 56 

G4«CMPLX( Y(7 ),Y(P ) ) C 57 

GS^Cr'PLXC Y(0) ,Y(io) ) C 58 

G6-rMPLX(Y(n )>Y(12) ) C 59 

G7.rMPLX(Y(13),Y(14) ) C 60 

GP*CMPLX(Yn5),Y(16) ) C 61 

TF (IN0FX.EO.3.nR.IN0EX.EQ.4) GO TO 2 C 62 

YP(l).vn) c 63 

YP(?)«Y(A) C 64 

71»A21*Gi+A?2*G2+A23*G3+A24*G4+A25*G5+A2b*G6 C 65 

YP(3)«Z1 C 66 

YP{<t).-XT*71 C 67 

7?»A31*G1+A33*G3+A34*G4+A35*G5+A37*G7 C 68 

YP(5)»Z2 : 69 

YP(6)*-XT*7? C 70 

73*AA1*G1+4A2*G2 + A43*G3 + A4<»*G4 + A45*G5*-A46*G6 + A47*G74-AA3*G8 C 71 

YP(7)»73 C 72 

YP(8)«-Xr*73 C 73 

YP(0)»Y(11) C 74 

YP(10)»Y(1?) C 75 

74,«Af,?*GP + 4b3»G3 + A6^*G4 + A65*G5 + A66*G6 + A6e+68 C 76 

YP(11)=74 C 77 

YP(1?)»-XT*ZA C 78 

YP(13)-Y(15) C 79 

YP(1<*)=Y(16) C 00 

7 5»A83*G3+ARA*G4+A85*G5+A86+G6+A87*G7+A 88+G8 C 81 

YP(15)»Z5 C 82 

YPn6)=-XT*Z5 C 83 

PFTUPN C 84 

: C 85 

P continue C 86 

Z1--A?1*G2-A31*G3-A<^1 + G4 C 87 

YP(!)»71 C 88 

YP(?)«-XT*Z1 C 89 

72»-GI-A22*G2-A42+G4-A62+G6 C 90 


vO 



Ln 

O 


♦ tidy* 


SUBROUTINE FMAT ( Xi Y> YP > I60F X, S# S P ) 


YP(’) »Z? 

YP ( ^ ) .-XT *7? 

7 3*-/'P3*G2-A3 3*G3-A<>3*G4-A63*G6-A8 3*Ge 
YP (f ) «Z3 
YP(6) «-XTt73 

74 .-A 2 4>ltG?-A 3A*G3-A4't*G<»-A6^*G6-A8«t*G8 

YP(7)«Z4 

YP(8) x-XT*Z^ 

Z 8--A P'5*G2-A3 5*G3-A^5*G^-A55 + G6-A8 5*G8 

YPiqi-ZB 

YP(10)x-XT*75 

Z<^=-A ?6’*'G2-A A6*G4-G5-A66*G6-A86*G8 

YPdl )-76 

YP(l?).-XT*7ft 

Z7.-A3 7*G3-A4 7*G^-A8 7t'G8 

YP (13 )-Z7 

YP(1A ).-XTYZ7 

ZR.-AA8*GA-A68*G6-G7-A68*G8 
YP (15 ) «Z‘^ 

VP(16).-XT*78 

RETURN 

END 


C 91 
C 92 
C 93 
C 94 
C 95 
C 96 
C 97 
C 98 
C 99 
C 100 
C 101 
C 102 
C 103 
C 104 
C 105 
C 106 
C 107 
C 108 
C 109 

: 110 

C 111 

c 112 - 



♦ TIDY* 

Sl'PPOUTINE CFITN ( E PS I j E PS IP » I FL AG» HNB» MVC» D>WN» EPS » IC ASE » BET AR» IT 
SUBROUTINE CFITN ( £ PSI » E P S I P» I FL AG> MNB» M VC# D> WN> EPS» I C AS E» BET AR» IT D 1 


IR ) 0 2 

COMPLEX KX,KZ#OMEGA»0#KXXi KZZ» XI 0 3 

COMflON /AAA/ XSAVE#KL# INDEX#NIT 0 A 

COHMON /3B3/ KX#KXX#KZ#KZ2#0NEGAi Xl» P#Q#GAMMA, QB D 5 

DEL-.0005 0 6 

GO TO (1#3,1#3)# INDEX 0 7 

1 PAR-n 0 8 

DAI«-XT*D 0 9 

FPSIA*EPST 0 10 

ALFAIA— XI+KX 0 11 

DEPSI.EPSI+nEL D 12 

IF (DEPSI.EO.O. ) DEPSI*DEL 0 13 

FPSI-EPSI+PPPSl 0 14 

IF (ICASE.E0.5) GO TO 2 0 15 

BFTAR«WN*STN( FPSl ) D 1(> 

2 ALFAP»3ETAR/TAN(EPSI) 0 17 

PETAI«ALFAIA*TAN(EPSIP) 0 18 

KX«CMPLX( ALFAR# ALFAIA) 0 19 

KZ«CMPLX( 3ETAR# BETAI ) 0 20 

INDEX-2 0 21 

INOnx-1 0 22 

RETURN 0 23 

3 continue D 24 

IF (IND0X.E0.2) GO TO 5 0 25 

DFPS- ( REAL (0 )-DAP )/DEPSI 0 26 

0FIS-{ AIMAG(0)-DAI)/0EPSI 0 27 

EPSI-EPSIA 0 28 

DALFAI -ALF AT A*DEL 0 29 

ALPAI- ALFAIA+OALFAI 0 30 

IF (ICASE.EO.5) GO to 4 0 31 

BFTAP-WN + SIN(EPSI ) 0 32 

4 ALFAR-BETAR/TAN(EPSH 0 33 

BETAT»ALFAT*TAN(EPSIP) 0 34 

<X-C‘1PLX( ALFA(<#ALFAI) D 35 

KZ»CMPLX(BETAR#BFTAI) 0 36 

NTT-NTT+1 D 37 

IF (NTT.GE.ITP) GO TO 11 0 38 

INDEX-2 0 39 

INDOX-2 0 40 

RFTIJPN 0 41 

5 OFRAI-(RFAL(D)-DAR) /OALFAI 0 4? 

DFIAI*(AIMAG(0)-CAI) /DALFAI 0 43 

SAL»nFRS*DFTAI-DFRAI*DFIS 0 44 

IP (SAL.FO.O, ) GO TO 8 D 45 


I 


T I D Y + 

SU'^PGUTINF rpiTN (EPSI^F.PSIP>IFLAG»MNB»MVC#0,WN,EPS,i:4SE,9ETAP,IT 

DEPSI» (-DAP*DFIAI*DAI*DFRAI)/SAL D A5 

DALFAI*(DAP*DFIS-DAI*0FPS ) /SAL D 47 

WRITE (6,13) NIT, KX,KZ,EPSI*57. 29577, DEPSI+57. 29577, DALFAI,IFLAG,M D 48 

1NB,MVC D 49 

IF ( ABS(DE°S I ) ,LE .EPS. AND. ABS (OALFAI ) .L E.EPS ) GO TO 10 D 50 

5 IF (ABS(DFPSI ) ,LE . ( ,02*A9S ( EPSI ) ) ) GOTO? 0 51 

DEPSI-.5+DEPSI 0 52 

GO TO 6 0 53 

7 IF ( A3S(D ALF AI ) .LE . ( .5*ABS ( ALF AI ) ) ) GO TO 8 0 54 

DALFAI=.5*DALFAI 0 55 

GO TO 7 0 56 

8 CONTINUE 0 57 

e PSI-E PST+DFPSI 0 58 

ALTAI- ALFAT+DALFAI D 59 

IF (NIT.FO.TTR) GO TO 11 0 60 

IF (ICASE.E0.5) GO TO 9 0 61 

RFTAP»WN*SIM( EPSI ) D 62 

9 ALTAR -BETAR/TAN ( EPS I ) 0 63 

BETAI»ALFAI*TAM (EPSIP) D 64 

Ky«rNPLX( ALFAR, ALFAI) D 65 

K7-CMPLX( BETAR, BETAI ) 0 66 

INDFX-1 0 67 

PFTUPN D 68 

10 INDFX-10 0 69 

GO TO 12 D 70 

11 WPITF (6,14) 0 71 

STOP 0 72 

12 RETURN D 73 

: 0 74 

13 format (1M , IX, I2,7E13.4,1X,3( 13, IX) ) 0 75 

14 FORMAT (IH ,27HEXESSIVE NO. OF ITTERATIONS) 0 76 

END 0 77- 





♦ TIDY* 

SUBROUTINE PNITN (EPSI^EPSIP^IFLAG#MNB» M VC » D* I SPTM» IC A SE , WN, EPS ♦ IT 
SUBROUTINE PNITN ( E P S I , E PS I I EL AG , MNB» MVCt D» I SPTH» I C A SE » WN, E PS » IT E 1 


IP) E 2 

COMPLEX D,nA,OK,KX>KZ#OMEGA,DPHOK,K,XI,KXX,KZZ E 3 

COMMON /AAA/ XS A VE » KL , INDE X , N I T E ^ 

COMMON /BRB/ KXf K XX, K Z » K Z Z » OM E G A > X I » P» 0 » 6AMM A> OB E 5 

C E 6 

DEL*. 005 E 7 

GO TO (1,7, 1,7), INDEX E 3 

1 DA«D E 9 

IF (ISPTM.E0.2) GO TO 6 E 10 

nK»KX*DEL E 11 

KXxKX+OK E 12 

ALFAR-KX E 13 

ALFAI--XT+KX E 

RF.TAP»KZ E 15 

BETAI»-XI*KZ E 16 

GO TO (A, 3, 2), ICASE E 17 

2 RFTAI»ALFAI*TAN (EPSIP) E 10 

EPCI.aTAN(RETAP/ALFAR) E 19 

GO TO 5 E 20 

3 BETAP=ALFAR*TAN(EFSI ) E 21 

GO TO 5 E 22 

A RETAR-ALF AR^TAN (EPS!) E 23 

BFTAI *ALFAI*TAN (EPSIP) E 2A 

5 KZxCMPLX(RFTAR,BETAI) E 25 

NIT-NIT+1 E 26 

INDEX»2 E 27 

RETURN E 28 

6 CONTINUE E 29 

DK»OMEGA+DEL E 30 

OMEGA-OMEGA+DK E 31 

NIT=NIT+1 E 32 

INDEX*2 E 33 

RETURN E 3A 

E 35 

CONTINUE E 36 

DPHDK* (D-DA) /DK E 37 

PK«-n/OPMDK E 38 

IF (TSPTM.EQ.I) K»KX E 39 

IF (TSPTM.F0.2) K«OMEGA E AO 

IF (CARS(DK/K).IE.EPS) GO TO 15 E A1 

R IF (GABR(DK) ,LE.(0.2*CARS(K) ) ) GO TO 9 E A2 

DK*0.5*0K E A3 

GO TO 8 E AA 

9 CONTINUE E A5 


I 



* T T n Y 

SUBROUTINE pNITN ( E P S I# E PS I P* I FL AG» MNB , MVC» D» I SPTM, 1C ASE » WN> EPS » I T 


WPITF (6»18) NIT>K,0,DK»IFLAG>MN8»MVC 

IF (ISPTM.EO.Z) GO TO lA 

KX»KXfOK 

ALFAP=KX 

ALFAI--XT+KX 

PETAP*KZ 

BFTAI*-XI*KZ 

GO TO ICASE 

10 RFTAI«ALPAT*TAN(EPSIP) 
FPSI«ATAN(BETAR/ALFAR) 

GO TO 13 

11 RETAR*ALPAP*TAN(EPSI ) 

GO TO 13 

12 BFTAP*ALFAR*TAN(EPSI) 
BETAI»ALFAT*TAfl(EPSIP) 

13 CONTINUE 
KZ»CMPLX(BETAR^ BETAI ) 

R A = 0 

IF (NTT.EO.TTR) GO TO 16 

NIT-NIT+1 

INDEX.? 

RETURN 

lA CONTINUE 

nr'EGA=ONEGA+DK 
D A-n 

IF (NIT.EO.TTR) GO TO 16 

NIT.NT T+1 

INDFX.2 

RETURN 

C 

15 CONTINUE 
INDEX. 10 
GO TO 17 

16 WRITE (6,19) 

STOP 

17 RETURN 

C 

IB FORMAT (1H,1X,I2,6E13.A,1X,3(I3,1X) ) 

19 FORMAT (1H,27HEXESSIVEN0.QFITTERATI0NS) 

FND 


E A6 
E A7 
E A8 
E A9 
E 50 
E 51 
E 52 
E 53 
E 5A 
E 55 
E 56 
E 57 
E 58 
E 59 
E 60 
E 61 
E 62 
E 63 
E 6A 
E 65 
E 66 
E 67 
E 68 
E 69 
E 70 
E 71 
E 72 
E 73 
E 7A 
E 75 
E 76 
E 77 
E 78 
E 79 
E 80 
E 81 
E 82 
E 83 
E 8A 
E 85 
E 86- 


* T I f) Y * 


SUHRni.'TINE GIMTR ( G I M , F PS I » E PSI P » I C ASE, NI TB , INOOX, WN, EPS ) 


SUBROUTINE GIMTR ( GI N# EPS I , E PS I P# IC ASE# N 1 TB^ IND0X> WN, E PS ) F 1 

COMPLEX KZ»KX,XI,0MEGA»KXX,KZZ F 2 

COMMHN /BBR/ KX>KXX,KZ»KZZ^0MEGA,XI,P,0,GAMMA,0B F 3 

C F 

0FLX-.OOO5 F 5 

ITP-10 F 6 

GO TO (1»6)» INDOX F 7 

1 GIMAeGTM F 8 

RETAP*KZ F9 

RETAI— XI*KZ F 10 

ALFAI*-XI*KX p XI 

GO TO (13.'»»3>2), ICASE F 12 

2 DWM-WN*OELX F 13 

IF (PWN.EQ.O.) DWN-DELX F 14 

WN«WN+OWN F 15 

ALFAR*WM+COS (EPSI ) F 16 

BETAR-WN^STNIEPSI ) F 17 

KX-CMPLX(ALPAR»ALFAI) F is 

GO TO 5 F 19 

DEPSI-EPST*DELX F 20 

IF (DEPSI.EO.O. ) DEPSI-DELX F 21 

EPSI«EPSI+OEPSI F 22 

ALFAR*WN*COS (EPSI ) F 23 

BETAR»WN*STN(EPSI ) F 24 

KX«CMPLX( ALFAP, ALFAI) F 25 

GO TO 5 F 26 

3 DBFTP«8ETAR*DELX F 27 

IF (DBETR.FO.O. ) DBETR»OELX F 28 

BETAP«BETAP+DBETR F 29 

GO TO 5 F 30 

4 PBETI-RETAI+DELX F 31 

IF (DRETI .EQ.O. ) DBETI-DELX F 32 

RFTAI«8ETAI+0BETI F 33 

5 KZ»CHPLX(3ETAR,BETAI) F 34 

NTTB«NITB+1 F 35 

INOOX*2 F 36 

RETURN F 37 

6 CONTINUE F 38 

GO TO (13,8,9,7), ICASE F 39 

7 nPHDK«{GIM-GIMA)/DWN F 40 

DWN«-GIM/nPHDK F 41 

IF (A9S(r)WN/WN),LE.EPS) GO TO 11 F 42 

WRITE ( 6 , 16 ) NITB,WN,GIM,DPHDK,OWN F 43 

uj WN«WN+r)WN F 44 

ALFAP -WN*Cn3 ( EPSI ) F 45 



+ T I D Y + 

SDRPniJTINE GIMTR ( G I M , E P S I » E PS I P » I C ASE , NI TB » INOHX » WN» E P S ) 



PFTAR«WM*STN( EPSI ) 

F 

A6 


alfai»-xi+kx 

F 

A7 


petai»-xi*kz 

F 

A8 


KX-CMPLX ( ALFAR» ALFAI ) 

F 

A9 


GO TO 10 

F 

50 


DPHDK* (GTM-GIMA ) /DEPSI 

F 

51 


ofpsi*-gtm/ophdk 

F 

52 


IF (ABS(OFPSI/EPSI).LE.EPS) GO TO 11 

F 

53 


write (6»16) NITB,EPSI*57.2g577*GIM»DPHDK,DEPSI*57. 29577 

F 

5A 


EPSI=EPSI+0FPSI 

F 

55 


ALFAP*WN*COS (EPSI ) 

F 

56 


BFTAP»WN*STM(EPSI ) 

F 

57 


ALFAI»-XT*KX 

F 

58 


BETAI «-XT*KZ 

F 

59 


KX«CMPLX( ALF AR> ALFAI ) 

F 

60 


GO TO 10 

F 

61 

9 

OPHDK» (GIM-GIMA) /DBETI 

F 

62 


OBETI»-GIM/CIPHDK 

F 

63 


IF (ABSCORETI/BETAI ).LE.EPS) GO TO 11 

F 

6A 


WRITE (6,16) NITB,KZ,GIM,DPH0K»0BETI 

F 

65 


BETAR«KZ 

F 

66 


betai*b«:tai+pbeti 

F 

67 


GO TO 10 

F 

68 

9 

OPHDK« (GIM-GIMA ) /DBETP 

F 

69 


DPETP *-GTM/DPHDK 

F 

70 


IF { ABS(DBETR/eETAR ) .LE.EPS ) GO TO 11 

F 

71 


WPITF (6,16) MITB,KZ,GIM»DPH0K,0BETR 

F 

72 


BETAT--XI+KZ 

F 

73 


BFTAP«BETAR+P8ETR 

F 

7A 

10 

KZ«CMPLX ( PETAP, BETAI ) 

F 

75 


GIMA-GIM 

F 

76 


NITB»NITB+1 

F 

77 


IF (NITB.EO. ITR) GO TO 12 

F 

73 


INOOX-2 

F 

79 


PFTURN 

F 

80 

11 

CONTINUE 

F 

81 


TNDPX»10 

F 

82 


GO TO lA 

F 

83 

12 

WPITF (6,15) 

F 

8A 

13 

STOP 

F 

85 

lA 

RFTUPS 

F 

86 

C 


F 

87 

15 

format (1H ,26HEXESSIVE NO. OF ITERATIONS) 

F 

88 

16 

FORMAT (1H,1X,I2,2E13.A,E13.A,AE13.A) 

F 

09 


END 

F 

90- 



♦ T I D Y + 

SUBROUTINE SONG (N»X,YM,Z) 

SUBROUTINE SONG (N>X,YM,Z) G 1 

COMPLEX X(N),Z G 2 

DIMENSION GMKlODf GM2(101), YM(lOl), Y(lOl) G 3 

DO 1 I»1,N G 4 

Y (I )«YM(N+1-I ) G 5 

GM1(I)»PEAL(X(N+1-I)) G 6 

1 GM2H )-AIMAG(X(M+l-I ) ) G 7 

VRI«0. G 8 

VII«0. G 9 

DO 2 I«2*N G 10 

VPI»VRI+(GMl ( D+GMK I-l) )*(Y( I )-Y( I-D) /2. G 11 

VII-VII+(G«2 (I)+GM2(I-1) )♦ (Y(I )-Y(I-D) /2. G 12 

2 CONTINUE G 13 

Z»CMPLX(VRI, VII ) G 14 

RETURN G 15 

END G 16- 


U1 

•*vj 



Ul 

00 


♦ TIDY* 


SUBROUTINE CFA (EPSICF,N) 


SUBROUTINE CFA (EPSICF^N) H 1 

PEAL '^U»MUo,MUPP H 2 

dimension IJ(101)> W(101)» UPPdODf WPP(lOl) H 3 

COMMON /CCC/ Y( 101 )*Ul( 101 )»UP( 101 )»UlPP( 101 )»Wm 01 )»WP( 101 )»WlPP H 4 

l(101),T(101),TP(101),TPP(101),MU(10n»MUP(101)»MUPP{101),DVDT(101) H 5 

2»nVDTP(101 ), PR(lOl) H b 

no 1 J«1,N H 7 

I«N-J+1 H 8 

U(I)«U1(J) H 9 

W(I)»W1(J) H 10 

UPP(I1«U1PP( J) H 11 

WPPd ) »W1PP( J ) H 12 

1 CONTINUE H 13 

Bl-l.E+06 H lA 

DO 2 I»2»N H 15 

IP (I) ( I) .GT.0.999) GO TO 3 H 16 

A1«U(I)/W(T) H 17 

A2«UPP (I )/WPP (I ) H 18 

R2*A1-A2 H 19 

IF ( (B2.GE.0. .AND.Bl.LE.O. ) .OR. (B2.LE.0..AND.B1.GE.0. )) ISAVE-I H 20 

B1«R2 H 21 

2 CONTINUE H 22 

3 CONTINUE H 23 

EPSICF>-ATAN(UnSAVEWW(ISAVE) >+3.1^15977 H 2A 

WRITE (6, A) EPSICF+57. 29577 H 25 

RETURN H 26 

C H 27 

A FORMAT (IH . THE PS IC F ■, El2 , 5 ) H 28 

END H 29- 



+ T I D Y * 

SUBROUTINE PROF ( YARG» UL » UP# UP P» WL » WP » W PP»TL »TP» TPP» MU» MUP» MUPP» AL 


SUBROUTINE PROF ( YARG> UL * UP# UP P, WL » WP» WPP» TL »T P» TPP» MU» MUP» MUPP» AL I 1 

1FA»AI.FAP» PRANDL#KL» INDEX) I 2 

PEAL INTER I 3 

REAL NU»MUP,MUPP#NACH»NO I 

COMPLEX KX»KZ#OMEGA#KXX»KZZ.XI I 5 

COMMON /BBB/ KX» K XX , K Z , KZ Z # DME G A » X I , P , Q, 6AMM A* OB I 6 

COMMON /BBl/ XC#CHL#UFS»R#MACH»ETA,PRA#ND#NS#IE I 7 

COMMON /CCC/ Y(101)»U(101)»DU(101)»DDU(101)»W(101)»0W(101)»0DW(101 I 8 

1) #T{101)#DT(101 )#DDT(101)» VS(101)»VSP(101)»VSPP(101 )»DVDT( 101) »DVD I 9 

?TP(101)#PR(101) I 10 

C I 11 

DO 1 J«KL»IF I 12 

I«J I 13 

IF (YARG.GT.Y(J) ) GO TO 2 I 14 

IF (YARG.EQ. Y(J ) ) GO TO 3 I 15 

1 CONTINUE I 16 

? MIN-I-3 I 17 

IF (I.LE.3) MIN»1 I 18 

IF (I.GE.(IE-2) ) MIN. IE-6 I 19 

UL«INTER( Y»U# YAPG#6# MIN) I 20 

UP’INTER( Y,DU#YARG#6»MIN) I 21 

UPP. INTER (Y# DDU# YARG»6#MIN ) I 22 

Wl«INTER(Y#W,YARG#6»MIN) I 23 

WP.INTER(Y,0W#YARG#6»MIN) I 24 

WPP.INTER(Y#0DW#YARG»6#MIN) I 25 

TL.INTER(YiT,YARG#6,MIN) I 25 

TP. INTER! Y,DT#YARG,6#MIN) I 27 

TPP. INTER (Y, DDT# YARG»6#MIN ) I 28 

PPANDL.PRA I 29 

MU«INTER( Y,VS#YARG#6#MIN) I 30 

MUP. INTER (Y, VSP# YARG#6»MIN ) I 31 

MUPP. INTER (''#VSPP»YARG»6# MIN) I 32 

ALFA.TNTEP(Y,DVDT#YARG#6»MIN) I 33 

ALFAP-INTER! Y#DVDTP# YARG»6#MIN) I 34 

Kl«T I 35 

RETURN I 36 

3 UL.U(I) I 37 

UP.DU(I) I 38 

UPP.OOU(I) I 39 

WL-Wd) I 40 

WP-DWd) I 41 

WPP.DDWd) I 42 

o. TL«Td) I 43 

TP-OTd) I 4^ 

TPP.ODTd) I 45 


» T I 0 Y * 

o SUBRPUTINE PROF ( Y ARG , UL » UP, U P P , W L , W P, W P P, TL » T P, TPP, MU» MUP » MU»P, AL 

PRANOL-PRA r 46 

MU«VS(I) I ^7 

HUP«VSP(I) I 48 

HUPP-VSPPd) I 

ALFA=DVDT(I) I 50 

ALFAP»OVOTP( I ) r 51 

K I - T 15 2 

RETURN I 53 

END I 5A- 



♦ TIDY* 


REAL FUNCTIONINTER( X, Y,XARG»IDEG»MIN) 


PEAL FUNCTiaNINTER(X,Y,XARG»IDEG,rilN) 
DIMENSION X(101)» Y(lOl) 

1 FACTPR-1.0 
MAX-MTN+TOEG 
PD 2 J-MIN,MAX 

IF (XARG.NE.X(J) ) GO TO 2 
INTER. Y(J> 

RETURN 

2 FACTOR«FACTOR*(XARG-X( J) ) 

YFST.0.0 

on A I«MTN,MAX 

TFRM«Y(I)*FACT0R/(XAR6-X(I )) 

DO 3 J«min,MAX 

3 IF (I.NE.J) TERM. TERM/ (X(I)-X(J) ) 

A YEST-TEPM+YEST 

INTER. YEST 

RETURN 

FNp 


J 1 
J 2 
J 3 
J A 
J 5 
J 6 
J 7 
J 8 
J 9 
J 10 
J 11 
J 12 
J 13 
J lA 
J 15 
J 16 
J 17 
J 18- 


O' 



c^ 

to 

* T r n Y * 

SUBROUTINE GVEC (X>C»S) 




SUBROUTINE OVEC 

(X,C.S) 

K 

1 


niNENSION C( 1 )> 

S(1 ) 

K 

2 


RETURN 


K 

3 


END 


K 




♦ TIDY* 


SUBROUTINE COMINV (A>N»D) 


SUBROUTINE CDMINV (A»N/D) L 1 

COMPLEX A(N»N)»D»BIGA>HOLO L 2 

DIMENSION L(8)# M(8) L 3 

0-(1.EO,O.FO) L 4 

no 18 K*1,N L 5 

L(K)«K L 6 

M(K)«K L 7 

BIGA«A(K,K) L 8 

DO ? J-K»N L 9 

no 2 I«K»N L 10 

IF (C ABS ( BIGA)-CABS ( A( I> J ) ) ) 1,2»2 L 11 

1 BIGA«A(I,J) L 12 

L(K)»I L 13 

M(K)«J L lA 

2 CONTINUE L 15 

J-L(K) L 16 

TF (J-K) 5»5»3 L 17 

3 DO A I«1.N L 18 

HOLD«-A(KW) L 19 

A(K,I )-A( J,I ) L 20 

A A(J#I)»HOLn L 21 

5 I«M(K) L 22 

IF (I-K) 8,8i6 L 23 

6 DO 7 J»ltN L 2A 

HOLO»-A(J»K) L 25 

A( J,K)«A( Jtl ) L 26 

7 A(J»I)*HnLn L 27 

3 IF (CABS(BIGA)) 10>9»10 L 28 

9 D-(O.EO»O.FO) L 29 

RETURN L 30 

10 DO 12 I-1,N L 31 

IF n-K) 11,12^11 L 32 

11 A(I,K)«A{I,K)/(-BIGA) L 33 

12 CONTINUE L 3A 

DO 15 I»1,N L 35 

HOLD»A(I,K» L 36 

DO 15 J*1»N L 37 

IF (I-K) 13»15»13 L 38 

13 IF (J-K) 1A,15»1A L 39 

lA A(I,J)»HOLD*A(K,J)+A(I,J) L AO 

15 CONTINUE L A1 

DO 17 J*1»N L A2 

IF (J-K) 16,17fl6 L A3 

16 A (K» J ) -A(K» J ) /BIGA L AA 

17 CONTINUE L A5 


o^ 

4 > 


♦ TIDY* 


SUBROUTINE CDMINV (A,N,D) 



[)«D*BIGA 

L 

46 


A(K»K)«1.F0/BIGA 

L 

47 

18 

CONTINUE 

L 

48 


K«N 

L 

49 

19 

K« (K-1 ) 

L 

50 


IF (K) 26»26,20 

L 

51 

20 

I «L (K ) 

L 

52 


IF (I-K) 23.23^21 

L 

53 

21 

DO 22 J«1,N 

L 

54 


H01D«A(J>K) 

L 

55 


A ( J»K )»-A ( J» I ) 

L 

56 

22 

A ( J» I ) »HOLD 

L 

57 

23 

J ■M(K ) 

L 

58 


IF (J-K) 19,19,24 

L 

59 

24 

DO 25 I«1,N 

L 

60 


HOLD» A(K, I) 

L 

61 


A(K,I )— A( J, I ) 

L 

62 

25 

A( J,n*MOLD 

L 

63 


GO TO 19 

L 

64 

26 

RETURN 

L 

65 


END 

L 

66- 


* T I D Y * 

SUBROUTINE GRVEL ( ALAN, W2 » W3> W4 » W6» W8» TATA, E PGR» SM ) 

SUBROUTINE GRVEL ( AL AM , W2, W3, W4 , W6, W8,T AT A, E PGR, SM) M 1 

PEAL M ACH,MU,MUP,MUPP,ND M 2 

rOMPLEX Z1{101),Z2(101),Z3(101),ZA(101),Z5(101),Z6(101>,Z7(101),Z8 M 3 

1(101) M <» 

COMPLEX W2(101),W3(101),W4(101),W6(101),W8(101) M 5 

COMPLEX G1(101),G2( 101) M 6 

COMPLEX GZ1,GZ2,G23,G24,GZ5,GZ6,GZ7,G28,GW2,GW3,GW4,GW6,GW8 M 7 

COMPLEX VIGl , VIG2, VIGIO, VIG20,TATA,VIHX, VIHZ M 8 

COMPLEX XI,KX,KXX,KZ,KZZ,A111, A15,A102, AlO M 9 

COMPLEX A12,G,ALAM, OMEGA M 10 

COMMON /BBB/ KX , K XX, K2 , KZ 2 , OME G A , X I , P,Q , G AMM A, OB M 11 

COMMON /BBl/ XC,CHL,UFS,R,MACH,ETA,PRA,ND,NS,N M 12 

COMMON /CCC/ Y(101),U(101),UP(101),UPP( 101),W(101),y<P(101),WPP(101 M 13 

1),T(101),TP(101),TPP(101),MU(101),MUP(101),MUPP(101),ALFA(101),ALF M 14 

PAPdOl ),PR (101) M 15 

COMMON /FFF/ 21, 22, Z3, Z4, Z5, 26, 27, 28 M 16 

: M 17 

DO 1 J»l,N M 18 

I»J M 19 

A111«KX*U(I)+K2*W(I )-OMEGA M 20 

A15«-KXX-KZ2+Q*MUP( I )+TP(I )/( MU( I )*T(I) )+QyTPP(I ) /T{ I ) -XI*R*4lll M M 21 

lMU(r)*T(D) M 22 

AO.MUP(I) /MU( I)+TP(I )/T( I ) M 23 

A102»KX+UP( T )+KZ+WP ( I ) M 24 

A10«A9*A111+A102 M 25 

A121«Q/T(I)+ALFA(I)/MU(I) M 26 

A12-A121*A102+0*MUP(r)+Alll/(MU(I)+T{n ) M 27 

A13*2.*MUP(I)/MU(I)+0*TP(I)/T(I) M 28 

AP»PP ( D’MGAMMA-l. )*MACH*MACH M 29 

G«P/MU(I )+XT*0*0B*Alll M 30 

GZ1«(XI*0+IJ(I )/(T(I)’CMU(I) )+2.fKX)*Zl(J) M 31 

GZ3*-XI*(MUP(I)/MU(I)+PYTP(I)/T(I))*23( J) M 32 

G2A.( XI+P/MU ( I)-P*QB*(2.*KX*U( I )+KZ*W(I )-OMEGA ) )*Z4( J ) M 33 

675«P*(2.*KX*U(I)+KZ*W{I)-0MEGA)*Z5(J)/T(I) M 34 

GW2-(G21 + G73 + GZ4 + GZ5)*W2( J ) M 35 

G71«-XT*71(J ) M 36 

GZ4»-XT*0RtU(I)*Z4( J) M 37 

GZ5*XI*U(I)*Z5(J)/T(I) M 38 

GW3«(GZ1+GZ4+GZ5)*W3( J ) M 39 

G71«-XI*(1.-XI*KX*QB*Q+U(I)/G)*A13+Z1(J)/G M 40 

G7?=XIt(XI*KX*CB+0*U(I)/G-l.)*Z2( J)/G M 41 

GZ3« (-?.*RX-XI’4>P*U( I ) / (MUd )*T ( I ) ) -X 1*0 *QB*U ( I )* A15 /G ) +73 U ) /G M 42 

g; GZ4«-XI+0*0R+(A9+U(I)+UP(I)-XI+Q+QB+U(I )*A10/G)+Z4(J)/G M 43 

GZ5*(XI+{ A121+UP(I) + 0*MUP( I )+U(I) / (MU(I )+T{D) ) /G + 0*0B*U( I) +A12/ (G M 44 

1*G))+Z5(J) M 45 



nil nil 


0 ^ 

a^ 


♦ TIDY* 


SUBROUTINE GRVEL ( AL AM , W2» W3> W4 » W6» W 8, TATA» EPGR» SM ) 


G7f)»0 + U(n*(yi + Q*CP + Alll/G)*Z6(J )/(G*T( I) ) M A6 

GZ7»-A13t'K7*Q*0B*U( I )*Z7( J )/ (G*G ) M ^7 

GZe»-KZ*0+0B+U(I)*Z8( J)/(G+G) M 48 

GW4«(GZ1+GZ2+GZ3+GZ4+GZ5+GZ6+GZ7+GZ0)*W4(J) M 49 

GZ3«-2.*XT*A8*UP(I)*Z3(J ) M 50 

GZ4»-XI*R*AP*U(I )*Z4( J)/MU(I ) M 51 

GZ5»(X I + P + PR m>fU( I )/(MU( I )*T ( I ) )+2.*KX )*Z5( J) M 52 

GW6«(GZ3+G74+GZ5)*W6( J ) M 53 

GZ4»-P*K7*9B*U{I )*Z4( J) M 54 

G75«P*KZ*U( I )*Z5U ) /T( I) M 55 

GZ7« I XT*P*U( I )/ IMUl I)*T(1) )+2.*KX)tZ7IJ ) M 56 

GWe» ( GZ4+GZ5+GZ7) *W8( J ) M 57 

G1(J)*GW2+GW3+GVI4 + GW6 + GW8 M 58 

GZ1« ( XT*R*W( I ) / (MU( I )*T( I ) )+2 ,+KZ )*Z1( J ) M 59 

GZ<,,-p*oR*Ky*W(I )*Z4( J) M 60 

GZ5»P*KX*W(I )*Z5( J)/T(I) M 61 

GW2-( GZ1+GZ4+GZ5 l*W2(J ) M 62 

GZ4«-XI*aB*W(I)*Z4( J) M 63 

GZ5*XI*W(I)*Z5U)/T(I) M 64 

' GZ7«-XI*Z7( J ) M 65 

GW3» (GZ4 + GZ5 + GZ7)*W3 ( J > M 66 

GZ1«-0*0B*W( I )*KXt«Al3*Zl( J )/IG*G) M 67 

GZ2«-0*QB*KX*w( I)*Z2(J)/(G*G) M68 

‘ GZ3»(-2.*K7-XI*R*W( I ) /(MU( I )*T( I n-XI*Q*OB*W(I )*A157G)*Z3( J) 7G M 69 

G 74»-X I*Q*O0f ( A9*W ( I )+WP ( I )-X I*0*0B*W(I )*A10/G )»Z4( J ) /G M 70 

GZ5«(XIt( A121*WP{1)+Q*MUP (I )4W(I ) / (MU(1 )*T(1 )) )/G+O4‘0B*Wn)*A12M6 M 71 

1*G))*Z5(J) M 72 

GZ6»0 + wm + ( XI + 0*QB*A111/G)*Z6U) / (G4T( I) ) M 73 

= GZ7»-A13+(K74Q*QB*WI I) /(Gt=G) + XI/G)*Z7(J ) M 74 

GZB— (KZ*0*0B»W(I)/G + XI)*Z8(J)/G M 75 

: GV4» (GZ1+GZ2+GZ3+GZ4+GZ5+GZ6+GZ7+GZ8 )4M4( J) M 76 

GZ3»- 2.*XI*A8*WP(I )*Z3(J ) M 77 

s GZ4 — XI*R*AB*W( I)4Z4( J )/MU(I ) M 78 

GZ5»(XT*RtPP(I)*W(I)/(MU(I)*T(I) )+2.*KZ)*Z5(JI M 79 

= GW6»(GZ3+GZ4+GZ5)*W6(J) M 80 

GZ3«-XI*(MUO(I)/MU(I) + P*TPm/T(n) + Z3(J) M 81 

GZ4.(XI*P/MIJ(I)-P40B*(KZ*W(I)+A111))*Z4(J) M 82 

GZ5*P*(KZ*W(I)+A111)+Z5(J)/T(I) M 83 

GZ7«(XI*R+W(I)/(MU(I)*T(I) )+2.*KZ)*Z7(J) M 84 

GW8»<GZ3 + GZ4 + 6Z5 + GZ7)*W8( J) M 85 

G2(J) -GW2+GW3+GW4+GW6+GW8 M 86 

1 CONTINUE M 87 

: M 88 

CALL SONG (M,G1#Y, VIGl) M 89 

CALL SONG (N,G2>Y,VIG2) M 90 

i 


U* IN) 


♦ TIDY* 

SU'^ROUTINE GRVEL ( AL AM , W2 » W3# W<t # W6» W8» TATA , EPGR> SM ) 

Alll=KX+NO*KZ-OMEGA M 91 

A1^*-KXX-KZ7-XI*R*A111 M 92 

A1?1«0+ALFA(1) M 93 

AP-PR ( l)*(r,AMMA-l. )*MACH*MACH M 94 

6«P+X I*Q*OR*Alll M 95 

GW2«( (XT*R+2.*KX)*Z1(1)+(XI*R-P*QB*(2.*KX+ND*KZ-0MEGA) )*Z4(1)+P*(2 M 96 

l.*KX+N0*K7-nHEGA)*Z5(l) )*W2(1) M 97 

GW3»XI*(-Z1( 1 )-QB*Z4(l )+Z5 ( 1) )*W3( 1 ) M 98 

GW4«(X I*(XI*KX*QB*Q/G-1. )*Z2( 1 ) /G+ (-2. *K X-X I *R-XI*0*0B *A 15 /G ) *Z 3 ( 1 M 99 

D/G+O* (XT+0*QB*A111/G)*Z6(1)/G-KZ*0*QB*Z8(1 )/(G*G) )*W4(1) M 100 

GW6«(-XI*R*A8*Z4(1)+(XI*R*PP(1)+2.*KX)*Z5(1) )*W6(1) M 101 

GW8«{-p*KZ*0B*Z4(1)+P+KZ*Z5(1)*(XI*R+2.*KX)*Z7(1))*W8(1) M 102 

VIGin*-(3W2+GW3+GW4+GW6+GW8)/(?.*ALAM) M 103 

GW2«( (XI*R*ND + 2.*KZ )*Z1(1 )-P*QB*ND*KX*Z4(l) + P*KX*ND*Z5(l) )*W2m M 104 

GW3»(-XI*OB*NO*Z4 (1 )+XI*ND*Z5{ 1 )-XI*Z7( 1 ) )*W3( 1) M 105 

GW4«(-KX*0*QB*ND*Z2( 1) / (G*G)+(-2.*KZ-XI*R*N0-XI*Q*0B*N0*A15ZG)*Z3( M 106 

11 )/G+0+ND*(XI+0*0B*Alll/G)*Z6(l)/G-(KZ*0*QB*ND/G*XI )*Z3(l)/S)*W4(l M 107 

2) M 108 

GW6«(-XI*R*A8*ND*Z4a)+(XI*R*PR(l)*ND + 2.*KZ)*Z5(l))*W6(l) M 109 

GWB»( (XT*R-P*KZ*QB*ND-P*QB*A111)*Z4(1)*P*(KZ*ND+A111 )*Z5(1 )+( XI*R* M 110 

1ND+?.*KZ)*Z7(1) )*W8(1) Mill 

VIG?D--(GW2+GW3+GW4+GW6+GW8)/(2.*ALAM) M 112 

VIG1«VIG1+VIG10 M 113 

VIG2»VIG2+7IG2D M 114 

TATA-VIG2/VIG1 M 115 

SM.TATA M 116 

EPGR*ATAN(SM)*57. 29577 M 117 

WRITE (6»2) M 118 

WRITE (6^3) VIG1^VIG2,TATA,EPGR H 119 

RETURN M 120 

'• M 121 

FORMAT (IH ,/>2X»43HINFDRMATI0NS FROM 6RVEL VIG1» VIG2» TATA,EPGR/ ) M 122 

FORMAT (IH ,5X,1C(011.4,1X),F5.2//) M 123 

END M 124- 


tr> 


APPENDIX I 


THE MEAN FLOW 


The meanflow solution is an input to HADY-I. The boundary layer solution 
is calculated using as input the airfoil geometry, pressure coefficient dis- 
tribution, and suction requirements. 

The boundary layer program used here was adapted from the program of 

4 

Kaups and Cebeci for laminar, compressible boundary layers with adiabatic 
wall and wall suction boundary conditions. Extensive modifications and 

4 

additions was necessary to this program , to suit the need of HADY-I stability 
program. 

A. Input/Output Files 
The program card is 

PROGRAM MFLOW (INPUT, OUTPUT, TAPES = INPUT, TAPES = OUTPUT, TAPED, 

TAPE 10, TAPEll) 

T/PED A file for internal use 

TAPE 10 An output file that contains boundary layer profiles, used 
for parallel stability calculations 
TAPEll An output file that is used with TAPE 10 for nonparallel 
stability calculations 

B. Control Cards 

The following control cards can be used to execute the program 

JOBS, Tt, CM. 

USER, USERNO, PASSWRD. 

CHARGE, CHARNO, LRC. 

GET, MFLOW. 

ATTACH (FT^ML IB/UN f LIBRARY) 

NOEXIT. 

LDSET (LIB = FTNMLIB, PRESET = ZERO) 

MFLOW. 

REPLACE, TAPEIO = TAPEN. 

EXIT. 

7/8/9 end of record 
Input Cards 
S/7/8/9 


68 



C. Program Input 


The input is through data cards 


Card 1 

8A10 


TITLE 

Description of the case 


Card 2 

411, 31, 2X, 2F15.1 


IPRPJT 

Define type of output printed 
= 1 Long print 
= 2 Short print 


IPANPA 

Define type of analysis desired 
= 1 Parallel flow 
= 2 Nonparallel flow 


MK 

Number of input stations before minimum x/c (.not including 
the minimum) 

NK 

Number of input stations before x/c = .001 (not including 
x/c = .001) 

NM 

Number of output stations starting from x/c = .001 

(NM = 11 

RCR 

Streamwise chord Reynolds number where input is given 
(tunnel conditions and suction distribution) 

RCU 

Streamwise chord Reynolds number where solution is 

needed 

Card 3 

213, 3F10.0 


NI 

Number of input stations for the streamwise airfoil 


NZT 

Number of input stations where suction coefficient 
speci fied 

Cs is 

ETAE 

Estimated value of maximum n at the first station 


DETA1 

First An - step size 


VGP 

Variable grid parameter 


Card 4 

8F10.0 


X 

Chord length in feet for the streamwise airfoil 


SWLE 

Leading-edge sweep in degrees 


SWTE 

Trai ling-edge sweep in degrees 



CMACH Freestream Mach number 

UREF Freestream velocity in ft/sec (only if = 0) 

2 

TPRES Freestream static pressure, in Ib/ft 

TT Freestream static temperature, in degrees Rankin 

PR Prandtl number 

Card 5 2F15.10 

XLE >./c of the leading edge 

YLE y/c of the leading edge 

Card 6 4F15.10 

Total of NI cards, one per station. These cards contain informations 
of the normal to the leading edge airfoil, (airfoil data output from 
Garabedian program is used directly). 


XA x/'C yal ue 

YA y/c value 

CM local Mach number 

CP 2 

Card 7 8F10.0 


Input Cp values (two-dimensional) 


BLP Input suction coefficient values, defined as C,, 
Total of NZT points 


(-pV)^/(pU)^ 


D. The Program MFLOW 
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+ T I D Y * 

PROGRAM MFLOW( INPUTfOUTPUT, TAPE5«INPUT»T4PP6-0UTPUT,TAPE9»TAPE10>T 


PROGRAM MFLOW(INPUT»OUTPUT,TAPE5»rNPUT,TAPE6-OI)TPUT»TAPE9,TAPElO,T A 1 

lAPEll) A ? 

COMMON /BLCO/ NZT» NZ> NP» I T» X, ROFS» C MACH» TT, ET A ( 1 51) » DETA ( 1 51 ) » A ( 15 A 3 

11) » Y(151 ),IPANPA A A 

COMMON /BLCl/ HE, PR » C MUF S> UF S» CEL ( 1 51 ) » BETA 1 ( 151 ) , UF ( 151 ) , WE( 151 ) , A 5 

1Z(151),PE(151),PHI(151),RH0E( 1 5 1 ) , XC ( 1 5 1 ) , CMUE ( 1 5 1 ) , PI (151),P3(151 A 6 

2),PA(151),RR(151),BLP(151),DDW{151) A 7 

COMMON /PROF/ DELV ( 1 51 ) , F ( 101, 2 ) , U ( 101 , 2 ) ,V ( 101, 2 ) , G ( 101, 2 ) , W ( 1 01, A 8 

12) ,T(101,2),Ba01,2),C(151),BG(101,2),E(10],2),DENR(101, 2),CA1( 101 A 9 

2,2),CA2(101,2) A 10 

COMMON /PAR/ Al, A2, A3,VGP,NNN,IPRINT A 11 

COMMON /BAB/ H A 12 

C A 13 

CALL INTIAL A lA 

NZ«1 A 15 

IS0LV2-0 A 1(S 

ITMAX-20 A 17 

1 IGROW-0 A 18 

2 IT«0 A 19 

3 IT-IT+1 A 20 

IF (IT.LE.ITMAX ) GO TO A A 21 

WRITE (6,10) A 22 

60 TO 9 A 23 

A IF (IS0LV2.EQ.1) CALL FLUID A 2A 

IF (H.LT.0.0) GO TO 9 A 25 

CALL COEF A 26 

CALL S0LV6 A 27 

IF (ABS( DELVd ) ).LE. 0.0001) GOTO 5 A 28 

IF (IS0LV2.EQ.1) CALL S0LV2 A 29 

GO TO 3 A 30 

5 IF (IS0LV2.EQ.0) 60 TO 6 A 31 

CALL S0LV2 A 32 

GO TO 7 A 33 

6 IF (CMACH.EQ.0.0) 60 TO 7 A 3A 

IS0LV2-1 A 35 

GO TO 2 A 36 

7 IF (ABS(T(NP,2) ).LE.l.E-0) GO TO 8 A 37 

IF (NP.EQ.lOl) GO TO 8 A 38 

IGROW-IGROW+1 A 39 

IF (I6R0W.GT.1) GO TO 8 A AO 

LL*1 A Al 

CALL PROFTL (LL) A A2 

GO TO 2 A A3 

0 CALL OUTPUT A AA 

IF (NZ.LE.NZT) GO TO 1 A A5 


♦ TIDY* 

PROGRAM MFLQWl INPUT, OUT PU T, T APE 5 « I NPUT , T APE 6-OUT PUT, T A PE9, TAP El 0, T 


<J CONTINUE 

IF (IPANPA.EQ. 1) STOP 
CALL XZDER 
STOP 
C 
C 

10 FORMAT (1H0,23HITERATI0NS EXCEED ITMAX) 

END 


A ^6 
A 

A <t8 
A 49 
A 50 
A 51 
A 52 
A 53- 



* T I D Y * 


SUBROUTINE INTUL 


SUBROUTINE INTIAL B 1 

COMMON /BLCO/ NZT, NZ» NP, I T^ X. ROF S* C MACH, TT, ETA ( 151 » > OETA ( 151) , A ( 15 B 2 

11) ,Y(151),IPANPA B 3 

COMMON /BLCIZ HE, PR, C MUF S, UF S, C EL ( 1 51 ) , BETA 1 ( 151 ) , UE ( I 51) , WE ( 151 ) , B A 

1Z(151),PE(151),PHI(151),RH0E(151),XC(151),CMUE{151),P1(151),P3(151 B 5 

2),P4(151 ),RP(15l),BLP(151),OOW(151 ) B 6 

COMMON /PROF/ D EL V ( 1 51 ) , F ( 10 1, 2 ) , U { 101, 2 ) , V ( 1 01, 2 ) , G( 101, 2 ) , W ( 101, B 7 

12) ,T(101,2) ,B(101,2),C(151), BG( 101,2 ),E( 101,2), DENR (101,2 ),CA1( 101 B 8 

2,2 ),CA2( 101,2) B 9 

COMMON /PAR/ A1,A2, A3,VGP,NNN,IPRINT B 10 

DIMENSION DUE(151), DWE(151), DPR(15l), TITLE(8) B 11 

DIMENSION XA(151), YA(151), CM(151), CP2(151), CP3(151) B 12 

DIMENSION DUMK151), DUM2(151) B 13 

DIMENSION SARA(151,3), SARBdSl), SARBN(151), SAR A0( 151, 3 ) , WK(200 B lA 

10), Y0(3) B 15 

C 8 16 

READ (5,35) TITLE B 17 

READ (5, A3) I PR INT , IP AN P A, MK, NK, NM, RC R, RC U B 18 

READ (5,36) NI,NZT,ETAE, DETA1,VGP B 19 

READ (5,37) X , S WLE, S WTE, C MAC H, UR E F, TPRE S , TT, PR 8 20 

READ (5,3A) XLE,YLE B 21 

READ (5,3A) ( X A ( I ) , YA ( I ) , C M ( I ) , C P 2 ( I ) , I -1 , NT ) B 22 

READ (5,37) ( BL P( I ) , I»l, NZ T ) B 23 

SARBN(1)«0.0010 B 2A 

DO 1 I«2,10 B 25 

1 SAR8N(I)«SARBN(r-l)+,0010 B 26 

DO 2 I«ll,28 B 27 

2 SARBN(I)«SARBN(I-l)+.005 B 28 

DO 3 I«29,NM B 29 

3 SARBN(I)«SARBN(I-1 )+.01 B 30 

MKl-MK+1 B 31 

DO A I»MK1,NI B 3? 

M«I-MK1+1 B 33 

SARA(M,1)»YA(I ) B 3A 

SARA(M,2)-CP2( I) B 35 

A SARB(M)-XA(I ) B 36 

MKK«MK-2 B 37 

DO 5 I«MKK,NZT B 38 

M-I-MKK+1 B 39 

5 SARA(M»3)-BLP(I) B AO 

IW«-1 B A1 

N«NI-MK1+1 B A2 

DO 7 I«1,NM B A3 

XO»SARBN(I) B AA 

CALL lUNI (151,N,SARB,3,SARA,2,X0,Y0,IW,IERR) B A5 



4 > 


* T I D Y ♦ 


SUBROUTINE INTIAL 



IF (lERR.EQ.O) GO TO 6 

B 

A6 


WRITE ( 6 , 33 ) lERR 

B 

A7 


STOP 

B 

A0 

6 

SARAOd, 1)»Y0( 1) 

B 

A9 


SARA0(I»2)«Y0(2) 

B 

50 


SARAQ(I,3)-YO(3) 

B 

51 

7 

CONTINUE 

B 

52 


NKl-NK+1 

B 

53 


DO 8 I«1^NM 

B 

5A 


M»NK1+I-1 

B 

55 


XA(M)-SAR0N (I) 

B 

56 


YA{M)»SARAO(Id) 

B 

57 

0 

CP2(M)»SARA0(I*2) 

B 

58 


DO 9 I»1»NM 

B 

59 


M«MK1+I-1 

B 

60 

9 

BLP (M)-SARAO(I ,3) 

B 

61 


NI»NM+NK 

B 

62 


NZT-NI-3 

B 

63 


PI»4.*ATAN( 1. ) 

B 

6A 


CD-COS(PI*SWLE/10O. ) 

B 

65 


DO 10 I«1»NI 

B 

66 


XA ( I )«XA (I )-XL E 

B 

67 


YA<r)»(YA(I)-YLE)*CD 

B 

68 

10 

CP3(I)»CP2{I)*CD*CD 

B 

69 


DO 11 I»1,NI 

B 

70 


A( I)»XA( I) 

B 

71 

11 

Y(I )-YA(I ) 

B 

72 


DO 12 I»1,NZT 

B 

73 


XC(I)»XA(I+3) 

B 

7A 

12 

PA (I )»CP3(I+3) 

B 

75 


SORCR»SORT( RCR ) 

B 

76 


SORCIJ«SORT(RCU) 

B 

77 


DO 13 I-1»NZT 

B 

78 

13 

BLP(I)«BLP(I)*SORCR/SQPCU 

B 

79 


TPRES-TPRES*RCU/RCR 

B 

80 

C 


B 

81 


WRITE (6#39) TITLE 

B 

82 


WRITE (6»A1) { I»A (I ) , Y( I ) ,I»1,NI ) 

B 

83 

C 


B 

8A 


XINPUT-X 

B 

85 


CMSO-CMACHf+Z 

B 

86 


ROFS«TPRES/ (1716.+TT) 

B 

87 


UFS«CMACH*SQRT( 1.A+1716.+TT) 

B 

88 


IF (CNACH.EO.0.0) UFS-URFF 

B 

89 


CMUFS-2.27E-08+TT++1, 5/( TT+198.6) 

B 

90 



♦ TIDY* 


SUBROUTINE INTIU 



TTT»TT*(1.0+0. 2*CMS0) 

B 

91 


HE«TTT*6006.0 

B 

92 


RE Y«UFS*ROFS*X INPUT /CMUFS 

B 

93 


Al-l.+VGP 

B 

94 


A2» A1 + V6P** 2 

B 

95 


A3»A2+VGP**3 

B 

96 

c 


6 

97 


DEL«AC0S<1.0-A<1)) 

B 

9B 


ETA( 1)«DEL 

B 

99 


IF (A(l) .GT.A(2)) ETA(D— DEL 

B 

100 


DO 16 I«2»NI 

B 

101 


PHANG«ACQSa.O-A (T ) ) 

B 

102 


IF (A(I).LT.A(I-D) GO TO 14 

B 

103 


ETAT I )«PHANG 

B 

104 


GO TO 15 

B 

105 

lA 

ETA (I )— PHAN6 

B 

106 

15 

IF (AID.EO.O.O) ETA(I)«0.0 

B 

107 

16 

CONTINUE 

B 

108 


CALL SPLINE ( Y, ETA, NI, DELV » 

B 

109 


TLE»TAN(0.0174533*SWLE) 

B 

110 


TTE»TAN{0.0174533*SWTE) 

B 

111 


CB-TLE-TTE 

B 

112 


X»X*SQRT (1. 0+TLE**2) /CB 

B 

113 


DO 17 I»1,NI 

B 

114 


SF»SIN(ETA( I) ) 

B 

115 


TCS«TLE-CB*A{I ) 

B 

116 


FF»1,0+(CB*Y(I ) )**2+TCS*TCS 

B 

117 


DF-2.0*CB*(-TCS*SF+CB*Y( I )*DELV( I ) ) 

B 

118 


XF»(CB*SF+TCS*0F/FF/2.0)**2 

B 

119 


YE- 0,25*(DF/FF )**2 

B 

120 


ZF»(CB*(DELV(I )-Y(I )*DF/FF/2.0) )**2 

B 

121 


DETA(I)»SORT( (XF+YF+ZF)/FF) 

B 

122 

17 

CONTINUE 

B 

123 


CALL INTEG ( ETA, DET A, C, NI ) 

B 

124 

C 


B 

125 


DO 23 I»1,NZT 

B 

126 


IF. (CMACH.EQ.O.O) GO TO 18 

B 

127 


DPR{ I)»1.0+0.7*P4(I )*CMSQ 

B 

128 


PE ( I)-1.0+(1.0-DPR( I) **0,285714 )/ (0.2*CNS0) 

B 

129 


GO TO 19 

B 

130 

18 

PE( I)»1.0-P4(I ) 

B 

131 


DPR ( I)«l. 0+UREF*UREF*P4( I )/TT/343 2.0 

B 

132 

19 

IF (I.GT.l) GO TO 20 

B 

133 


DEL-ACOS (l.O-XC(I) ) 

B 

134 


P3(l )«DEL 

B 

135 


ON 


♦ TIDY* 


Sl>BROUTINE INT14L 


IF (XC(1).GT.XC(2) ) P3(1)«-DEL B 136 

GG TO 23 B 137 

20 PHANG»ACOS( 1.0-XCm ) B 138 

IF (XC(I).LT.XCn-l)) GO TO 21 B 139 

P3(I)«PHANG B 1<»0 

GO TO 22 B lAl 

21 P3(I)«-PHANG B 142 

22 IF (XC(I).EQ.O.O) P3(I)-0.0 B 143 

23 CONTINUE B 144 

CALL CUBIC (C»ETA»NI,P3>NZT* Z ) B 145 

0'JE(1)-S0RT (PE (D) B 146 

DWE(1)»0.0 B 147 

NUM«1 B 148 

DZ1»Z(1) B 149 

Z( 1)«0.0 B 150 

CR«CB*CB/ (1.0+TLE++2) B 151 

CC-2.0+TLE/CB B 152 

RR ( 1) «0. 0 B 153 

DO 24 I«2»NZT B 154 

Z(I )-Z(I )-DZl B 155 

DZ-Z(I)-Z(I-1) B 156 

RR(I)»RR(I-1)+DZ*X B 157 

G1--DWE(I-1 )*DZ B 158 

P1(1)»Z(I-1)+0.5*DZ B 159 

CALL CUBIC (PE>Z»NZT,P1,NUM,UE) B 160 

G2TRM— (DUE(I-l)+Gl/2.0)* + 2+UE(l) B 161 

IF (G2TRM.LT.0.0) G2TRM-0.0 B 162 

G2«-SQRT( G2TRM )*DZ B 163 

G3TRM«-(DUE(I-H+G2/2.0)**2+UE( 1) B 164 

IF {G3TRM.lt. 0.0) G3TRM-0.0 B 165 

G3«-S0RT(G3TRM )*DZ B 166 

G4TRM-PE (I )-(DUE(I-l)+G3)**2 B 167 

IF ( G4TRM.L T. 0. 0) G4TPM«0.0 B 168 

G4»-S0RT(G4TRM )*DZ B 169 

DUE( I)-DUE(I-l)+(Gl+2.0*G2+2.0+63+G4)/6.0 B 170 

ALA«PE(I )-DUE (I )**2 B 171 

IF (ALA. LT. 0.0) ALA-0.0 B 172 

DWEd)-SQRT(ALA) B 173 

24 CONTINUE B 174 

CALL SPLINE ( 0 W E» Z> NZ.T, 00 W) B 175 

DOW(D— 2.0*(DUE(2)-DUF(1))/Z(2)/Z(2) B 176 

C B 177 

DO 25 I»1,NZT B 178 

DUMl (I )»OUE (I )+UFS B 179 

DUM2(I)-DWE (D + UFS B 180 



* T I 0 Y * 


SUBROUTINF INTIAL 


25 CONTINUE B 181 

WRITE (6,40) CMACH,UFS,TPRES,TT,PR,ROFS,CMUFS,REY,XINPUT,X,SWL£,SW B 182 

1TE,NI,NZT,ETAE,DETA1,VGP B 183 

WRITE (6,42) (I,XC(I),Z(I),RR(I),P4(I),BLP(n,DUE(n,DWE(I),DDW(I) B 184 

1,DPR(I),I-1,NZT) B 185 

UFS2«UFS* + 2 B 186 

DO 28 J»1,NZT B 187 

UE( J) -UFS*OUE ( J ) B 138 

WE ( J )»UFS*DWE( J ) B 189 

BETA1(J)«DWE(J)/DUE(J) B 190 

PE( J)«OPR(J)*TPRES B 191 

IF (CMACH.EQ.O, 0) GO TO 26 B 192 

TE-TT+(1.0-0.2*CMSQ*(DUE( J)++2+DWE(J)**2-1.0)) B 193 

S«(-DWE( J)*(DDW(J)-DUE( J) ) * ( UFS2 / ( 171 6. *TE ) ) * ( 1. 0+ ( 198 . 6-TF ) / ( 7 . 0* B 194 

1(198. 6+TE)) )) B 195 

RHOE(J)«PE( J)/(1716.*TE) B 196 

GO TO 27 B 197 

26 TE»TT B 198 

S-0.0 B 199 

RHOE(J)«ROFS B 200 

27 CIUE (J )«2.27E-08*(TE**1 .5/( TE+198.6 ) ) B 201 

Pl( J)-D0W(J)/DUE(J) B 202 

P4( J)»BETA1(J)**2 B 203 

P3( J)«0.5*(2.0*DDW(J)/DUE(J)+P4(J)+S*BETA1(J)) B 204 

BLP( J)»SQRT(UE( J)*RHOE( J)*X/CMUE( J ) ) *BL P( J ) *UFS*ROFS /RHOE ( J ) / UE ( J ) B 205 

IF (J.EQ.l) GO TO 28 B 206 

BETA1B»0.5*(BETA1( J)+BETA1( J-1) ) B 207 

CEL(J)-BETA1B/(Z(J)-Z(J-D) B208 

28 CONTINUE B 209 

CEL(1)-0.0 B 210 

C B 211 

DETA (1 )-DETAl B 212 

ETA(1)«0.0 B 213 

IF ((VGP-l.O).LE. 0.001) GO TO 29 B 214 

NP-ALOG( (ETAE/DETAd) )*(VGP-1.0)+1.0)/AL0G(VGP) + 1.001 B 215 

GO TO 30 B 216 

29 NP«ETAE/DETA( 1 ) +1.001 B 217 

30 IF (NP.LE.lOl) GO TO 31 B 218 

WRITE (6,38) B 219 

NP-101 B 220 

31 DO 32 J-2,101 B 221 

DETA(J)»DETA(J-1)*VGP B 222 

ETA( J)«ETA(J-1)+DETA(J-1) B 223 

32 A( J)»0.5*DETA(J-1) B 224 

LL«0 B 225 


CO 


♦ TIDY* 


SUBROUTINE INTIAt 


CALL PROFIL (LL) B 226 

RETURN B 227 

C B 228 

C B 229 

C B 230 

33 FORMAT ( 1H1» 10X» IIHIERR NE 0 »5X»I5) B 231 

34 FORMAT (4F15.10) B 232 

35 FORMAT (8A10) B 233 

36 FORMAT (2I3»3F10.0) B 234 

37 FORMAT (8F10,0) B 235 

38 FORMAT (1H0»36HNP EXCEEDED DIMENSIONS -- SET TO 101) B 236 

39 FORMAT (IHI^SAIO) B 237 

40 FORMAT (IHOWHMACHN E 14 . 6, 3 X, 7HUF S E14. 6» 3X, 7HPF S «,E14.6» B 238 

13X,7HTFS »,E14.6,3X,7HPR », E14 .6/lH0,7HRQFS « » E 1 4 . 6 » 3X, 7HMUF B 239 

2S «#E14.6»3X»7HREC E14. 6/lH0» 7HCH0R0 F14, 6» 3X>7HRADIUS»> E14 B 240 

3.6»3X,7HLESW - , € 14 . 6» 3X , 7HTE SW E14. 6/ IHO, 7HNI - , I 3 » 14X , 7HN Z B 241 

4 «»I3#14X,7HETAE E14,6» 3X,7H0ETA1 « »E 14. 6, 3X » 7HVGP »»E14.6 B 242 

5/) B 243 

41 FORMAT ( //1H0MX^30HSTREAMWI SE AIRFOIL COOROINATES/lHO, 3H NI,8X,3H B 244 

IX/C, 16X, 3HZ/C/ ( IH »I3#3X,E14.6,5X,E14.6)) B 245 

42 FORMAT ( IHO# 58 X, 12HS T AT ION D ATA / 1H0» 1 X, 2HNZ » 5X, 3HX / C» lOX, 5HTHET A, 1 B 246 

IIX. IHSf 13X, 2HCP»11X» 3HCQL» lOX , 5HUEUFS, 9X, 5HWEUFS » 9X » 6H0WEUFS . flX, 5H B 247 

2PEPFS/(1H ,I3»9E14.6)) B 248 

43 FORMAT ( 4 1 1» 1 3, 2X» 2F 1 5 , 1) B 24<> 

END B 250- 



* T I 0 Y * 

SUBROUTINE PROFIL (L) 

SUBROUTINE PROFIL (L) C 1 

COMMON /BLCO/ NZ T» NZ ^ NP» I T» X, ROF S > C MAC H, TT, ET A ( 1 51 ) » DE T A U 51 ) » A ( 15 C 2 

1U»Y{ 151)*IPANPA C 3 

COMMON /BLCl/ HE> P R» CMUFS » UF S > C E L (1 5 1 ) » BETA 1 ( 1 51 ) , UE ( 1 5 1 ) » WE ( 15 1 ) , C A 

1Z(151),PE(151),PHI{151),RH0E( 151),XC(151)»CMUE(151),P1(151)>P3( 151 C 5 

2)»PA(151)»RR(151)»0LP(151),DDW(151) C 6 

COMMON /PROF/ DE L V ( 1 5 1 ) » F (10 1, 2 ) » U( 101, 2) » V ( 1 01» 2 ) » 6 ( 101» 2) , W ( 101» C 7 

12),T(101f 2),B( 101,2 ),C (151), BG( 101,2 ),F( 101,2 ),DENR( 101,2 ),CA1( 101 C 8 

2,2),CA2( 101,2) C 9 

IF (L.EQ.l) GO TO 2 C 10 

C C 11 

E(l,2)»0.0 C 12 

BG(1,2)»1.0 C 13 

BGl«A.0’t'(B6(l,2)-l«0) C lA 

BG2-A .0* (1.0-BG (1,2 ) ) C 15 

DO 1 J»1,NP C 16 

ETAB«ETA(J)/ETA (NP) C 17 

F( J,2)»0.5*ETAB*£TA( J) C 18 

U(J,2)-ETAB C 19 

V( J,2)«1.0/ETA (NP) C 20 

G( J,2)»F( J,2) C 21 

W( J,2)»U(J,2) C 22 

T( J,2)«V(J,2) C 23 

DENR ( J,2 )»1.0 C 2A 

B(J,2)«1.0 C 25 

C(J)»1.0 C 26 

BG(J,2)-1.0 C 27 

1 CONTINUE C 28 

RETURN C 29 

C C 30 

2 NPl-NP+1 C 31 

NPll-NPl-1 C 32 

NP«NP+3 C 33 

IF (NP.6T.101) NP-101 C 3A 

KK-1 C 35 

IF (NZ.EQ.l) KK«2 C 36 

DO A K»KK,2 C 37 

DO 3 J«NP1,NP C 38 

DENR (J,K)-OENR (NP11,K) C 39 

C(J)»1.0 C AO 

FU,K)-ETA( J) + F (NP11,K)-ETA(NP11) C A1 

U(J,K)«1.0 C A2 

V( J,K)»V(NP11,K) C A3 

2 G( J,K) «ETA( J ) + G (NPll, K)-ETA( NPll ) C AA 

W(J,K)«1.0 C A5 



00 

o 


♦ TIDY* 


SUBROUTINE PROFIL (L) 


TU^K)»T(NP11#K) 

C 


B( J,K)«B (NP11»K) 

C 

<»7 

IF (CMACH.EQ.0.0) GO TO 3 

C 

<»8 

BG( K)«1.0 

C 

49 

E( J»K)»E(NP11»K) 

C 

50 

CA1(J,K)»CA1(NP11»K) 

C 

51 

CA2( J,K)»CA2(NP11,K) 

C 

52 

CONTINUE 

C 

53 

CONTINUE 

C 

54 

RETURN 

C 

55 

END 

C 

56 



♦ TIDY* 


SUBROUTINE CUBIC ( YL » XL» IN» FI> NR» PR ) 


SUBROUTINE CUBIC (YL> XL» IN> F I ,NR» PR ) D 1 

DIMENSION YL(1), XL(1)» FI(1), PR(1) D 2 

COMMON /BLCO/ N ZT ,N Z »NP» I T, X, POF S» CMACH» TT, ET A ( 1 51 ) » OE TA ( 1 51 ) » AS( 1 D 3 

151 ), YK151) > IPANPA 0 A 

DO 8 I-lfNR 0 5 

DO 2 J»1»IN D 6 

IF ((FI(I)-XL(J)).LE.O.O) GO TO 1 D 7 

GO TO 2 D 8 

1 K2«J D 9 

GO TO 3 D 10 

2 CONTINUE 0 11 

K2-IN D 12 

3 IF (I.EQ.l) Kl-lOO D 13 

IF <K2.E0.K1) GO TO 7 D lA 

IF (K2.GT. 2. AND.K2.lt. IN) GO TO 5 D 15 

IF (K2.E0.IN) GO TO A D 16 

L-3 0 17 

60 TO 6 0 18 

4 L»IN-1 D 19 

GO TO 6 0 20 

5 L»K2 0 21 

6 CONTINUE 0 22 

A«-(XL(L-l)-XL(L-2) )*(XL(L)-XL(L-2) )*(XL(L+l)-XL(L-2)) 0 23 

B»(XL(L-1)-XL(L-2))*{XL(L)-XL(L-1) ) + (XL{L+l)-XUL-l)) D 2A 

C»-(XL(L)-XL(L-2))*(XL(L)-XL(L-1) )*(XL{L + 1)-XL(D) D 25 

D«(XUU-1)-XL(L-2))*(XL (L+l)-XL(L-l)) + (XL(L+l)-Xl(L)) 0 26 

7 A1-(FI(I)-XL(L))*(FI(I)-XL(L+1)) D 27 

A6»(FI(I)-XL(L-2))*(FI(I)-XL(L-1)) D 28 

C IFUHNA.EQ.A) WRITE(6»15) A # B^ C » D» A1 > A6 D 29 

C15 FQRMATdH » 3 X» 6 ( E 11 . A, 2X ) ) 0 30 

PR(1)-(FI (I )-XL (L-l))*Al*YL(L-2)/A+(Fni)-XL(L-2))*Al*YL(L-l)/B + (F D 31 

lI(I)-XL(L+l))*A6*YL(L)/C+(FI(I)-XL(L))*YL(L + n*A6/D D 32 

K1«K2 0 33 

8 CONTINUE 0 3A 

RETURN 0 35 

END D 36- 


00 



CO 


♦ TIDY* 


SUBROUTINE SPLINE (X»FI»IN,XP) 


SUBROUTINE SPLINE {X»FI»IN»XP) 

DIMENSION X(l), FKl), XP(1), QJ<131)f UJ(131) 

QJ(D— 1.0 

UJ ( l)-2,0*( X(2)-X(l) ) /(FI (2)-FI (1 ) ) 

DO 3 I-2»IN 
AJ«FI(I)-FI(I-1) 

IF (I.EQ.IN) GO TO 1 
BJ«FI(I+1)-FI( I) 

CJ«AJ/( AJ+BJ ) 

DJ»3.0*(CJ*(X( I+1)-X(I ) )/BJ + (l,0-CJ )*(X(n-xn-D) / AJ) 
GO TO 2 

1 DJ-2.0*(X(I )-X( I-l) ) /AJ 
CJ-0.0 

2 PJ»(1.0-CJ)*QJ(I-l)+2.0 
IF (I.EQ.IN) PJ-PJ-1.0 
QJ( n— CJ/PJ 

UJ(I )»(DJ-(1.0-CJ )*UJ(I-1) ) /PJ 

3 CONTINUE 
XP(IN)»UJ(IN) 

INMl-IN-1 

DO I-1#INM1 
NR- IN-I 

XP( NR)-QJ(NR)*XP(NR*1 )+UJ ( NR) 

A CONTINUE 

RETURN 
END 


E 1 
E ? 

E 3 
E 

E 5 
E 6 
E 7 
E 8 
E 9 
E 10 
E 11 
E 12 
E 13 
E lA 
E 15 
E 16 
E 17 
E 18 
E 19 
E 20 
E 21 
E 22 
E 23 
E 2A 
E 25 
E 26- 



♦ TIDY* 

SUBROUTINE INTEG {X»Y#TAB#NPT> 


SUBROUTINE INTEG ( X, Y> T A8> NPT ) F 1 

DIMENSION XU), Y(l), TAB(l) F ? 

IF (NPT.LT.A) GO TO 6 F 3 

DO 5 I»1^NPT F A 

TAB(I)»0.0 F 5 

K«I-1 F 6 

IF (1-2) 5»3,1 F 7 

1 IF II.LT.NPT) GO TO 2 F 8 

K»K-1 F P 

2 K«K-1 F 10 

3 A-X(T) F 11 

B»X(I-1) F 12 

L-K+1 F 13 

M»K+2 F lA 

N«K+3 F 15 

DO A J-1,A F Ife 

XL-X(L) F 17 

XN»X(N) F 18 

XM«X(M) F 19 

XK«X(K) F 20 

YK«Y(K) F 21 

S'JM«YK/( ( XK-XM )♦ (XK-XN)*(XK-XL) ) F 22 

SUM1»((A**A)-(B**A) )/A.O F 23 

SUM2»(XL+XM+XN)*((A**3)-(B**3))/3.0 F 2A 

SUM3»(XM*XN + XM*XL + XL*XN)*(A**2-B**2 )/2.0 F 25 

SUMA«(XM*XN*XL )*(A-B) F 26 

SUM»SUM*(SUM1-SUM2+SUM3-SUMA) F 27 

TAB(I)«TAB(I)+SUM F 28 

ITFMP»K F 29 

K«N F 30 

N«M F 31 

M»L F 32 

A L-ITEMP F 33 

TAB(I)»TA0(I )+TAB(I-l) F 3A 

5 CONTINUE F 35 

6 RETURN F 36 

END F 37- 


00 

Oo 



00 


♦ TIDY* 


SUBROUTINE FLUID 


SUBROUTINE FLUID G 1 

COMMON /BLCO/ NZT, NZ, NP, I T, X, ROF S, C MAC H, TT, E TA( 15 1 ) » DE T A ( 1 51 ) , A ( 15 G 2 

11) > Y( 151)»IPANPA G 3 

COMMON /BLCl/ HE» PR» CMUFS ^ UF S ,C EL ( 1 51 ) » BE TA 1 (151 ) » UE ( 151 ) » W E ( 15 1) , G 

1Z(151),PE(151)» PHI(151)»RH0E(151)*XC(151)»CMUE(151)»P1(151)»P3(151 6 5 

2),P4(151),PR(151),BLP(151),DDW(151) G 6 

COMMON /PROF/ DEL V (1 5 1 ) > F (1 01 » 2 ) ♦ U (101, 2) , V ( 1 01 » 2 ) , G ( 10l» 2 ) » W ( 1 01 » G 7 

12) ,T(101,2),B(101,2),C(151>,BG(101,2),E(101,2),DENR(101,2),CA1(101 G B 

2»2),CA2(101,2) G 9 

COMMON /BOX/ VIS(101),DMUDT(101) 6 10 

COMMON /BAB/ H G 11 

WW-0.0 G 12 

IF (IT.GT.l) GO TO 1 G 13 

PE35»3.5*PF (NZ) G 1<» 

UE2H«0.5*UE (NZ )**2 G 15 

C G 16 

1 DO 2 J-1,NP G 17 

IF (NZ.GT.l) WW«W(J,2) G 18 

H.HE*BG( J,2 )-UE2H*(U(J,2)+*2 + P<r(NZ)*WW**2) 6 19 

IF (H.LT.0.0) RETURN G 20 

TTT-H/6006.0 G 21 

CMU«2.27E-O8*TTT**1.5/(TTT+190.6) G 22 

VIS(J)«CMU G 23 

SP»1.5/TTT-1./(TTT+198.6) G 2A 

OMUDT( J)«CMU*SP G 25 

DENR( J,2)-RH0E (NZ )*H/PE35 G 26 

C( J)»CMU/(CMUE (NZ)*OENR( J,2) ) G 27 

2 CONTINUE G 28 

C G 29 

UE2HE«UE (NZ )**2/HE 6 30 

RPR» 1.0-1.0/PR G 31 

DO 3 J»1,NP G 32 

CAK J,2)-C( J )/PR G 33 

CA2( J,2)»(C ( J)*UE2HE)*RPR*(U( J,2)*V(J,2)+P4(NZ)*W(J,2)*T( J,2) ) 6 3A 

B(J,2)-C(J) G 35 

3 CONTINUE G 36 

RETURN G 37 

END G 38- 



♦ TIDY* 


SUBROUTINE COEF 


SUBROUTINE COEF H 1 

COMMON /BLCO/ NZT» NZ» NP, I T, X, ROF S» C MAC H, TT, ETA( 151 ) » DE T A ( 1 51 ) * A ( 15 H 2 

11) ,Y{151)»IPANPA H 3 

COMMON /BLCl/ HE» PR > C MUFS » UF S , C E L ( 1 51 ) » B ETA 1 ( 151 ) » UE ( 151 ) , W E ( 151 ) » H 4 

1Z(151)»PE(151), PHI(151),RHOE( 151)>XC(151),CMUE(15l),Pl{151)»P3( 151 H 5 

2)»P4(151),RR{15l)^BLP{151),0DW(151) H 6 

COMMON /PROF/ DEL V ( 1 5 1 ) » F ( 1 01 > 2 ) , U ( 101» 2 ) » V (101» 2 ) » G ( 101 » 2 ) > W ( 101» H 7 

12) »T(101»2)>B(101»2)»C(151)» BG( 101 » 2 ) » E ( 1 01» 2 ) » DENR ( 101» 2 ) ^ C A 1 ( 101 H B 

2,2),CA2(101,2) H 9 

COMMON /BLC8/ Bl ( 101 ) » B2 ( 101 ) ,B 3 ( 101) » B4 (101 ) » B 5 ( 101 ), B6 ( 101 ) , B7( 1 H 10 

101),B8( 101) >B9( 101) fB10( 101), Rl( 101 )»R2( 101 )»R 3(101) »R 4(101 ),R5( 10 H 11 

21),R6(101), SKlOl), S2(101 ),S3 (101),S4(101),S5(101),S6(101),S7(101) H 12 

3, SB(lOl), 59(101), SlO(lOl) H 13 

P1P»P1(NZ)+CEL(NZ) H 14 

P3P«P3(NZ)+CEL(NZ) H 15 

P4P»P4(NZ)-CEL (NZ) H 16 

P1T2-2.0*P1(NZ) H 17 

P4T2«2.0*P4(NZ) H 18 

P1P2-2.0+P1P H 19 

DO 4 J«2,NP H 20 

UB-0.5*(U(J,2)+U(J-1,2) ) H 21 

V3»0,5*( V(J,2)+V(J-1,2) ) H 22 

GB»0.5*(G(J,2)+G(J-1,2) ) H 23 

We«0.5*(W(J,2)+W(J-l,2) ) H 24 

TB»0.5*(T(J,2)+T(J-1,2) ) H 25 

DENRB-0.5*(DENR(J,2)+DENR(J-1,2)) H 26 

FVB«0.5*(F(J,2 )*V(J,2)+F(J-1,2)*V(J-1,2)) H 27 

FTB»0,5*(F( J,2)*T(J,2)+F( J-1,2)*T(J-1, 2)) H 28 

UWB«0.5*(U(J,2)*W(J,2)+U(J-1,2)*W(J-1,2)) H 29 

GVB«0.5*(G(J,2)*V(J,2)+G( J-1,2)*V(J-1,2)) H 30 

GTB-0.5*(G( J,2)*T(J,2)+G(J-1,2)*T(J-1,2)) H 31 

WSB-0.5*(W( J,2)**2+W(J-1,2)**2) H 32 

IF (NZ.GT.l ) GO TO 1 H 33 

C H 34 

CUB-0.0 H 35 

CVB-0.0 H 36 

CGB-0.0 H 37 

CWB-0,0 H 38 

CTB-0.0 H 39 

CFTB-0.0 H 40 

CUWB-0.0 H 41 

CGVB-0.0 H 42 

CGTB-0.0 H 43 

CWSB-0,0 H 44 

CDENRB-0.0 H 45 


00 

Ul 
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SUBROUTINE COEF 




H 

^6 

- ATTACHMENT-LINF FLOW 


H 

A7 

DEFINITIONS OF COEFFICIENTS IN DIFFERENCED X-MOM 

EG. 

H 

<t8 

Sl( J)«B(J»2) +A(J )♦(-!. 5*F(J, 2)+Pl{NZ)*G(J,2)-BLP(NZ)) 

H 

A9 

S2( J)— B(J-1»2)+A(J)*(-1.5*F(J-1*2)+P1(NZ)*6( J“l» 

2)-BLP(NZ) ) 

H 

50 

S3( J)— 1.5*A(J )*V(J»2) 


H 

51 

S4( J)»-1.5*A»J)*V(J-1>2) 


H 

52 

S5( J) «0.0 


H 

53 

S6( J)«0.0 


H 

5A 

S7( J)»A( J)*P1(NZ)*V( J,2) 


H 

55 

S8( J)«A( J)*P1(NZ)*V(J-1»2) 


H 

56 

S9(J)«0.0 


H 

57 

S10( J)-0.0 


H 

58 

DEFINITIONS OF COEFFICIENTS IN DIFFERENCED Z-MOM 

EG. 

H 

59 

B1(J)-S1( J) 


H 

60 

B2(J )«S2(J) 


H 

61 

B3( J)»-1.5*A(J)*T(J^2) 


H 

62 

B<t( J)»-l .5*A(J)*T(J-1»2) 


H 

63 

B5( J)— A(J)*(-U(J>2) + P1T2*W( J»2)) 


H 

66 

B*>( J) — A( J)*(-U(J-1»2) + P1T2*W( J-lf 2 )) 


H 

65 

B7( J)-A( J)*W( J>2) 


H 

66 

B9( J)«A(Jj)*W(J-l,2) 


H 

67 

B9( J)«A( J)*P1(NZ)*T(J,2) 


H 

68 

B10( J)-A(J)*P1 (NZ)*T(J-1, 2) 


H 

69 

SO TO 2 


H 

70 



H 

71 

GENERAL FLOW 


H 

72 

CJ3»0.5*(U( J»1) + U(J-1^ 1 ) ) 


H 

73 

CVB-0.5*(V( J»1)+V(J-1,1) ) 


H 

76 

CSB«0.5*(G( J»1)+G(J-1»1) ) 


H 

75 

CWB»0.5*(W( J,1)+W(J-1,1)) 


H 

76 

CTB»0.5*(T( J»1)+T{J~1» 1) ) 


H 

77 

CFTB»0.5*(F (J,1)*T(J,1)+F(J-1,1)*T(J-1,1) ) 


H 

78 

CFVB«0,5*(F ( J,1 )*V( J*1) + F{J-1>1) + V(J-1» in 


H 

79 

CUWB«0 ,5*(U( J» 1 )*W(J, 1 )+U( J-1»1)*W( J-l»l ) ) 


H 

30 

CSVB-0.5* (G(J,1 )*V( J,l)+G(J-l,n*V( J-l» 1) ) 


H 

81 

CGTB«0,5*(G( J»1 )*T(J»1)+G{J-1,1)*T( J-l,in 


H 

82 

CWSB»0.5*{W(4,1)**2+W(J-1»1)**2) 


H 

83 

CDENRB-0.5+ (DENR( J»1)+DENR(J-1* 1 ) ) 


H 

86 



H 

85 

DEFINITIONS OF COEFFICIENTS IN DIFFERENCED X-MOM 

EG* 

H 

86 

SI ( J)«B( J»2 ) + A ( J )'M-1.5*F( J,2) + P3P*G(J»2)-CEL(NZ)*CGB-BLP(NZ) ) 

H 

87 

S2(3 )— B (J-l,2)+A( J )*(-!. 5*F(J-1»2 )+P3P*G(J-l, 2)- 

-CEL(NZ)*CGB-BLP(N 

H 

88 

IZ) ) 


H 

89 

S3( J)— 1.5*A( J)*V( J,2) 


H 

90 
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00 


S'KJ)— 1.5*A(J)*V( J-1,2) H 91 

S5( J)«A( J)*(PAP*W(J>2)-CEL(NZ)*CWB) H 92 

S6( J)-A(J)*(P^P*W(J-1»2)-CEL(NZ)*CWB) H 93 

S7(3)«A(J)*(P3P*V(J,2)+CEL(NZ)*CVB) H 9A 

S8( J)»A( J)*(P3P*V(J-1»2)+CEL(NZ)*CVB) H 95 

S9( J)«A( J)*(P4P*U( J>2)-PAT2*W(J#2)+CEL(NZ)*CUB) H 96 

S10(J)"A(J)*(PAP*U(J-1,2)-P<»T2*W{J-1»2)+CEL(NZ )*CUB) H 97 

C H 98 

C DEFINITIONS OF COEFFICIENTS IN DIFFERENCED Z-NO« EQ, H 99 

B1(J)»S1(J) H 100 

B2(J)»S2(J) H 101 

B3( J)«-1.5*A(J)*T(J,2) H 102 

B'H J)— 1.5*A(J)*T(J-1»2) H 103 

B5(J)— A(J)*(P1P2*W{J»2)-U(J>2)) H 104 

B6( J)«-A( J)*(P1P2*W(J-1»2)-U( J-l»2) ) H 105 

B7(J)«A(J)*W(J#2) H 106 

B8('J)-A(J)^W(J-1»2) H 107 

B9( J)«A( J)*(P3P*T(J,2)+CEL(NZ)*CTB) H 108 

B10( J)«A( J)*(P3P*T(J-1»2)+CEL (NZ)*CTB) H 109 

H 110 

DEFINITION OF RJ H 111 

R1U)»F(J-1»2)-F(J»2)+DETA( J-D+UB H 112 

R2( J)-U(J-1>2)-U(J»2)+DETA(J-1)*VB H 113 

R3( J)«G(J-1#2)-6(J»2)+0ETA( J-1)*WB H 114 

R4( J)-W( J-1»2)-W(J»2)+DETA( J-D + TB H 115 

IF (NZ.GT.l ) GO TO 3 H 116 

R5( J) — (B(Jf 2>*V( J»2)-B( J-1, 2)*V( J-1»2)+DETA( J-l)*(-l. 5*F V8 +P 1 ( NZ ) H 117 
1*3VB-BLP(NZ)*VB)) H 118 

R6( J )»-DETA(J-l)*DENRB*(-1.0 + Pl(NZ) )-(B(J»2)*T(J»2)-B(J-l»2)*T( J-1 H 119 
l/2)+0ETA(J-l)*(-1.5*FTB+Pl(NZ )*GTB+ UWB-Pl (NZ ) +WSB-BL P ( NZ ) ♦TB ) ) H 120 

60 TO 4 H 121 

3 DERBV»(B(J, 1)*V(J»1)-B( J-l»l )*V(J-1,1) )/DETA(J-l) H 122 

CL 5B-DERBV>1.5*CFVB+P3«NZ-1)*CGVB+P4(NZ-1)'MCUWB-CWSB)-BLP(NZ-1)*C H 123 

IVB H 124 

CR5B«-CL5B+CEL(NZ)*(C6VB-CUWB) H 125 

R5( J) -DETA( J-1 )*CR5B-(B(J>2)*V(J»2)-B( J-1,2)*V( J-1#2)+DETA(J-1)*(- H 126 
11.5*FVB+P3P*GVB+P4P4UWB-P4(NZ)*WSB-CEL(N7 )♦ (CWB* UB-CUBTWB-CVB*GB+C H 127 
2GB*VB)-BLP(NZ)*VB)) H 128 

DERBT«(B{J»1)*T(J,1)-B( J-1,1)*TU-1,1)»/0ETA(J-1) H 129 

CL6B«DER8T-1.5*CFTB+P3(NZ-1 )*CGTB+P1(NZ-1 )'MCDFNRB-CWSB)+CUWB-CDEN H 130 
1RB-BLP(NZ-1 )+CTB H 131 

CR6B«-DENRB*(P1(NZ)-1.0)+CEL( NZ )*(CGTB-CWSB)-CL6B H 132 

R6( J )-DETA( J-l)*CR6B-( B( J, 2) *T( J, 2 ) -B ( J-1 »2 ) *T { J-1, 2 ) +DETA ( J-1 ) * (- H 133 
11,5*FTB + P3P*GTB-P1P*WSB+UWB-CEL (NZ)'f(CGB*TB-CTB*6B)-BLP(NZ)*TB) ) H 134 

4 CONTINUE H 135 
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H 136 
H 137- 


RETUPN 

END 
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SUBROUTINE SQLV6 


SUBROUTINE S0LV6 I 1 

COMMON /BLCO/ NZ T, NZ » NP, I T, X , ROF S, C M AC H, TT, ET A ( 1 51 ) » DE TA { 1 51 ) » A ( 15 I 2 

11) > Y(151 )>IPANPA I 3 

COMMON /BLCl/ HE ^ PR^ C MUFS , UF S ^C EL { 1 51 ) , BE TA H151 ) , UE ( 1 5 1 ) , WE ( 15 1) , I A 

1Z(151)»PE(151)>PHI(151)>RH0E( 151)»XC(151)>CMUE(151)»P1{151),P3(151 I 5 
2)»P<f (151),RR(151),3LP(151),DDW( 151) I 6 

COMMON /PROF/ DEL V ( 1 5 1 ) > F ( 101 , 2 ) » U ( 101, 2 ) , V ( 101, 2 ) , G ( 101 » 2 ) » W ( 1 01 , I 7 

12) ,T(101,2),B(101»2),C(151),BG(101,2),E{101,2),DENR(101,2),CA1(101 I 8 

2,2),CA2(101,2) I 9 

COMMON /BLC8/ B 1 ( 101 ) , B2 ( 101 ) ,B3 { 10 1 ) , B A { 101 ) , B 5 ( 1 01 ), B6 ( 101 ) , B7( 1 I 10 
101) ,B8( 101) ,B9(101),B10( 101) , R1 (101 ),R2 (101 ),R3( 101 ),RA(101), R5(10 I 11 
21),R6(10l),Sl(101),S2(101),S3(101 ),SA(10l),S5(101),S6(101),S7(101) I 12 
3,S8(101),S9(101 ),S10( 101) I 13 


DIMENSION AlKlOl), A2K101), A3K101), 441(101), A5K101), A6K10 I 14 
11), A12(101), 422(101), 432(101), 442(101), 452(101), 462(101), B1 I 15 
21(101), 321(101), B3K101), 841(101), B5K101), B6K101), B12(101) I 16 
3, 622(101), 632(101), 642(101), 852(101), 862(101), DELF(lOl), DEL I 17 
4U(101), DELT(lOl), DELG(lOl), DELW(lOl), WKlOl), W2(101), W3(101) I 18 


5, W4(101), W5(101), W6(101) I 19 

C CALCULATION OF GAMMA (AI1,AI2 1-1,6) VECTOR FOR J-2 I 20 

C FIRST All I 21 

A11(2)-(S5(2)+S1(2)/A(2)+S3(2)*A(2) )/(S2(2)-Sl(2)) I 22 

A21 (2)-( B7( 2)+B3(2)*A(2) )/(B2(2)-Bl (2) ) I 23 

431(2) — A(2) I 24 

A4K2)— All (2)-l./A(2) I 25 

A51(2)«0.0 I 26 

A6K2)— A21 (2) I 27 

C THEN AI2 I 28 

A12(2)»(S7(2)*A(2)+S9(2))/(S2(2)-S1(2)) I 29 

A22(2)-(B5(2)+A(2)*B9(2)+B1( 2)/A(2) )/(B2(2)-Bl(2)) I 30 

A32(2)«0.0 I 31 

A42(2)«- A12(2) I 32 

A52(2)«-A(2) I 33 

A62(2)— A22 (2 )-!♦/ A(2) I 34 

C I 35 

C CALCULATION OF WKI-1,6) A1 ( V EC TOR ) ( VECTOR ) -R ( VECTOR ) , AT J»2 I 36 

Wl(2)-(R5(2)+(R2(2)tSl(2))/A(2)-S7(2)+R3(2)-S3(2)*Rl(2))/(S2(2)-Sl I 37 

1(2)) I 30 

W3(2)»R1(2) I 39 

W4( 2)— Wl(2)-R2(2)/A(2) I 40 

W5(2)»R3(2) I 41 

W2(2)-(P6(2 )-Bl (2)*R4(2)/A(2)-B9(2)*R3(2)-B3(2)*R1(2) )/(B2(2)-Bl(2 I 42 

D) I 43 

W6(2)— W2(2)-R4(2)/A(2) I 44 

C I 45 



VO 

O 
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C CALCULATION OF ALFA COEFFICIENTS BI1»BI2 WITH I-l»6 I ^6 

C NOTE-THE SUBSCRIPT FOR THESE COEF. START FROM 11. I A7 

DO 1 J«3»NP I 48 

Bll ( J)»-A(J )+A31(J-l) I 49 

B21( J )»-1.0+A( J )*A41( J-1) I 50 

B31 ( J)«A51( J-1 ) I 51 

B41( J)-A( J) *A61 (J-1 ) I 52 

B5K J)«S6(J )-S4 (J )*A31( J-l)-S2( J )*A4UJ-1)-S8(J )*A51(J-1) I 53 

B61( J)«B8(J)-B4(J)*A31(J-1)-B10(J)*A51(J-1)-B2(J)*A61(J-1) I 54 

B12( J)-A32( J-1 ) I 55 

B22( J)»A42( J-1)*A(J) I 56 

B32( J)-A52( J-1 )-A( J ) 1 57 

B42{J)»A62(J-1 )*A(J)-1.0 I 58 

B521J )»-(S4(J)*A32( J-l)+S2( J )*A42( J-1)+S8(J )* A52 { J -1 ) ) +S10 ( J ) I 59 

B62 ( J)-B6(J)-B4(J)+A32( J-1)-B10(J)*A52(J-1)-B2( J)*A62( J-1) I 60 

C I 61 

C CALCULATION OF AI1»AI2 WITH 1-1,6 I 62 

CCA1»B51(J)-S3(J)*B11(J) + S1(J)*B2UJ)/A(J)-S7(J)*B31(J) I 63 

CB1«B52( J)+S1(J )*B22(J)/A( J)-S7( J)*B32(J)-S3(J)*B12( J) I 64 

C:A2-B61 ( J) -B3 (J )*B11(J )-B9( J)*B31( J) +B1 ( J) *B41( J) /A(J ) I 65 

CB2>iB62U)-B3(J)*B12(J)-B9(J)*B32(J)+Bl(J)*B42(J)/A(J) I 66 

CC1«S5(J)+S3(J)*A(J)+S1{J)/A(J) I 67 

CC2»B7(J)+B3(J)*A(J) I 68 

DEN-CCAl+CB 2-CB1+CCA2 I 69 

A11(J)-(CC1*CB2-CB1*CC2)/DEN I 70 

A2K J)»(CCA1*CC2-CC1*CCA2)/DEN I 71 

A3K J)»-A(J)-B11(J)*A11(J)-B12(J)*A2KJ) I 72 

A41(J)«(-1.0+B21(J)*A11(J)+B22(J)*A21(J))/A(J) I 73 

A51 (J)— B31(J)*A11( J)-B32(J)*A21(J ) I 74 

A61(J)-(B41(J)*A11(J)+B42(J)*A21(J))/A(J) I 75 

CC1-S7(J)*A(J)+S9(J) I 76 

CC2»B5(J)+B9(J )*A(J)+B1( J)/A( J) I 77 

A12 ( J )• (CC1+CB2-CB1+CC2) /DEN I 78 

A22( J)»(CCA1*CC2-CC1*CCA2)/DEN I 79 

A32U)— Bll (J)*A12( J)-B12(J)*A22(J ) I 80 

A42(J)-(B21(J)*A12(J)+B22(J)*A22(J ) )/A(J) I 81 

A52( J)«-A(J)-B31(J)*A12( J)-B32( J)*A22(J) I 82 

A62( J)»(-1.0+B41(J)*A12(J)+B42(J)*A22(J))/A(J) I 83 

01-R1(J )+W3(J-l) I 84 

D2-R2( J)+A(J )*W4(J-1) I 85 

03-R3( J)+W5 (J-1 ) I 86 

D4»R4( J )+A( J )7W6( J-1 ) I 87 

D5»R5(J)-{S4(J)*W3(J-1)+S2(J)*W4(J-1)+S8(J)*W5(J-D) I 88 

D6-R6( J)-(B4 ( J )*W3(J-1)+B10( J)*W5( J-1 )+B2 (J )*W6( J-1 ) ) I 89 

CC1-D5-D1*S3(J )+(Sl ( J )*D2)/A ( J)-S7( J)*03 I 90 
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CC2«D6-D1*B3( J) + (B1( J)*0^»)/A ( J)-B9( J)*D3 
Wl( J)»(CC1*C82-CB1*CC2)/DEN 
W2( J)»(CCA1*CC2-CC1*CCA2)/DEN 
W3( J)-01-B11(J )*W1(J)~B12(J)*W2(J) 

W<*(J )»{-D2+B21( J)*W1( J)+B22( J)*W2(J ))/A(J ) 

W5( J)-D3-B31(J)*W1(J)-B32(J)*W2(J) 
W6U)»(-D4+B41( J)*W1( J) + B42( J )*W2(J ))/AU) 
CONTINUE 

CALCULATION OF PERTURBATION QUANTITIES 
DELU(NP)«0.0 
DELW(NP)»0.0 
0ELF(NP)«W3(NP ) 

OELVlNEj-W^tlNP) 

DELG(NP)«W5(NP) 

DELT(NP)«W6(NP ) 

DELU(NP-l)-WKNP) 

OELWINP-l )«M2(NP) 

J»NP 

2 J-J-1 

DELU( J-1)-WI{J)-A11( J)*DELU( J )-A12(J)*DELW( J) 
DELW(J-1)«M2(J )-A21( J)*DELU(J )-A22(J)*DELW(J) 
DELF( J)»W3( J)-A31(J)*0ELU( J)-A32(J)*DELW(J) 
0ELVU)«W^( J )-A41 { J )*0ELU( J )-A<t2(J )*DELW{ J) 
0ELG( J)»W5( J)-A51(J )*DELU(J )-A52(J)*DELWU) 
DELT( J)»W6( J)-A61(J)*0ELU(J)-A62U)*DELW{ J) 

IF (J.GT.3) GO TO 2 

DELF(2)»W3(2)-A31(2)*DELU(2)-A32(2)+DELW( 2) 

DELV(2)»WA(2)-AA1(2)*DELU(2)-A42(2)*DELW(2) 

DELG(2)«W5(2)-A52(2)*DELW(2)-A51(2) *DELU{2) 

0ELT{2)»W6(2)-A61(2)*DELU(2)”A62(2)*0ELW(2) 

DELV(1)-W1(2)-A11(2)*DELU(2)-A12(2)*DELW(2) 

DELT(l)«W2(2)-A21(2)*DELU(2)-A22(2)*OELW(2) 

DELF(1)«0.0 

0ELG(l)-0.0 

0ELU(1)»0.0 

DELW(l)-0.0 

C IF(IT.EO.l) WRITE(6»A) 

C WRITE(6,5) IT,V(1»2),DELV(1),T(1,2)»DELT(1) 

DO 3 J»1#NP 
F{ J,2)»F(J,2)+DELF(J) 

U( J,2)»U(J,2)+0ELU(J) 

V( J>2)»V( Ji2)+DELV(J) 

G( J*2)»G(J»2)+DELG( J) 

W( J*2)«W(J«2)*DELW(J ) 


I 91 
I 92 
I 93 
I 94 
I 95 
I 96 
I 97 
I 98 
I 99 
I 100 
I 101 
I 102 
I 103 
I 104 
I 105 
I 106 
I 107 
I 108 
I 109 
I 110 
I 111 
I 112 
I 113 
I 114 
I 115 
I 116 
I 117 
I 118 
I 119 
I 120 
I 121 
I 122 
I 123 
I 124 
I 125 
I 126 
I 127 
I 128 
I 129 
I 130 
I 131 
I 132 
I 133 
I 134 
I 135 



VO 

K> 
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T( J»2)»T(Ji2)+DELT{ J) 
3 CONTINUE 
RETURN 
C 
C 

END 


I 136 
I 137 
I 138 
I 139 
I 140 
I 141- 
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SUBROUTINE SOLV2 


SUBROUTINE S0LV2 J 1 

COMMON /BL:0^ NZT,NZfNP,IT,X,R0FS»CMACH,TT,ETA(151)»DETA(151),A(15 J 2 

11) »Y(151),IPANPA J 3 

COMMON /BLCl/ HE» PR» CMUFS^ UFS# C EL ( 1 51 ) » BE TA I ( 151 ) , UE ( 151 ) , WE ( 15 1 ) , J ^ 

lZ(151)f PE(151)> PHI(151)fRHQE(151)»XC(151),CMUE(l51)»Pl(151)>P3(15i J 5 

2)»P<»(151),RR(151),BLP(151),DDW(151) J 6 

COMMON /PROF/ DEL V ( 151 ) » F ( 101# 2 ) » U ( 101» 2 ) » V ( 101» 2 ) » G ( 1 01# 2 ) > W ( 10l> J 7 

12) ,T(101»2)#B(101,2)»C{151)#BG(101,2)#E{101»2)#DENR{101#2)#CA1{101 J 8 

2#2)#CA2(101#2) J 9 

DIMENSION SiaOD# $2(101)# S3(10D# RKIOD# R2(10D# YKlOl), Y2 J 10 

1(101), BIKIOD# B12(10D# AIKIOD# A12(101) J 11 

BG(NP,2)»1.0 J 12 

ALFAO-O.O J 13 

ALFAl-1.0 J 1^ 

GAMMAO-0.0 J 15 

BETA0»1.0 J 16 

BTAl-0.0 J 17 

C J 18 

DO 2 J-2#NP J 19 

FB«0.5*(F(J#2)+F(J-1#2)) J 20 

GB»0.5*(G(J#2)+G(J-1#2)) J 21 

WB»0. 5*(W(J,2)+W(J-1,2) ) J 22 

IF (NZ.GT.l) GO TO 1 J 23 

C J 2^# 

CFB-0.0 J 25 

CGB«0.0 J 26 

CWB-0.0 J 27 

CEB-0.0 J 28 

C - ATTACHMENT-LINE FLOW J 29 

S1(J)-CA1(J#2)+A(J)*(-1.5+FB+P1(NZ)+GB-BLP(NZ)) J 30 

S2(J)--CA1( J-1#2)-CA1(J#2)+S1(J) J 31 

S3(J)-0.0 J 32 

Rl( J)»CA2(J-1#2)-CA2(J#2) J 33 

R2(J)»0.0 J 3A 

GO TO 2 J 35 

C J 36 

1 CFB»0.5*(F( J,l)+F( J-l#l) ) J 37 

CGB-0.5* (G( J#1)+G(J-1# 1) ) J 38 

CWB-0,5*(W(J#1)+W(J-1#D) J 39 

CEB»0.5*(E( J#1 )+E (J-1# 1) ) J ^0 

CBGB»0.5*(BG(J#1)+BG(J-1#1)) J ^1 

C J ^2 

SI ( J )«CA1 ( J#2)+A( J)*(-l. 5*FB+P3 (NZ )*GB+CEL(NZ )*(GB-CGB )-BLP(NZ ) ) J 43 

S2(J)--CA1( J-1#2)-CA1(J#2)+S1(J) J 44 

S3(J)»- A (J)*CEL (NZ)*(WB+CWB) J 45 



VO 


♦ T I D Y * 


SUBROUTINE SOLV? 


C J <»6 

0ERCA1«(CA1(J,1)*E(J*1)-CAI( J-1,1)*E(J-1»1))/DPTA(J-1) J 

0ERCA2- { (CA2(J#2)-CA2( J-li2> ) /DE TA ( J-1 ) ) ♦ ( ( C A2( 4» I ) -C A2 ( J-l» 1 ) ) / DE J 
ITA(J-D) J 

CLBE»0ERCAl-1.5'fCFB’*'CEB+P3<NZ-l)*CGB*CEB-BlP(NZ-l)*CEB J 50 

Rl( J) ■0ETA( J-1 )*(-CLBE+CEL(NZ )*(-CBGB*(WB+CWB)-(GB-CGB )*CEB)-DERCA J 51 

12) J 52 

R2(J)-0.0 J 53 

2 CONTINUE J 54 

R2(NP)«1,0 J 55 

C J 56 

Rl(l)»GAMMAO J 57 

R2(l)»0,0 J 58 

B11(1)»ALFA0 J 59 

812(1)-ALFA1 J 60 

Y1(1)»R1(1) J 61 

Y2(1)»R2{1) J 62 

DO 3 J«2»NP J 63 

All { J )« ( $2( J )-A (J )*S3( J ) ) /(B12( J-1)-A(J )*B11C J-1 ) ) J 64 

A12 ( J)»B11( J-1)*A11 ( J )-S3 ( J) J 65 

C CALCULATION OF ALFA COEFFICIENTS J 66 

B11(J)«S3(J)-A12(J) J 67 

B12(J)»S1(J)+A12(J)*A(J) J 68 

Y1(J)-R1(J)-A11(J)*Y1(J-1)-A12{J)*Y2(J-1) J 69 

Y2(J)»R2(J) J 70 

3 CONTINUE J 71 

8G(NP#2)«R2(NP) J 72 

E( NP»2)« {Y1(NP)*BETA0-B11(NP) *Y2(NP) )/(B12(NP)*BETA0-Bll(NP)*BTAl) J 73 

J-NP J 74 

4 J«J-1 J 75 

PAR1«Y2(J)-BG(J+1»2)+A(J+1)*E(J+1»2) J 76 

E( J>2)»(Y1( J)+Bll(J)*PARl)/(-A(J+l)*Bll{J)+B12(J)) J 77 

BG( J#2)»-A( J+1)*E( J,2)-I»AR1 J 78 

IF (J.GT.l) GO TO 4 J 79 

RETURN J 80 

END J 81- 



♦ TIDY* 

SUBROUTINE OUTPUT 

SUBROUTINE OUTPUT K 1 

DIMENSION UP(101)» UPP(101)» WP(lOl), WPP(lOl)# TP(lOl), TPP(101)» K 2 
1 CMUPUOl), CMUPP(lOl)^ ALFAP(lOl) K 3 

COMMON /BLCO/ NZT, NZ, NP» I T, X, ROES, C MAC H, TT, ET A ( 151 ) » OE T A { 151) » A ( 15 K A 

11) >Y{151),IPANPA K 5 

COMMON /BLCl/ HE» PR# C MUF S^ UF S »CE L ( 1 51 ) » BETAl ( 1 51 ) , UE ( 1 51 ) > WE ( 151 ) » K 6 

1Z(151)»PE(151)»PHI(151),RH0E (151)»XC(151),CMUE(151),P1(151),P3(151 K 7 
2)»PA{151)»RR(151)»BLP(151)»D0W(151) K 8 

COMMON /PROF/ DEL V ( 151 ) . F ( 10 1» 2 ) > U ( 101» 2 ) » V ( 101» 2) » G( 1 01» 2) » W ( 101» K 9 

12) »T(101*2),B{101,2),C(151)^BG(101,2)fE(101,2)»DENR(101,2),CAl(101 K 10 

2»2) »CA2( 101> 2) K 11 

DIMENSION YKIOD* Ul(101)» Wl(101)» TKIOD* CMUKlOl), ALFA(lOl) K 12 
DIMENSION AKIOD# A2(101)» A3(101)» AA(10l)» A5(101)» A6(101)» A7 K 13 
1(101)» ABdOl), A9(101)» AlO(lOl), AlKlOl), A12{101), A13(101)» A K lA 
21A(10D* A15(10D# A16(10l), RAA(lOl) K 15 

COMMON /PAR/ A111,A21# A31»VGP»NNN, IPRINT K 16 

COMMON /BOX/ VIS(101)»DMU0T(101) K 17 

COMMON /WAW/ RjXO»CO*ZST»CHLENTH>THETA,WEEA»SO»NEDGE»TEEA»CMUEEA K IB 

C K 19 

RX«RHOE(NZ)*UE (NZ)*X/CMUE(NZ) K 20 

SQRX-SQRT(RX) K 21 

PAR3«X/SQRX K 22 

SUM-OtO K 23 

F1«DENR(1,2) K 2A 

Y(1)«0.0 K 25 

DO 1 J«2»NP K 26 

F2«DENR(J,2) K 27 

SUM»SUM+<F1+F2 )*A(J) K 28 

F1«F2 K 29 

Y( J)«SUM*PAR3 K 30 

1 CONTINUE K 31 

DUMl-0. K 32 

DUM2«0. K 33 

DJM3-0. K 3A 

RES-0. K 35 

DO 2 J«2»NP K 36 

DUM1»1.-W(J-1»2) K 37 

DJM2«1.-W(J»2) K 38 

RES«RES+(DUM1+DUM2)/2.*(Y(J )-Y( J-1) ) K 39 

2 CONTINUE K AO 

DSTZINC-RES K A1 

RDSTZ-RHOE(NZ)*WE(NZ)*DSTZINC/CMUE (NZ) K A2 

CID-SUM K A3 

w DELSTX«PAR3*(-F (NPi2)+CID) K AA 

DELSTZ«PAR3*{-G (NP,2 )+CID) K A5 



v£> 

ON 


♦ TIDY* 


SUBROUTINE OUTPUT 


SUM»0.0 

SUM2-0.0 

F1«U(1,2)*U(1»2) 

F11«W(1»2)*W(1»2) 

DO 3 J-2»NP 
F2»U( J»2)*U(J»2) 

F22«W(J»2)*W{J»2) 

SUM-SUM+(F1+F2)*A(J) 

SUM2-SUM2+(F11+F22)*A( J) 

F1«F2 

F11-F22 

3 CONTINUE 

THETAX»PAR3*(F (NP^Z)-SUM) 

T^ETAZ»PAR3♦(G( NP,2)-SUN2) 

CFX«2.0*C(1 )*V(1»2)/SQRX 

HX-DELSTX/THETAX 

HZ-DELSTZ/THETAZ 

IF (CMACH.EQ.O.O) GO TO A 

TE-PE (NZ)/RH0E(NZ)/1716.0 

TW»TE+DENP(1,2) 

RHOW«RHOE(NZ)/OENR(l,2) 

GO TO 5 
A TE-TT 

TW-TT 
RHOW-ROFS 

5 VW-BLP(NZ )*SQRT(UE (NZ )*CNUE (NZ)*RHOE(NZ)/X) /RHOW 
IF (NZ.GT.l) GO TO 6 

CFZ»0.0 

SQUIG-BLP(l)/SORT(Pl(l>) 

GO TO 7 

6 CFZ«2.0*C(1)*T(1»2)*UE(NZ )/WE (NZ )/SQRX 
S3UIG-BLP(NZ)*S0RT(RR(NZ)*UE(NZ)/WE(NZ)/X) 

7 IF (NZ.EQ.l) GO TO lA 

R-SQRT(RHOE (NZ )*WE(NZ )*RR (NZ ) /CMUE (NZ) ) 

CHL ENTH-SQRT(CMUE (NZ) *RR (NZ ) / (RHQE (NZ )*WE (NZ) ) ) 
DO 8 J«1»NP 
Yl( J)»Y( J)/CHLENTH 
C Yl( J)-Y(J)/DSTZINC 

U1 ( J )«W( J>2) 

W1(J)»U( J»2)*UE(NZ)/WE(NZ) 

T1 ( J)«DENR( J,2 ) 

CMUK J)-VIS( J) /CMUE(NZ) 

ALFA( J)»DMUOT( J )*TE/CMUE (NZ ) 

8 CONTINUE 
C 


K 46 
K 47 
K 48 
K 49 
K 50 
K 51 
K 52 
K 53 
K 54 
K 55 
K 56 
K 57 
K 58 
K 59 
K 60 
K 61 
K 62 
K 63 
K 64 
K 65 
K 66 
K 67 
K 68 
K 69 
K 70 
K 71 
K 72 
K 73 
K 74 
K 75 
K 76 
K 77 
K 78 
K 79 
K 80 
K 81 
K 82 
K 83 
K 84 
K 85 
K 86 
K 87 
K 88 
K 89 
K 90 



♦ TIDY* 

SUBROUTINE OUTPUT 

DO 9 I»1,NP K 91 

A16(I)»PR K 92 

9 CONTINUE K 93 

IF ( IPANPA. EQ.2) 60 TO 12 K 94 

NPM1«NP-1 K 95 

DO 10 J«2»NPM1 K 96 

DY1«Y1(J)-Y1(J-1) K 97 

DY2»Y1( J+1)-Y1(J ) K 98 

UP( J)-(0Y1*U1(J+1 )/0Y2-DY2*Ul (J-D/OYl) / (DY1+0Y2)-U1( J )* (0Y1-DY2)/ K 99 

1(DY1*DY2) K 100 

WP ( J )«(DY1*W1( J + 1)/DY2-0Y2*W1(J-1)/DY1)/(0Y1 + DY2)-W1( J)*(DY1-DY2)/ K 101 
1(0Y1*DY2) K 102 

TPU)»(DY1*T1(J + 1)/DY2-0Y2*T1(J-1)/0Y1)/(0Y1 + 0Y2)-T1( J )*(0Y1-DY2)/ K 103 
1(0Y1*DY2) K 104 

CMUP(J)» (0Y1*CMU1(J+1)/0Y2-0Y2*CMU1(J-1»/0YI)/(0Y1+0Y2 »-CMUl{ J)*(0 K 105 
1Y1-DY2)/ (DY1*DY2) K 106 

ALFAP( J)«(DY1*ALFA(J+1)/OY2-OY2*ALFA(J-1)/OY1)/(OY1+DY2)-ALFA(J)*( K 107 
1DY1-DY2) / (0Y1*DY2) K 108 

C K 109 

UPP ( J)»( 0Y1*U1 t J + 1 )+DY2*Ul ( J-D-Ul ( J)*(DY1+0Y2) ) / ( . 5*DY1 *DY2* ( 0Y1+ K 110 
1DY2)) K 111 

WPP { J )-(DYl*Wl ( J + 1 )+DY2*Wl( J-D-Wl ( J)*(DY1+DY2) ) /( . 5*0 Y1*D Y2* ( 0 Y1+ K 112 
1DY2)) K 113 

TPP(J)«(DY1*T1( J + 1)+DY2 + T1( J-D-TK J)*(DY1+0Y2>)/(.5*DY1*DY2*(DY1+ K 114 
10Y2)) K 115 

CMUPP(J)-(DY1*CMU1( J+1)+DY2*CNU1( J-1)-CMU1(J)*(DY1+0Y2) )7( .5*0Y1*D K 116 
1Y2* (DY1+DY2 ) ) K 117 

10 CONTINUE K 118 

C K 119 

UP(NP)«(U1(NP)-U1(NP-1))/0Y2 K 120 

WP(NP)-(W1(NP)-W1(NP-1))/DY2 K 121 

TP(NP)«(T1(NP)-T1(NP-1) )/DY2 K 122 

: M'JP( NP) »(CNU1 (NP)-C MUK NP-1) ) /DY2 K 123 

ALFAP(NP)»(ALFA(NP)-ALFA(NP-1))/DY2 K 124 

C K 125 

UPP(NP)» (U1 (NP-1 J-UKNP) ) /0Y2**2 K 126 

WPP (NP)»( W1 (NP-l)-Wl(NP) ) /0Y2**2 K 127 

TPP (NP>« (T1 (NP-l)-Tl (NP) ) /DY2**2 K 128 

CMUPP(NP)»(CHU1(NP-1)-CMU1(NP))/DY2**2 K 129 

C K 130 

DY1»Y1(2)-Y1(1) K 131 

DY2»Y1(3)-Y1(1) K 132 

UP(1)-(DY2*U1(2)/DY1-DY1*U1(3)/DY2)M0Y2-0Y1)-U1(1)*(DY1 + DY2)/(DY1 K 133 
1*3Y2) K 134 

WP(1)» (DY2*W1 (2) /0Y1-DY1*W1 (3 )/DY2 )/(DY2-DYl)-Wl(l)*(DYl + DY2) /( DYl K 135 



♦ TIDY* 

S SUBROUTINE OUTPUT 

1*DY2) K 135 

TP(1)«(DY2*T1(2 )/0Yl-DYl*Tl(3)/DY2)/(0Y2-DYl)-Tm)*(DYl + DY2) / (DYl K 137 

1*0Y2) K 138 

CMUP(1)«{DY2*CMU1(2)/DY1-DY1*CMU1(3)/DY2) M 0Y2-DY1 ) -CMUl ( 1 )♦ ( DYl+D K 139 

1Y2) /(DY1+DY2) K l<rO 

ALFAP(l) «{DY2*ALFA(2) /0Y1-DY1*ALFA(3)/0Y2)/ (DY2-DY1)-ALFA(D* (DY1+ K lAl 

1DY2)/ (DY1*DY2) K 1<»2 

C K 1<^3 

X1»Y1(2) K lAA 

X2-YK3) K 145 

X3«Y1(4) K 146 

X4«Y1(5> K 147 

B2»X2*X3+X2*X4+X3*X4 K 148 

DL1»X1*(X1-X2)*(X1-X3)*(X1-X4) K 149 

C2«X1*X3+X1*X4+X3*X4 K 150 

0L2«X2*(X2-X1)*(X2-X3)*(X2-X4) K 151 

D2-X1*X2+X1*X4+X2*X4 K 152 

DL3«X3*(X3-X1)*(X3-X2)*(X3-X4) K 153 

E2*X1*X2+X1+X3+X2*X3 K 154 

DL4«X4*(X4-X1)*{X4-X2)*(X4-X3) K 155 

UPP(1)«2.*B2*U1 (2)/DLl+2.*C2*Ul(3)/DL2+2,*02*Ul(4)/0L3+2.*E2*Ul (5) K 156 

1/DL4 K 157 

WPP(1)»2.*B2*W1(2) /DL1+2.*C2*W1 (3 )/DL2+2.*D2*Ml( 4) /DL3+2.*E2*W1 ( 5) K 158 

1/DL4 K 159 

TPP(1)»2.*B2*T1{2)/0L1+2.*C2*T1{3)/DL?+2.*02*T1(4)/DL3+2.*E2*T1(5) K 160 

1/0L4 K 161 

CMUPP(1)«2.*B2*CMU1 (2 )/DLl+2.*C2*CMUl{3)<0L2+2.*02*CNUl(4)/0L3+2.* K 162 

1E2 + CMUK5 )/0L4 K 163 

C K 164 

TPP(1)»TPP(2) K 165 

CMUPP(1)-CMUPP(2) K 166 

C K 167 

DO 11 J«1#NP K 168 

I-NP-J+1 K 169 

AKD-YKJ) K 170 

A2(I)«U1{J) K 171 

A3( I)-UP( J) K 172 

A4(I)«UPP<J) K 173 

A5(I)«W1(J) K 174 

A6( I)«WP( J) K 175 

A7(I )-WPP(J ) K 176 

A8(I)-T1(J) K 177 

A9(I)«TP(J) K 178 

A10(I)-TPP(J) K 179 

A11(I)«CMU1(J) K 180 



* T I D Y * 


SUBROUTINE OUTPUT 


A12(I)-CMUP( J) K 181 

A13(I )»CMUPP(J) K 182 

Al<f(I)»ALFA(J) K 183 

A15(I)«ALFAPU) K 184 

11 CONTINUE K 185 

CALL CFA (A2»A4»A5»A7f EPSICF, NP ) K 186 

C K 187 

C K 188 

WRITE (10) NZ»XC(NZ)*CliACH>R»CHLENTH»WE(NZ)*Al(l),PR»A5(l),NP K 189 

WRITE (10) (Aid), I-1,NP) K 190 

WRITE (10) ( A2(I),I«1,NP) K 191 

WRITE (10) (A3(I),I-1,NP) K 192 

WRITE (10) ( A4 (I),I«1,NP) K 193 

WRITE (10) (A5(I),I«1,NP) K 194 

WRITE (10) (A6(I),I»1,NP) K 195 

WRITE (10) (A7(I),I-1,NP) K 196 

WRITE (10) (A8(I),I«1,NP) K 197 

WRITE (10) (A9(I),I-1,NP) K 198 

WRITE (10) (AIO(I) ,I«1,NP) K 199 

WRITE (10) ( All (I ),I«1,NP) K ZOO 

WRITE (10) (A12(I),I»1,NP) K 201 

WRITE (10) (A13(I ),I«1,NP) K 202 

WRITE (10) ( A14(I),I»1,NP) K 203 

WRITE (10) (A15 (I ), I-l, NP) K 204 

WRITE (10) (A16(I),I«1,NP) K 205 

C K 206 

IF ( IPRINT.E0.2) GO TO 14 K 207 

WRITE (6,16) NZ,XC (NZ ),NP,CHL ENTH,R,EPSICF*57. 29577 K208 

WRITE (6,17) K 209 

WRITE (6,19) (I,A1(I),A2(1),A3(I),A4(I),A5(I), A6(I),A7(I),A8(I),A9 K 210 

1(1) , A10( I), I»1,NP,5) K 211 

WRITE (6,19) NP,A1(NP),A2(NP),A3(NP),A4(NP),A5(NP),A6(NP),A7(NP),A K 212 

18(NP), A9(NP), AIO(NP) K 213 

WRITE (6,18) K 214 

WRITE (6,20) (I,A1(I),A11(I),A12(I),A13(I),A14(I),A15(I),A16(I),I» K 215 

11,NP,5) K 216 

WRITE (6,20) NP,A1(NP), AIKNP), A12(NP),A13(NP),A14(NP),A15(NP),A16 K 217 

KNP) K 218 

60 TO 14 K 219 

12 CONTINUE K 220 

RA-SORT(RHOE(NZ)*CMUE(NZ)/UE(NZ)) K 221 

00 13 I«1,NP K 222 

13 RAA(I)»RA*G(I,2) K 223 

WRITE (9) NZ,XC(NZ),RR(NZ),Z(NZ),RHOE(NZ),CMUE(NZ),UE(NZ),WE(NZ),P K 224 

1E(NZ),DDW(NZ),R,CHLENTH,RH0W,VW,TE,X,UFS,PR,CMACH,NP K 229 


* T I D Y ♦ 

SUBROUTINE OUTPUT 

o 

o 

WRITE (9) {Y1(I)#I«1»NP) K 226 

WRITE (9) (U1 ( I )» I» 1» NP ) K 22T 

WRITE (9) (W1(I),I«1»NP) K 228 

WRITE (9) (TKI )»I*1»NP) K 229 

WRITE (9) (CMU1(I)»I«1»NP) K 230 

WRITE (9) ( ALFA (I ),I«1,NP) K 231 

WRITE (9) (F(I,2),I«1»NP) K 232 

WRITE (9) (U(I>2)»I«1,NP) K 233 

WRITE (9) (RAA( I) ,I«1,NP) K 23A 

WRITE 19) lDENRa,2)»I«l»NP) K 235 

C K 236 

CONTINUE K 237 

DO 15 J-1»NP K 238 

F( J,1)«F ( J,2) K 239 

U( J»l)-U( J,2) K Z<t0 

V( J.l)»V(Ji2) K 241 

G(J,1)-G (J,2) K 242 

W( J,1)»W(J,2) K 243 

T(Jf l)«T(Ji2) K 244 

B(J, 1)«3 (J>2) K 245 

DENR( J, 1)«DENR (J,2) K 246 

IF (CMACH.EQ.0.0) go to 15 K 247 

E( J»1)«E (J»2) K 248 

BG( J,1)-BG( J»2) K 249 

CA1(J,1)»CA1(J,2) K 250 

CA2{J>1)«CA2(J»2) K 251 

15 CONTINUE K 252 

C K 253 

NZ-NZ+1 K 254 

RETURN K 255 

C K 256 

C K 257 

C K 258 

16 FORMAT ( A//1X,3HNZ«» I5,3X,4HX/C«»E13.6»3X,3HNP»,I5»3X»8HCHLENTH«>E K 259 

116.6,3X,2HR-tE16.6^ 3X»4HCFA-»E16.6) K 260 

17 FORMAT (IH ,5X, 30H(Y,U,UP,UPP,W,WP»WPP,T,T!>,TPP)/) K 261 

18 FORMAT (IH , /» 5 X, 31 H ( MU» MUP » MUPP , AL F A» ALF AP» PR ANDU / ) K 262 

19 FORMAT (IH » 2 X , I4» 2 X, F6 . 3> 9 E 1 2. 4 ) K 263 

20 FORMAT (IH , 2X » 14, 2X, F6 . 3, 6E 1 2. 4 ) K 264 

END K 265- 



♦ TIDY* 

SUBROUTINE XZDER 


SUBROUTINE XZDER L 1 

REAL LST L 2 

DIMENSION X (3)» Xl(3)# S(3), LST(<»), R(3)» THETA(3)» NSA(A), RHOE( L 3 

13), CMUE(3), UE(3), WE(3), PE(3), DDW(3), RH0W(3), VW(3), TE(3), X L A 

21(3), UFS(3), NZN(3) L 5 

COMMON /BXT/ D6TH (101 ) , DEN (1 01 ) , YSTS ( 101 ) ,YS ( 101 ) , F ( 101 ) , U ( 101 ) , G ( L 6 

1101),V(101),VP(101),VPP(101),N L 7 

COMMON /BLCO/ NZ T, NA, NB , I T, C , D, C M AC H, TT, FTA ( 1 51 ) , DE T A ( 1 51 ) , A ( 15 1 ) , L 8 

lYK 151),IPANPA L 9 


DIMENSION AKIOI), A2(101), A3(101), A<*(101), A5(101), A6(101), A7 L 10 
1(101), ABdOl), A9(101), AlO(lOl), AlKlOl), A12(101), A13(101), A L 11 
21^(101), AlSdOl), A16(101), A17(101), A18(101), A19(101), A20(101 L 12 
3), A2K101), A22(101), A23(101), A2^(101), A25(101), A26(101), A27 L 13 
4(101), A28(101), A29(101), A30(101), A3K101), A32(101), A33(101) L 14 
DIMENSION DUX(lOl), DUPX(lOl), DUPPX(lOl), OHX(lOl), OWPX(lOl), DW L 15 
IPPX(lOl), DTX(lOl), DTPX(lOl), DTPPX(lOl), OMUX(lOl), DMUPX(lOl), L 16 


2DMUFPX(101), DALFX(lOl), DALFPX(lOl), PRANDL(lOl) L 17 

DIMENSION UP0(3,101), UPP0(3,101), WP0(3,101), WPP0(3,101), TP0(3, L 18 

1101), TPP0(3,101), CMUPO( 3, 101), C MUPPO ( 3, 1 01 ) , ALFAPO ( 3, 101) , YST L 19 

2(4,101) L 20 

DIMENSION Y(4,101), U0(4,101), W0(4,101), T0(4,101), CMU0(4,101), L 21 

lALF A0(4,101 ), FA0(4,101), GA0(4,101), UA0(4,101), DENR AO ( 4, 101 ) L 22 

COMMON /TAXI/ Dl, D N, D3, D4 , D 5, D6, 07, D8 , D9, D1 0, Dll, D 1 2 L 23 

COMMON /PAR/ AAl, AA2, AA3, VGP, NNN,IPRINT L 24 

C L 25 

REWIND 9 L 26 

DO 1 I«l,3 L 27 

READ (9) NZN(I),X(I),S(I),THETA(I),RHOE(I),CMUE(I),UE(I),WE(I ),P£( L 28 

II), DDW(I)>R(I)» LST(I),RHOW(I ) ,VW ( I ) , TE ( I ) » X I ( I ) , UFS ( I ) , PR ,C MACH, NS L 29 

2A(I) L 30 

NP«NSA(I) L 31 

READ (9) (Y(I,J),J-1,NP) L 32 

READ (9) (UO(I, J ), J«1,NP) L 33 

READ (9) (WO(I,J)» J«1,NP) L 34 

READ (9) (TOd, J),J«1,NP) L 35 

READ (9) (CMUO(I,J), J«1,NP) L 36 

READ (9) (ALFAO (I,J ), J-1, NP) L 37 

READ (9) (FA0(I,J),J«1,NP) L 38 

READ (9) (UAOd, J ), J»1,NP) L 39 

READ (9) (GAO(I,J), J«1,NP ) L 40 

READ (9) (DENRAOd , J ) , J» 1,NP ) L 41 

C WRITE(6,31) ( Y(I, J),FAO(I, J),UAO(I,J ),GAO(I, J),J»1,NP) L 42 

C WRITE(6,13)NZN(I),Xd),S(I),THETA(I),RH0Ed),CMUEd),UE(I), L 43 

C *WE(I),PE(I),ODWd),Rd),LST(I ),RHOW(n,VW(I),TE(n,XI(I), L 44 

C *UFS(I),PR,CMACH,NSA(I) L 45 
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♦ TIDY’* 


SUBROUTINE YZOE"? 


C WRITE(6»12) (Y( I, J),UO( I» J )»WO(I»J TO(I» J),CMUO(I, J ), ALFAOd, 

C *FAO( J )»UAO(I> J))GAO(If J )>DENRAO(I f J ), J-1,NP) 

1 CONTINUE 
K»A 

2 CONTINUE 
NP-NSA (3 ) 

NSA (A)»NSA( 3) 

LST(4)»LST(3) 

DO 3 J-1>NP 
Y(A,J )-Y (3»J ) 

UO(A» J)»UO( 3» J ) 

W0(4»J)>W0(3«J ) 

T0( A# J )-T0( 3,J ) 

CMUO(A» J)«CMU0(3»J ) 

ALFAO(-V^ J)-ALFA0(3» J) 

FA0(4,J )«FA0(3>J) 

UAO(A, J)«UA0(3» J) 

GA0(A>J)-GA0(3» J) 

3 DENRA0(A»J)»DENRA0(3»J) 

DO A I«1,A 

NP-NSA(I ) 

DO A J*1#NP 

A YST(I» J)«Y( I»J )*LST(I ) 

N«NSA(2) 

DO 5 I«l»3»2 
M«NSA( I) 

CALL PROFO ( YST,U0»I»2»N,M) 

CALL PROFO (YST,WO, I»2,N,M) 

CALL PROFO (YSTdOfl »2#N,M) 

CALL PROFO (YST,CMUOf I»2»N»M) 

CALL PROFO (YST, ALFAO, I» 2, N, M) 

CALL PROFO ( YST»FA0»I.2»N,M) 

CALL PROFO { YST^UA0,I»2,N,M) 

CALL PROFO (YST,GA0,I»2,N,M) 

CALL PROFO (YST,DENRA0,I»2»N,M) 

5 CONTINUE 

C 

DO 6 I»1»N 
YSTS(I)»YST(2,I ) 

YS( I)»Y(2,I ) 

DEN(I)«DENRA0(2»I) 

F( n-FA0( 2»I ) 

U(I)«UA0{2,I) 

6 G(I)-GA0(2i I) 

DO 0 1-1,3 


L A6 
L A7 
L AB 
L A<? 



L 75 
L 76 
L 77 
L 78 
L 79 
L 80 
L 81 
L 82 
L 83 


L 8A 
L 85 
L 86 
L 87 
L 88 
L 89 
L 90 
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♦ TIDY* 

SUBROUTINE XZOER 

NPMl-N-1 L 91 

DO 7 J«2»NPM1 L 92 

DY1»Y{2, J)-Y(2,J-1) L 93 

0Y2»Y(2j J+1)-Y(2#J) L 9^ 

UPO{ I, J ) - (OY1*UO(I » J + 1 ) /DY2-DY2tU0n» J-l)/DYl)/( DY1+DY2)-U0(I» J )♦( L 95 
1DY1-DY2)/(DY1*DY2) L 96 

WPO( J )■ (DY1 + W0( I» J+1) /DY2-DY2 + W0 ( I » J-1 ) /O Y1 ) / ( DY 1+0 Y2 ) -WO ( I # J )* ( L 97 
1DY1-0Y2)/(DY1*0Y2) L 98 

TPOdf J)-(DY1*T0(I» J + 1) /DY2-DY2*T0( I» J-1 ) /O Y1 ) / ( D Y 1+OY 2 )-T0 ( I »J)*( L 99 
1DY1-0Y2)MDY1*DY2) L 100 

CMUPO(I» J)» (DYl*CMUO( I> J+l)/DY2-DY2*CNU0n, J-D/DYl )/ (0Y1 + DY2)-CMU L 101 
10(I,J)*(DY1-DY2)/(0Y1*DY2) L 102 

ALFAP0(I»J)«(DYI*ALFA0{I» J+1 )/DY2-DY2tALPA0(I,J-l) /DYl ) / ( DY1+ DY2)- L 103 
1ALFA0(I,J)*(DY1-DY2)/(DY1*DY2) L 10^ 

C L 105 

UPPO( I» J)«(0Y1*U0{ If J + 1)+0Y2+U0{ I, J-1 )-U0 ( I , J > ( D Y 1 + OY 2 ) ) / ( .5+DY1+ L 106 
10Y2*(DY1+DY2 ) ) L 107 

WPPO(I»J )«(DY1*W0( I» J + 1) +0Y2 + W0 (I, J-1)-W0(I, J )*(DY1+0Y2) ) /( ,5*DY1* L 108 
1DY2*(DY1+0Y2)) L 109 

TPPO (I» J )«(DY1*T0(I» J+1)+DY2*T0(I, J-1 )-TO(I,J )*(0Y1+DY2) ) /( .5*DY1* L 110 
1DY2*(DY1+DY2) ) L 111 

CMUFPO(I» J)-(0Y1*CMU0(I, J + 1)+DY2*CMU0(I» J-1)-CMU0( J ) * ( DY1+0Y2 ) )/ L 112 
1(.5*0Y1*DY2*(DY1+0Y2) ) L 113 

7 CONTINUE L ll<^ 

C L 115 

UP0(I»N)»{U0(I»N)-U0(I»N-1) ) /0Y2 L 116 

WP0(I»N)-(W0(I,N)-W0( I»N-1) ) /DY2 L 117 

TP0n»N)»(T0<I»N)-T0(I,N-l))/DY2 L 118 

CNUPO(I» N)«(CMU0(I»N)-CriU0( I, N-1) )/DY? L 119 

ALFAPO(I»N)«(ALFAO(I,N)-ALFAO(dN-l ))/DY2 L 120 

C L 121 

UPPO(I,N)-(UO(I»N-l)-UOn»N) )/DY2**2 L 122 

WPP0(I,N)«(W0(I,N-1)-W0( I»N) )/DY2* + 2 L 123 

TPP0(I*N)«(T0{I.N-1)-T0(I»N) ) /DY2**2 L 12A 

CMUPPO(I»N)»(CMUO(I»N-1)-CMUO(I»N))/OY2**2 L 125 

C L 126 

DY1-Y(2»2)-Y(2»1) L 127 

DY2»Y(2,3)-Y(2^1) L 128 

UP0(I,1)« (DY2*U0(I»2) 7DY1-0Y1*U0(I#3)/0Y2)/ (DY2-0Y1)-U0(I,1 )♦ (DY1+ L 129 
1DY2)/(DY1*DY2) L 130 

WPO( I>1)« (DY2*W0(I,2) /DY1-DY1*W0(I,3)7DY2)/ (DY2-0Y1 )-W0(I,l)*(0Yl+ L 131 
1DY2)/(DY1*DY2) L 132 

TPOd, l)»(DY2*T0{If 2) /DY1-DY1*T0( I, 3) 70 Y2 ) / (D Y2-DY1 )-TOd , 1 )* (DY1+ L 133 
1DY2)/(DY1*DY2) L 134 

CMUPO( dl)« (DY2t'CMU0(I,2)/DYl-DYl*CMU0d»3) 7DY2) /( DY2-DY1 )-CMUO d, L 135 



o r> o <3> 


o 


♦ TIDY* 


SUBROUTINE XZOER 


11)*(DY1+DY2)/(DY1*DY2) L 136 

ALFAP0(I»1)«(DY2*ALFA0(I,2) /DY1-DY1+ALFA0(I»3)/DY2)/(DY2-DY1)-ALFA L 137 

10(I,1)*(DY1+DY2)/(0Y1*DY2) L 138 

C L 139 

X11«Y{2»2) L 140 

X2-Y(2,3) L 141 

X3»Y(2>4) L 142 

X4-Y(2»5) L 143 

B2« X2fX3+X2*X4+X3*X4 L 144 

DL1«X11*(X11-X2)*{X11-X3)*(X11-X4) L 145 

C2«X11*X3+X11*X4+X3*X4 L 146 

DL2»X2*( X2-X11)*(X2-X3)*(X2-X4) L 147 

D2« X11*X2+X11*X4+X2*X4 L 148 

0L3»X3*( X3-X11 )*(X3-X2)*(X3-X4) L 149 

E2- X11*X2+X11+X3+X2*X3 L 150 

DL4»X4^(X4-X11)*(X4-X2)*(X4-X3) L 151 

UI»PO( I>1)»2.*B2*U0(I»2 )/DLl + 2.*C2*U0(I»3) /DL2 +2 . *02*U0 ( I » 4 ) /D L 3 + 2. L 152 

l’*‘E2*U0(I,5)/DL4 L 153 

WPPOCIf 1)«2.+B2*W0( I#2 )/DLl+2.*C2*W0( 1,3) /0L2+2.*02*W0(I»4)/DL3+2. L 154 

l*E2*W0n,5)/DL4 L 155 

TPP0(I,1 )«2.*B2*T0( I,2)/DL1 + 2.*'C2*T0(I,3) /OL 2+2. *D2*T0 ( I, 4 ) /DL3 + 2. L 156 

1*E2+T0(I,5)/DL4 L 157 

CMUPP0(I,1)«2.*B2*CMU0(I,2) /OL1+2.*C2^CMUO( I,3)/DL2+?.*02*CMU0(I,4 L 158 

1)/OL3+2.*E2+CMUO(I,5)/DL4 L 159 

C L 160 

TPP0(I,1)-TPP0(I,2) L 161 

CMUPP0(I,1)»CMUPP0(I,2) L 162 

C L 163 

CONTINUE L 164 

WRITE(6,130) (Y(l,n,UO{l,I ),U0(2,I ),UP0(1,T), L 165 

*UPO (2, I ), I»1,N) L 166 

L 167 

DO 9 1*1,3 L 168 

9 Xl( I)«S(I)/(LST (2)*R(2)) L 169 

DX1»X1(2)-X1 (1 ) L 170 

DX2»X1(3)-X1( 2) L 171 

DTH1»THETA(2)-THETA(1) L 172 

DTH2»THETA(3)-THETA(2) L 173 

0ELX«DX1/DX2 L 174 

DELTH-DTH1/DTH2 L 175 

DELX1-DX1+DX2 L 176 

DELTH1-DTH1+DTH2 L 177 

DELX2»(DX1-DX2 )/(DXl*0X2) L 178 

DELTH2-(DTH1-DTH2 )Z (DTH1+DTH2) L 179 

no 10 1*1, N L 180 



105 


* T I D Y ♦ 

SUBROUTINE XZDER 

DUX(I)«{(DELX*U0(3,I)-U0(l,I)/DELX)/nELXl>0ELX2*U0(2»I)) L 181 

DUPXd )»( (DELX + UP 0(3,1 )-UPO( 1,1 ) /DELX ) /DELX1-0ELX2*UP0( 2, I ) ) L 182 

DUPPX( I) » ( ( DELX+UPPO (3,n-UPP0( 1, I ) /OELX) rOE LXl-OE LX 2 *UI»pO ( 2, 1) ) L 18 3 

DWX(I)»((0ELX*W0{3,I)-W0(1,I )/DELX)/DELXl-OELX2*WO(2,I)) L 184 

DWPX( I )» ( (DELX+WPO( 3,1 )-WPO( 1,1 ) /OELX) /DELX1-DFLX2*WP0( 2, I ) ) L 185 

DWPPX(I)» ((OELXi'WPPOO,! )-WPPO(l, I)/DELX)/DELX1-DELX2*WPP0(2, I ) ) L 186 

OTX(I)-( (DELX*T0(3,I )-TO(l,I )/DELX)/0ELXl-0ELX2 + T0(2,D) L 187 

DTPX{I)» ( (DELX*TP0(3,I )-TPO( 1,1) /DELX) /DELX1-0ELX2*TP0( 2, I ) ) L 188 

DTPPX( I)« ( (DELX+TPPO (3, I)-TPPO( 1, I ) / DE LX ) /OELX 1-DELX 2 + T P PO ( 2 , I) ) L 189 

DMUX(I)«( (DELX*CMU0(3,I )-CMUO(l,I )/DELX)/OELXl-DELX2*CMUO(2,I )) L 190 

DMUPX(I)- ((DELX»CMUP0(3,I ) -C MUP 0 ( 1 , I ) /O EL X ) T OE LX 1-D EL X 2*C MUPO ( 2, I ) L 191 

1) L 192 

DMUPPXd )«( (DELX>t'CMUPP0(3,I)-CMUPP0d,I)/DELX)/DELXl-0ELX2*CMUPP0( L 193 

12»D) L 194 

DALFXd)«((DELX4ALFAO(3,I)-ALFAO(l, n/OELX) ZDELXl-r)ELX2*ALFA0( 2, I ) L 195 

1) L 196 

DALFPX(I)»( (DELX*ALFAPO(3,n-ALFAPOd,I)/OFLX)/OELXl-DELX2*ALFAPO( L 197 

12, I>) L 198 

10 CONTINUE L 199 

C L 200 

DO 11 I»1,N L 201 

DGTH(I)« (DELTH*GA0(3,I)-GA0( 1,1) /DELTH)/0FLTH1-DELTH2*GA0(2,I ) L 202 

C WRITF(6,31) Fd ),Ud ),G(I),DGTHd) L 203 

11 CONTINUE L 204 

C L 205 

D1»RH0W(2) L 206 

DN«VW(2) L 207 

D3-UE(2) L 208 

D4»PE(2) L 209 

D5-XK2) L 210 

D6-WE(2) L 211 

D7-DDW(2) L 212 

D8»TE(2) L 213 

D9-CMUE(2) L 214 

D10«RH0E(2) L 215 

D11«UFS(2) L 216 

D12»R(2) L 217 

CALL VVEL L 218 

DO 12 I«1,N L 219 

12 PRANDLd)»PR L 220 

DO 13 J»1,N L 221 

I*N“J+1 L 222 

A1(I)«Y(2,J) L 223 

A2(I)-U0(2, J) L 224 

A3(I)»UP0(2,J) L 225 
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♦ TIDY* 


SUBROUTINE XZDEP 


A<t( I )-UPP0(2» J ) L 226 

A5m-WO(2>J ) L 227 

A6(I)-WP0(2, J) L 228 

A7{I)«WPP0(2»J ) L 229 

A8m»T0(2,J) L 230 

A9{ n«TP0(2, J) L 231 

A10(I)-TPP0(2, J) L 232 

All(I)-CriU0(2»J) L 233 

A12 ( n«CMUP0(2f J ) L 23A 

A13(I)-CMUPP0(2» J) L 235 

AlAm«ALEA0(2»J) L 236 

A15( I) -ALFAP0(2, J ) L 237 

A16(I)-PRANDL(J ) L 238 

A17(I)«DUX(J) L 239 

A18(I)-0UPX( J) L 2^0 

A19(I )»DUPPX ( J ) L 241 

A20(I )»DWX(J ) L 242 

A21{I)«DWPX( J) L 243 

A22(I)»DWPPX(J) L 244 

A23(I )-DTX( J ) L 245 

A24(I)-DTPX(J) L 246 

A25(I) «OTPPX(J ) L 247 

A26( I )»DMUX (J ) L 248 

A27(I)«0MUPX(J ) L 249 

A28 (I )»DMUPPX( J ) L 250 

A29(I )«DALFX(J ) L 251 

A30(I)»0ALFPX(J) L 252 

A31(I)»V(J) L 253 

A32(I )«VP(J ) L 254 

A33 ( I )-VPP( J) L 255 

13 CONTINUE L 256 

CALL CFA (A2»A4»A5iA7,EPSICF,N) L 257 

WRITE (10) NZN(2),X(2),CMACH,R(2)»LST(2),UFS(2),Y(2,N),PR,W0(2»N), L 258 

IN L 259 

WRITE (10) (Al( I),I*1,N) L 260 

WRITE (10) (A2( I)»I«1»N) L 261 

WRITE (10) (A3(I),I-1,N) L 262 

WRITE (10) (A4(I),I«1»N) L 263 

WRITE (10) (A5(I),I-1,N) L 264 

WRITE (10) (A6(I ),I«1,N) L 265 

WRITE (10) (A7(I)»I-1,N) L 266 

WRITE (10) (A8 (I )> I«1»N) L 267 

WRITE (10) (A9(I),I-1,N) L 268 

WRITE (10) ( A10(I),I»1»N) L 269 

WRITE (10) (A11(I)»I«1»N) L 270 
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♦ TIDY* 

SUBROUTINE XZOER 

WRITE (10) (A12 (I)>I«1^N) L 271 

WRITE (10) (A13(I)> I«1»N) L 272 

WRITE (10) (A1A(I)^I»1#N) L 273 

WRITE (10) (A15(I)»I-1»N) L 274 

WRITE (10) (A16(I)»I«1»N) L 275 

C L 276 

WRITE (11) NZN(2),X(2),N L 277 

WRITE (11) (A17(I )»I«1»N) L 278 

WRITE (11) (A18 (I )»I«1>N) L 279 

WRITE (11) (A19(I ),I«1,N) L 280 

WRITE (11) (A20(I)»I-1,N) L 281 

WRITE (11) (A21(I)»I«1»N) L 282 

WRITE (11) ( A22 (I )# I-1»N) L 283 

WRITE (11) (A23 (I)»I«1>N) L 284 

WRITE (11) (A24(I)»I»1,N) L 285 

WRITE (11) (A25(I ),I»1»N) L 286 

WRITE (11) ( A26(I),I«1,N) L 287 

WRITE (11) (A27(I)#I»1,N) L 288 

WRITE (11) (A28(I)»I«1#N) L 289 

WRITE (11) (A29(I)»I-1,N) L 290 

WRITE (11) (A30(I )»I«1#N) L 291 

WRITE (11) (A31 (I ) »I»1»N) L 292 

WRITE (11) (A32(I)»I«1^N) L 293 

WRITE (11) (A33(I)f I-1»N) L 294 

C L 295 

IF (IPRINT.E0.2) GO TO 14 L 296 

WRITE (6,20) NZN(2),X(2),N,LST(2),R(2),ERSICF L 297 

WRITE (6,21) L 298 

WRITE (6,25) (I,Y(2,I ),U0(2, I),UP0(2,I),'JPP0(2,I),W0(2,I),WP0(2,I) L 299 

1,WPP0(2,I),T0(2,I),TP0(2,I),TPP0(2,I),I«1,N,5) L 300 

WRITE (6,25) N,Y(2,N),UO(2,N),UPO(2,N),UPPO(2,N),WO(2,N),WPO(2,N), L 301 
1WPP0(2,N),T0(2,N),TP0(2,N),TPP0(2,N) L 302 

WRITE (6,22) L 303 

WRITE (6,25) (I,Y(2,I),CMU0(2,I),CHUP0(2,I),CMUPP0(2,I),ALFA0(2,I) L 304 

1 ,ALFAP 0 ( 2 ,I),PRAN 0 L(I),DUX(I ),DUPX(I),DU'»PX(n,I«l,N,5) L 305 

WRITE (6,25) N,Y(2,N),CNUO(2,N),CMUPO(2,N),CMUPPO(2,N),ALFAO(2,N), L 306 
1ALFAPO(2,N),PRANDL(N),DUX(N),DUPX(N),OUPPX(N) L 307 

WRITE (6,23) L 308 

WRITE (6,25) (I,Y(2,I),DWX(I),DWPX(I),0WPPX(I),DTX(I),DTPX(I),DTPP L 309 

1X(I ),DMUX(I),DMUPX(I),DMUPPX(I),I-1,N,5) L 310 

WRITE (6,25) N,Y(2,N),DWX(N),DWPX(N),f)WPPX(N),0TX(N),0TPX(N),DTPPX L 311 
1(N),DHUX(N),0MUPX(N),DMUPPX(N) L 312 

WRITE (6,24) L 313 

WRITE (6,26) (I,Y(2,I),DALFX(I),0ALFPX(I),V(I),VP(I),VPP(I),I»1,N, L 314 

15) L 315 
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♦TIDY* 


SUBROUTINE XZDER 


WRITE (6,26) N,Y(2,N),DALFX( N),DALFPX(N),V(N), VP(N), VPP(N) 
lA CONTINUE 
C 

DO 15 1-1,2 
M-I + 1 

NZN(I)«NZN(M) 

X(I)-X(M) 

S( I )-S(M) 

THETA(I)»THETA{f1) 

RHOE (I )-RHOE(M) 

CMUE(I)-CMUE(M) 

UE(I)-UE(M) 

WE (I )»WE (M) 

PE ( D-PE (M) 

OOW(I)»DDW(M) 

Rn)«R(fi) 

LST(I)-LST(M) 

RHDW(I)-RHOW(M) 

VW(I )-VW (M) 

TE (I )-TE (M) 

XI(I)-XKM) 

15 UFSd )»UFS(M) 

C 

00 18 1-1,2 
IF (I.E0.2) GO TO 16 
M-I + 1 
GO TO 17 

16 H-I+2 

17 CONTINUE 
NSAd )-NSA(N) 

NP-NSA (M ) 

DO 18 J-1,NP 
Y( I, J) -Y(M, J ) 

UO(I,J)-UO(M,J ) 

WO ( I,J )-WO(M,J ) 

Ton, J)»TO(N,J ) 

CMUO( I,J )-CMUO (M, J ) 

ALFAO(I,J )-ALFAO(M,J) 

FAOd, J)-FAO(M,J) 

UAOd,J) -UAO(M, J ) 

GAOd,J)»GAO(M,J) 

18 DENRAOd, J)-DENRAO(M, J) 

C 

1-3 

K-K+1 


L 316 
L 317 
L 318 
L 319 
L 320 
L 321 
L 322 
L 323 
L 32A 
L 325 
L 326 
L 327 
L 328 
L 329 
L 330 
L 331 
L 332 
L 333 
L 334 
L 335 
L 336 
L 337 
L 338 
L 339 
L 340 
L 341 
L 342 
L 343 
L 344 
L 345 
L 346 
L 347 
L 348 
L 349 
L 350 
L 351 
L 352 
L 353 
L 354 
L 355 
L 356 
L 357 
L 358 
L 359 
L 360 



♦ TIDY* 

SUBROUTINE XZDER 


IF (K.GT.NZT) GO TO 19 L 361 

READ (9) NZN(I)jX(I)»S(I)>THETA(I),RHOE{I)»CMUEm»UE(I)>WE(I),PE( L 362 

lI)#DDVMI)»Rn),LST(I)>RHOwm>VWn) »TE{I)»Xin)»UFS(I)»PR»CMACH»NS L 363 

2A(I ) L 36A 

NP-NSA(I) L 365 

READ (9) (Y(I,J ),J-1»NP) L 366 

READ (9) (UO(I, J), J«1,NP) L 367 

READ (9) (WO(I> J)»J«1^NP) L 360 

READ (9) (TO(I» J)»J-1»NP) L 369 

READ (9) (CMUOdfJ ),J«1,NP) L 370 

READ (9) (ALFAO(I>J)» J«1,NP) L 371 

READ (9) (FA0(I,J),J«1>NP) L 372 

READ (9) (UAO( I»J)# J«1>NP) L 373 

READ (9) (GAO(I» J)» J»1,NP ) L 37<* 

READ (9) (DENRAO(IfJ),J»l#NP) L 375 

C DO 11 I«l,3 L 376 

C NP«NSA(I) L 377 

Cll WRITE(6,3inYn,J),FA0(I,J).UA0(I,J),GA0(I,J),J-l,NP) L 378 

C WRITE (6> 13) NZN(I )» X(I ),S (I THETA( I ),RHOE{I ),CNUE( I).UF( I), L 379 

C *WE(I)iPE(I)»ODWn)»Rm»LSTm»RHOWm»VW(I)»TE(I)»XI(I)# L 380 

C *UFS(I)»PR»CMACH,NSA(I) L 381 

C WRITE(6»12) (Yn,J)»UO(I» J)»WOn»J)>TO(l,J)»CMUOn»J)»ALFAO(I»J)» L 382 

C *FAO(I^ J)»UAO(I»J)fGAO(I»J)»DENRAO(I»J)»J*l»NP) L 383 

GO TO 2 L 38A 

19 RETURN L 385 

C L 386 

C L 387 

20 FORMAT ( / 7/ IX, 3HNZ ■* I 5, 3X, X ZC E 16.6» 3X,3HNP», 1 5» 3X, 8HCHLEN TH»> E L 388 

116,6,3X,2HR«,E16.6,3X,AHCFA-,E16.6) L 389 

21 FORMAT (IH > 5X , 28HY, U, UP, UPP , W, W P , W PP» T , T P, TPP / ) L 390 

22 FORMAT (IH , 5X, <»2HY, MU, MUP, MUPP, ALF A, ALFA P, PR, DUX, DUPX, DUPPX/ ) L 391 

23 FORMAT ( IH , 5X , A7HY, OWX , D WP X , DWPP X, DT X, DT PX , DTPP X, DMUX, DMP X, DMPPX/ L 392 

1) L 393 

2<* FORMAT (IH , 5X, 23HY, ALFAX, ALF APX, V, VP, VP?/) L 39A 

25 FORMAT (IH , 2X, lA, 2X, F6 . 3, 9E 1 2. A ) L 395 

26 FORMAT ( IH , 2X, I A, 2 X, F6 . 3 , 5E 1 2. A ) L 396 

END L 397- 


o 

VO 
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♦TIDY* 

SUBROUTINE PROFO (Y,U»I »K»N»M ) 

SUBROUTINE PROFO ( Yt U» I » K# N, M ) M 1 

REAL INTER M 2 

DIMENSION Y(4,101), U(A»101)> Yl(101)» Y2(101), U2(101), U3(101) M 3 

DO 1 J«1»M M A 

U2(J)-U(I»J) M 5 

Y2(J)-Y(I,J) M 6 

1 CONTINUE M 7 

DO 2 J«1»N M 8 

2 Yl( Jl-YCKfU ) M 9 

IF (Y2(M) .GE.Yl (N) ) GO TO 3 M 10 

L«M+1 M 11 

Y2(L)-Y1(N) M 12 

U2(L)»U2(M) M 13 

GO TO ^ M l<f 

3 CONTINUE M 15 

L«M M 16 

CONTINUE M 17 

DO R JJJ«1»N M 18 

YARG-YKJJJ) M IP 

DO 5 J-1»L M 20 

JJ»J M 21 

IF (YARG.LT.Y2( J) ) GO TO 6 M 22 

IF (YARG.EQ.Y2( J) ) GO TO 7 M 23 

5 CONTINUE M 2A 

6 MIN-JJ-3 M 25 

IF (JJ.LE.3) MIN«1 M 26 

IF ( JJ.GE. (L-2 ) ) MIN-L-6 M 27 

U3( JJJ)«INTER(Y2#U2,YARG»6»MIN) M 28 

GO TO 8 M 29 

7 U3( JJJ)-U2(JJ) M 30 

8 CONTINUE M 31 

DO 9 J«1»N M 32 

U(I,J)-U3(J) M 33 

9 CONTINUE M 3A 

RETURN M 35 

end M 36- 



Ill 


♦ TIDY* 

REAL FUNCTIQNINTER(X,Y. XARG, IDEG,MTN) 

REAL FUNCTIONINTER(X»Y»XARG>IDEG#MIN) N 1 

DIMENSION XilOl), Y(lOl) N ? 

1 FACTOR-1. N 3 

HAX-MIN+IOEG N A 

DO 2 J«MIN»MAX N 5 

IF (XARG.NE.X(J)) GO TO 2 N 6 

INTER«Y(J) N 7 

RETURN N 8 

2 FACTOR«FACTOR*(XARG-X(J) ) N 9 

YEST-0. N 10 

DO A I-MIN»MAX N 11 

TERM«Y(I)*FACTOR/(XARG-X(I)) N 12 

DO 3 J»MIN,MAX N 13 

3 IF (I.NE.J) TERM»TERM/(X(I )-X(J ) ) N lA 

A YEST-TERM+YEST N 15 

INTER-YEST N 16 

RETURN N 17 

END N 18- 
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♦ TIDY* 

SUBROUTINE CFA ( U1,U1PP»W1,W1PP»EPSICF»N) 


SUBROUTINE CFA ( Ulf UIPP, Wl, W 1 PP^ E PS ICF, N ) □ 1 

DIMENSION U(lOl), W(lOl), UPP(lOl), WPP(101)» Ul(10l)» Wl(lOl)# U1 0 2 

1PP(101J» WlPP(lOl) Q 3 

DO 1 J«1»N 0 A 

I»N-J+1 0 5 

U(I)»U1(J) 0 6 

W(I)«W1(J) G 7 

UPP(I)«U1PP( J) 0 n 

WPP ( I )«W1PP ( J ) Q 9 

1 CONTINUE 0 10 

Bl»l.E+06 0 11 

DO 2 I»2»N 0 12 

IF (U( I) .GT, 0.999) GO TO 3 0 13 

A1*U(I)/W(I) 0 

A2«UPP(I)/WPP(I) 0 15 

B2«A1-A2 0 16 

IF ((B2.GE.O..AND.B1.LE.O.).OR.(B2.LE.O..ANO.B1.GE.O.)) ISAVE-I 0 17 

61«B2 0 18 

2 CONTINUE □ 19 

3 CONTINUE 0 20 

EPSICF— ATAN(U(ISAVE)/W(ISAVE ))+3.1A15977 0 21 

C WRITE (6#4) EPSICFt57. 29577 0 22 

RETURN 0 23 

C 0 2A 

C 0 25 

END 0 26“ 
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♦ TIDY* 


SUBROUTINE VVEL 


C 

C 

C30 


C 


1 


c 

c 

2 

c 


3 


SUBROUTINE VVEL 
DIMENSION ETA(lOl) 

COMMON /TAXI^ RHOW» VW» UE» PE# X,WE> DOW, TE>CMUE» RHOE» UFS» R 
C OMMON /BXT/ DGTH ( 101 ), D E N( 1 01) ^ YS T { 1 01 ) , Y ( 10 1 ) » F A ( 1 01 ) , UA ( 101 ) » GA 
1(101),V{101)» VP(101)» VPP(101)»N 

WRITE (6# 30) RHOWf VW»UEi PE»X, WEj DOW»TE>CMUE» RHOE»UFS 
FORMATdH > 5X# IKEIO. > ) 

SQRX«SQRT(RHOE*UE*X/CMUE) 

PAR3«SQRX/X 

SUM«0. 

F1»1/DEN(1) 

ETA(1)«0. 

DO 1 J«2>N 
F2-1/DEN(J) 

SUM-SUM+(F1+F?)*(YST( J)-YST( J-1) ) / 2. 

F1-F2 

ETA(J)«SUM*PAR3 
CONTINUE 
CON»RHOW*VW/WE 
DO 2 I«1>N 

RH0«PE/(1716,*TE*DEN(I) ) 

DPSIX« (CMUE*SQRX/2. )*( 3.*FA( I )-ETA( I)*UA{I) ) 

DPHITH«SQRT( X**3)* (WE*0GTH(I ) +GA ( I ) ♦DDW*UFS ) 

V( I )«CON/RHO+ (1 ./ (RHO*X*WE ) ) *(DPS I X-DPHITH/X ) 

V(I )»V(I )*R 

WRITE! 6»31)Y(I ) , F A ( I ) , UA ( I ) » G A ( I ) » DGTH ( I ) f D PS I X, 0 PH ITH, 

*V(I)»X 

CONTINUE 


N1«N-1 

DO 3 I«2#N1 

DYl-Y(I)-Yn-l) 

DY2«Y(I+1)-Y(I ) 

VP{ I)«(DY1*V(I+1)/0Y2-DY2*V( I-l ) /D Y1 ) / (OYl* DY2 )-V ( I ) ♦( 0Y1-DY2 ) / ( DY 
11*DY2 ) 

VPP (I )«(DY1*V( I+1)+DY2*V(I-1 )-V(I ) ♦ (0Y1+0Y2 ) ) / ( . 5*D Y1*D Y2*( D YI+DY2 
1 ) ) 

CONTINUE 

VP{N)«(V(N)-V{N-1))/DY2 
VPP(N)«(V(N-1)-V(N) ) /DY2**2 
DY1»Y(2)-Y(1) 

DY2»Y(3)-Y(1) 

VP(1)«(DY2*V(2 )/DYl“DYl*V( 3) /DY2 ) / ( DY2-DY1 )-V ( 1 ) * ( DY1 + DY2 )/ (DY1*DY 


P 1 
P 2 
P 3 
P A 
P 5 
P 6 
P 7 
P 8 
P 0 
P 10 
P 11 
P 12 
P 13 
P lA 
P 15 
P 16 
P 17 
P 18 
P 19 
P 20 
P 21 
P 22 
P 23 
P 2A 
P 25 
P 26 
P 27 
P 28 
P 29 
P 30 
P 31 
P 32 
P 33 
P 3A 
P 35 
P 36 
P 37 
P 38 
P 39 
P AO 
P A1 
P A2 
P A3 
P AA 
f A5 
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♦ TIDY* 


SUBROUTINE VVEL 


C 

C700 

c 


12) 

P 

^6 

X1»Y (2) 

P 

47 

X2» Y( 3) 

P 

48 

X3«Y(<t) 

P 

49 

XA-Y (5 ) 

P 

50 

B2«X2*X3+X2*X4 + X3*X<> 

P 

51 

DL1«X1*( X1-X2)*{X1-X3 )*(X1-X<^) 

P 

52 

C2»X1*X3+X1*X<> + X3*XA 

P 

53 

0L2-X2’M X2-XD+ (X2-X3)*(X2-X^) 

P 

54 

02«Xl*X2 + XltX4+ X2*XA 

P 

55 

OL3»X3*( X3-X1 )*(X3-X2 )*(X3-X4 ) 

P 

56 

E2-X1*X2+X1*X3+X2*X3 

P 

57 

DL^»X<**(X4-XD* (X^-X2 )*(X4-X3 ) 

P 

58 

VPP(l)«2.*B2*V(2)/DLl + 2.*C2t'V( 3) /DL 2+2 .*D2 * V( 4 ) / DL 3+ 2. ♦ E2*V (5 ) /DL<t 

P 

59 


P 

60 

VPP (1 )«VPP( 2) 

P 

61 

VPP(N)«VPP{N-1) 

P 

62 

RETURN 

P 

63 

FORMATdH >1X#8E13.A) 

P 

64 


P 

65 

END 

P 
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